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SECTION  1 


OVERVIEW 


This  is  the  final  report  of  a theoretical  and  experimental  investigation  into 
some  properties  of  holograms.  Optical  holograms  suffer  from  the  existence  of 
complex  conjugate  terms  and  nonlinearities  which  add  an  unavoidable  amount  of 
clutter.  Because  we  were  able  to  perform  our  experiments  using  a digital  com- 
puter (Data  General  Eclipse),  these  conjugate  terms  and  nonlinearities  were 
not  present.  This  provided  a certain  cleanness  in  the  mathematics  which  will 
be  evident.  It  also  showed  that  the  essential  properties  of  the  hologram  are 
not  restricted  to  a wave  propagation  environment.  Thus,  we  have  chosen  to 
provide  the  general  name  of  Machine  Holography  for  our  investigation.  The 
name  implies  that  the  principles  investigated  here  may  be  embodied  by  diverse 
types  of  machinery,  not  just  laser  based  optical  systems. 

The  investigation  was  divided  into  two  packages.  The  first  explored,  the  ulti- 
mate storage  capacity  of  what  we  have  come  to  call  Conventionally  Multiplexed 
Holograms.  These  are  holograms  which  are  multiplexed  by  a summation  process 
whereby  each  reference  pattern  participating  is  uncorrelated  with  the  other 
references  in  the  ensemble.'  This  subject  was  partially  explored  by  La  Macchia 
and  White*  for  certain  types  of  optical  holograms,  and  has  been  generalized 
here,  while  also  treating  important  new  cases. 

Theoretical  treatment  shows  that  the  amount  of  information  which  may  be  stored 
in  this  type  of  hologram  is  proportional  to  the  number  of  degrees  of  freedom 
available  in  the  output  of  the  holographic  system.  The  theory  also  predicts 
that  it  does  not  matter  how  the  information  is  divided  for  storage.  There  may 
be  only  a few  pages  of  data  with  many  bits  on  each  page;  conversely,  a few 
bits  may  be  stored  on  each  of  many  pages.  Only  the  total  number  of  bits  to  be 
stored  enters  the  calculation.  These  theoretical  results  were  confirmed  by 
means  of  computer  simulations. 

The  primary  goal  of  the  second  package  of  work  was  to  demonstrate  that  a holo- 
gram could  be  trained  to  reliably  distinguish  between  two  highly  correlated 
(reference)  patterns. ✓ This  effort  broke  entirely  new  ground.  In  an  effort  to 
design  a training  paradigm  to  ^demonstrate  this  prediction,  a new  kind  of  machine 
was  discovered.  This  new  machine  put  a negative  feedback  loop  around  the  holo- 
gram. The  hologram  then  became  part  of  a servo  mechanism  in  which  both  the 
controller  of  the  servo  and  the  object  of  control  was  the  stored  hologram, 

Theory  predicted  that  this  new  type  of  adaptive  hologram  is  capable  of  cleaning 
up  clutter  and  false  reconstructions.  It  also  indicated  that  this  type  of 
hologram  could  be  trained  for  differential  pattern  recognition  by  a remarkably 
simple  procedure.  Experiments  have  confirmed  both  of  these  predictions.  We 
report,  then,  the  pioneering  of  new  and  potentially  very  fertile  territory  for 
scientific  exploration  and  technological  exploitation. 


*La  Macchia,  J.  T. , and  White,  D.  L.  "Coded  Multiple  Exposure  Holograms," 
Applied  Optics,  Vol.  7,  No.  1,  (Jan  1968)  pp.  91-94. 
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MULTIPLEXING  THEORY  AND  EXPERIMENT 


2.1  HOLOGRAPHIC  RECORDING  AND  RECONSTRUCTION 

2.1.1  INTRODUCTION 

The  formation  and  reconstruction  of  a hologram  can  be  regarded  as  an  exercise 
in  encoding  and  decoding  an  image.  The  encoding  activity  involves  the  linear 
spatial  transformation  of  the  image  and  an  appropriately  configured  reference 
pattern.  The  transformed  reference  acts  as  an  encoding  key  and  is  mixed  with 
the  transformed  image.  The  recorded  result  is  the  hologram. 

Reconstruction  involves  inserting  the  reference  into  the  system,  transforming 
it  and  playing  it  against  the  hologram.  The  reference  then  decodes  th'*.  infor- 
mation contained  in  the  hologram.  After  an  appropriate  inverse  transformation, 
the  image  is  restored  to  something  like  its  original  form.  The  degree  of 
similarity  between  the  original  image  and  its  reconstruction  will  be  deter- 
mined by  the  properties  of  the  reference  pattern  used,  by  the  nature  of  the 
linear  transformation,  and  by  the  mode  of  multiplexing  which  may  be  employed. 

2.1.2  DERIVATION 

Let  us  consider  first  a simple  unmultiplexed  hologram.  For  the  following 
derivations  we  assume,  for  clarity,  continuum  images  and  transformations. 
Perhaps  the  most  typical  and  useful  linear  transformation  is  the  Fourier 
transform,  which  is  adopted  here.  All  functions  will  be  annotated  with  a 
single  independent  variable  to  avoid  undue  notational  complexity.  Generaliza- 
tion to  two  parameters  for  a full  description  of  two-dimensional  images  is  a 
straightforward  but  tedious  exercise. 

Suppose  we  have  a subject  image  s(x)  which  we  wish  to  record  and  reconstruct. 
This  will  be  paired  with  a reference  image  r(x) , which  will  be  used  to  encode 
and  decode  the  hologram.  These  two  patterns  are  connected  with  the  hologram 
domain  through  a Fourier  transform 


R(0  = . 

/ r(x) 

-i£x 

e dx 

(1) 

S(C)  = . 

f s(x) 

-i£x 

e dx 

(2) 

During  the  recording  process  both  the  reference  and  subject  images  are  trans- 
formed according  to  Equations  1 and  2.  Before  they  are  mixed  together  the 
reference  undergoes  a further  transformation  into  its  complex  conjugate  form. 
The  hologram  is  recorded  as  the  product  of  the  transformed  subject  and  complex 
conjugate  transformed  reference: 

H(U  = S(C)  R*(C)  (3) 
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Reconstruction  involves  the  reinsertion  of  the  reference  pattern,  r(x)  , into 
the  reference  channel.  The  reference  is  Fourier  transformed  and  delivered  to 
the  hologram  domain.  During  the  reconstruction  the  transformed  reference 
does  not  undergo  complex  conjugation.  Instead,  it  simply  multiplies  the 
hologram  resulting  in  a response  pattern: 

A(U  - S(U  R*(t)  RCO  (4) 

Finally,  this  response  undergoes  an  inverse  Fourier  transform  to  produce  the 
reconstructed  image 

a(x)  - F_  1 1 ACO  ] (5) 

The  degree  of  resemblance  between  the  reconstruction  a(x)  and  the  original 
image  s(x)  will  depend  on  the  properties  of  the  reference  r(x)  under  a Fourier 
transformation.  This  can  be  seen  by  noting  that  the  inverse  Fourier  transform 
of  the  triple  product  of  Equation  4 results  in  the  triple  convolution: 

a(x)  - s(x)  * (r(x)©r(x)j  (t>) 

Where  r(x)©r(x)  represents  an  autocorrelation  of  the  reference  pattern. 

This  autocorre lat ion  serves  the  purpose  of  a system  point  spread  function. 

The  system  output  can  therefore  be  written 

a(x)  » s(x)  * p(x)  (7) 

where 

p (x)  - r(x)  © r (x) 

It  is  within  this  context  that  the  significance  of  the  proper  choice  of 
reference  can  be  understood.  Suppose  that  the  reconstruction  point  spread 
function  is  delta  f unction- like;  then 

p ( x)  - S(x)  (8) 

but  by  definition 

a(x)  * 6(x)  * s(x)  - s(x) 

Thus,  a delta-like  reference  autocorrelation  will  produce  a perfect  recon- 
struction. In  general,  we  might  expect  a reference  to  be  delta-like  only  in 
approximation.  In  such  a case,  the  reference  acts  as  a filter  which  degrades 
the  reconstruction  somewhat.  For  example,  if  the  reference  autocorrelation  is 
finite  in  width  it  will  act  like  a low  pass  filter  to  smooth  or  blur  the 
reconstructed  imago. 
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Another  type  of  reference  pattern,  which  Is  of  particular  Interest  in  this 
investigation,  is  a rami  on  function.  The  autocorrelation  of  such  a function 
will  usually  show  a very  sharp  ami  narrow  spike  at  autoalignment. 


I 


This  spike  is  delta-like  and  leads  to  a high  resolution  recons t ruct ion  of  the 
subject  image.  When  the  reference  pattern  is  out  of  »utoal lgiunent , the  auto- 
correlation will  produce  a low  level  random  "noise."  The  convolution  of  this 
noise-like  pattern  with  the  subject  will  result  in  a scrambled  hare  of  clutter. 
Thus,  holographically  associating  a subject  image  with  a random  reference 
pattern  leads  to  a sharp  reconstructed  image  superimposed  on  a low  level  random 
clutter  background. 

Consider  next  what  happens  when  the  reference  used  for  reconstruction,  r'tx), 
is  not  identical  tv’  the  reference  used  for  recording  the  hologram.  The 
reconstruct  ion  has,  then  the  form 

alx)  - six)  * | rix)  © r'  lx)j  f^) 

The  system  point  spread  function 

p(x)  - rfx)  © r'ix)  (10) 


will  be  degraded  bv  the  degree  of  mismatch  between  r (x)  and  r'ix).  The  cor- 
relation coefficient,  taken  between  rfx)  and  r'(x),  provides  a measure  of  the 
similarity  between  these  two  functions  and,  therefore,  of  the  quality  of  the 
reconstruction.  If  r(x)  and  r'ix)  are  completely  uncorrelated  the  point  func- 
tion of  Equation  hi  will  usually  produce  a random  pattern  of  low  amplitude. 
This  random  function,  in  turn,  means  that  the  attempted  reconstruct  ion  o t 
Equation  *>  will  contain  onlv  low  level  random  noise. 

This  is  a useful  property  which  can  be  put  to  use  in  multiplexing  holograms. 
Suppose  we  establish  an  ensemble  ot  reference  patterns:  |l'k^x^}  • The  con- 
dition to  be  met  bv  these  patterns  is  that  each  member  is  completely 
uncorrelated  with  the  other  members  of  the  ensemble.  Thus 

r^lx)  © fj  lx)  - Six)  + lx)  (ID 


where 


" 0,  when  k f l. 

5.  , ■ l , when  k ■ L 

kL 

lx)  ■ a random  function  dependent  on  k and  1 


U 

I 
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If  the  set  of  references  has  this  property,  multiplexing  a set  of  subject 
images  can  be  achieved  by  addition.  Each  subject  is  paired  with  one  reference 
drawn  from  the  ensemble  of  references.  The  hologram  is  formed  by  adding  up 
all  of  these  pairings.  Thus,  the  form  of  the  multiplexed  hologram  is: 

H(f.)  - SL(C)  R*(0  + S2(U  R*(t.)  + S3  R*  + ...  (12) 

or 

H(0  - rk  Sk<£)  R^(f) 

If  we  write  this  in  the  output  image  domain  we  have 


h(x)  « sk(x)©rk(x) 


(13) 


Reconstruct  ion  of  one  particular  subject  will 
associated  reference  into  the  system.  In  the 
be  written  as: 


aL(x) 
a^  (x) 


rLU)  rkU)  * \(X)) 

2^k  (rL(x)©rk(x))  * s^fx) 


occur  when  we  reinsert  the 
image  domain  this  process  may 


(14) 

(15) 


Making  use  of  Equation  11  in  Equation  15  we  have: 


a (x)  * 5(x)  * s (x)  + Z P.  . (x)  * s (x) 
L L k kL  k 


a^(x)  - Sj (x)  + noise 


(16) 


The  reconstruction  of  the  subject  associated  with  the  chosen  reference  is 
accompanied  by  clutter  noise  generated  by  various  random  cross  terms.  As  we 
increase  the  number  of  holograms  stored  in  a multiplexed  ensemble  we  will  also 
get  a corresponding  increase  in  noise.  This  noise  ultimately  sets  a limit  on 
the  amount  of  information  which  may  be  stored.  This  limit  will  be  calculated 
in  detail  in  Paragraph  2.2.  Note  the  important  point  that  multiplexing  by 
this  means  (called  conventional  multiplexing)  depends  on  a correlation  process. 
This  correlation  is  the  consequence  of  recording  the  hologram  in  the  form  of  a 
sum  of  Fourier  transforms  of  the  contributing  Images. 


In  summary,  the  process  of  holography  depends  on  a convolutional  association 
between  a reference  and  a subject  pattern.  Reconstruction  results  from  the 
creation  of  a synthetic  point  spread  function  resulting  from  the  auto-  or 
crosscorrelation  of  two  references.  Multiplexing  is  possible  because  the  use 
of  an  uncorrelated  reference  in  the  reconstruction  results  merely  in  the 
generation  of  a haze  of  clutter  noise. 
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2.2  SIGNAL-TO-NOISE  RATIO  FOR  RECONSTRUCTION  OF  A MULTIPLEXED  HOLOGRAM 

2.2.1  BACKGROUND 

One  of  the  major  tasks  of  the  contract  was  to  explore  the  storage  capacity  of 
a conventionally  multiplexed  hologram.  A prime  goal  was  to  find  out  which 
factors  were  significant  in  determining  the  storage  limit.  We  also  wished  to 
know  what  the  functional  dependence  was  on  these  significant  parameters. 

A conventionally  multiplexed  hologram  is  defined  as  a multiply  exposed  hologram 
formed  by  adding  together  separately  recorded  subholograms.  The  result  is  a 
summed  ensemble.  Random  encoding  of  the  individual  references  paired  with  each 
subject  is  used  to  make  each  subject  in  the  ensemble  selectively  retrievable. 

A theoretical  investigation  of  this  task  had  previously  been  carried  out  for 
a type  of  optical  hologram*.  But  this  investigation  did  not  consider  the 
possibility  that  the  subjects  would  be  made  up  of  multiple  data  points.  Some 
questions  also  existed  as  to  the  effect  of  sampling  on  the  storage  capacity. 

These  issues  have  now  been  settled  with  the  following  principal  results: 

(a)  The  storage  capacity  of  a random  reference  multiplexed 
hologram  depends  directly  on  the  number  of  degrees  of 
freedom  of  the  image  output  of  the  system. 

(b)  The  storage  capacity  varies  inversely  with  the  energy 
signal-to-noise  ratio  (clutter)  required  at  the  system 
output . 

(c)  It  does  not  matter  in  what  configuration  the  data  is 

stored:  a few  frames  of  data  with  many  data  points  will 

produce  the  same  signal-to-clutter  ratio  as  many  frames 
with  few  data  points  each.  Only  the  total  number  of  data 
points  extracted  from  all  the  frames  matters. 

The  last  result  appears  to  be  new  and  is  important  in  assessing  the  performance 
of  holographic  data  storage  systems. 

2.2.2  DERIVATION 

Let  us  assume  that  the  hologram  plane  is  composed  of  n sample  points.  We 
wish  to  store  a total  of  H subholograms  within  this  domain.  Each  subhologram 
is  presumed  to  have  been  formed  with  a unique  random  reference  pattern,  each 
having  a narrow  autocorrelation  function.  We  suppose  that  at  each  sample  point 
in  the  hologram  the  contributed  amplitude  from  any  given  subhologram  in  the 
ensemble  has,  on  the  average,  a value  of  a. 

*La  Macchia,  J.  T. , and  White,  D.  L.  "Coded  Multiple  Exposure  Holograms," 
Applied  Optics,  Vol.  7,  No.  1,  (Jan  1968)  pp.  91-94. 
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We  will  find  that  the  reconstruction  Includes  clutter  terms  from  two  sources: 

o Cross  terms  within  the  proper  reference  function. 

o Cross  terms  between  uncorrelated  reference  functions. 

The  limits  of  data  storage  capacity  will  depend  on  the  amount  of  clutter  which 
creeps  Into  the  reconstruction  process  from  these  two  sources. 

The  proper  perspective  for  performing  signal  magnitude  calculations  is 
developed  by  considering  the  following  points: 

(a)  Every  sampled  Image  may  be  treated  as  a vector  array  of 
sample  values. 

(b)  The  vector  has  a length  derived  from  an  appropriate  sum  of 
these  values. 


(c)  A measure  of  the  magnitude  of  the  image  vector  is  the  square 
of  its  length.  This  measure  corresponds  to  the  intensity  of 
a conventional  optical  wave. 

(d)  If  we  examine  the  image  in  the  output  domain  we  find  that 
each  of  its  information  hearing  points  is  Isolated  and 
unaffected  by  its  associates.  Titus,  if  the  average  amplitude 
of  each  of  these  points  is  o,  the  square  magnitude  for,  p, 
points  will  Just  be  the  sum  of  the  squares  of  the  individual 
points : 

2 

* p o' 

Another  way  of  saying  this  is  that  the  magnitude  of  the 
output  image  vector  is  determined  hv  a root  sum  square  (RSS) 
addition  of  the  output  image  points. 

(e)  If  we  examine  the  discrete  Fourier  transform  (l)FT)  of  this 
image,  we  find  that  the  sample  values  in  this  new  space  are 
no  longer  Independent  and  do  not  combine  in  RSS  fashion. 

In  fact,  thev  are  coherently  linked.  This  means  that  their 
amplitude  values  add  directly,  and  the  square  magnitude  is 
the  square  of  the  sum  of  the  individual  amplitudes. 


Thus,  if  the  average  amplitude  of  a transformed  image  sample 
in  the  Fourier  domain  is  a,  and  if  there  are,  n,  such  samples, 
then  the  square  magnitude  of  the  image  will  he 
) 

Mj  - (no) ‘ 


Combining  Equations  18  and  19  we  find  that: 


The  effect  of  clutter  may  be  calculated  by  a detailed  examination  of  the 
reference  function  during  reconstruction.  We  recognize  that  the  reference  may 
be  considered  as  a sample  data  function,  in  the  Fourier  domain,  of  the  form: 

reference  - R,  (£) 
k 

Where  £ is  an  index  which  locates  a particular  sample  in  a space  of  n such 
samples.  The  index  k labels  the  particular  reference  function  from  the 
ensemble,  which  is  being  used. 

In  a conventionally  multiplexed  hologram  the  resulting  reference  involves  a 
linear  summation  of  many  terms.  If  we  ignore  the  associated  subject  functions 
we  may  write  a reference  ensemble  as: 

jR*  “ ^k  Rk(f')  (20) 

Let  us  now  go  through  the  reconstruction  process  by  inserting  a particular 
reference  Rj(£).  The  result  will  be  to  output  a function  of  the  form 

P<£)  - Rj(£)  \ (21) 

p(0  - |r.(0 I 2 + r.  r (o  r*(o 

k»*J  J k 

Let  us  next  express  these  reference  components  in  terms  of  their  image  domain 
representations.  The  connection  is  a discrete  Fourier  transform. 


P(0  ■ K r.U)  e~lU|2 


. (22) 


Let  us  consider  t irst  the  second  term  on  the  right,  which  contributes  pure 
clutter.  Rearranging  the  summations  in  this  term  we  get: 


C1  “ 1 Ex  ri(x)  rk(y)  °~1' (X_y> 

kf»J  x V -•  K 

This  is  a random  summation  evidenced  by  the  fact  that  the  factor  e~*s^X  is 
nonzero  and  ranges  uniformly  over  all  values  from  -1  to  +1  and  from  -i  to  +t . 
Thus  the  terms  must  be  summed  together  RSS  to  get  the  equivalent  vector 
magnitude.  Let  us  assume  first  that  the  average  amplitude  Independent  of 
phase  effects  will  be  a such  that 


a - <r  (x)  r ( v) > 

.1  k 
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Next,  It  Is  necessary  to  count  terms.  Wo  note  that  both  x and  y range  over  n 
samples.  So  the  double  sum  over  these  parameters  contributes  n-  terms.  1£ 
the  ensemble  has  H separate  subholograms,  H-l  of  these  will  contribute  to  Cj 
and  the  remaining  one  will  result  In  the  reconstruction  of  the  Image.  We 
therefore  have  a total  of  (H-l)n-'  amplitudes  contributing  to  clutter  from  this 
component  of  the  reconstruction.  Since  these  terms  are  randomly  phased,  their 
magnitude  will  be  the  result  of  an  KSS  combination.  The  clutter  square 
magnitude  from  this  component  will  therefore  be 

2 2 2 

C.  - (H-l)n  a (25) 

We  examine  next  the  first  term  on  the  right  In  Equation  22: 

U r lx)  e"Ux|2  - E |r  ,(x)|2  + Z I’.rdx)  r.(y)  e~  lf'(x_y)  (2b) 

x .1  x .1  Xfly  j J 

The  first  term  on  the  right  Is  the  reconstruction  term.  There  are  n of  these 
terms  each  represented  by  an  average  amplitude  it.  Note  that  there  is  no  phase 
variation  so  that  the  summation  of  amplitudes  (over  the  Index  x)  Is  coherent. 
This  leads  to  the  square  magnitude  result  given  by  Equation  18. 

The  second  term  on  the  right  Involves  an  incoherent  summation  and  contributes 
to  the  clutter.  (Note,  that  for  certain  types  of  reference  functions,  such 
as  a delta  function  In  the  image  plane,  this  term  Is  zero  and  results  in  a 
perfect  reconstruction  when  there  Is  no  multiplexing.)  We  recognize  that 
there  must  hi'  n--n  contributions  of  amplitude  it  from  this  term.  Since  the 
summation  of  these  contributions  is  Incoherent  they  will  be  KSS  added  and  will 
therefore  combine  in  square  magnitude  with  Equation  25.  The  total  clutter  is 
therefore 

2 ■>  y ■>  •> 

C • (H-l) n it*"  + (n*’-n)u*' 

2 > 2 2 

C ■ Hn  it*"  - na 

C2  - (Hn-l)nit2  (27) 


This  Is  the  clutter  as  measured  In  the  Kourler  domain.  Its  Impact,  however, 
will  be  on  the  reconstruct  Ion  In  the  image  domain.  Suppose  that  the  Image 
domain  contains  a total  of  m sample  points  or  degrees  ot  freedom.  It  the 
average  amplitude  ot  the  clutter  in  the  image  domain  Is  V,  then  the  square 
magnitude  of  this  clutter  will  be  (assuming  incoherent  summation) 


By  equating  Equations  27  and  28  we  determine  the  average  clutter  at  each 
sample  point  of  the  output  domain: 


n(Hn-l)a^  = rav^ 


P 


v * /—  (Hn-1)  a 


(29) 

(30) 


(Note:  The  assumptions  leading  to  Equations  29  and  19  involved  implicit  use 
of  Parseval’s  theorem.) 

Now  we  can  estimate  the  mean  signal-to-clutter  ratio  by  taking  the  ratios  of 
Equations  19  and  30. 


not 


7 


V / pn  (Hn-1)  a 
m 


(31) 


H 


/ 0 ) „ / s \ 

f nm 

\v/  ~ \ N j A J 

P(Hn-l) 

(32) 


represents  an  amplitude  signal-to-noise  ratio  where  the  noise  is  due 
A solely  to  crosstalk  clutter.  Rearranging  Equation  32  we  get 


m 

. P 

/ S\  2 

n 

lN)  A 

(33) 


Equation  33  is  the  primary  result  of  interest. 

The  total  amount  of  data  stored  in  the  type  of  multiplexed  hologram  is  P 
points  stored  in  each  of  H holograms  or  PH.  We  define  the  data  total  as: 

PH  = D ( 3A) 

Note  that  for  a sufficient  ratio  of  n to  p we  have,  approximately: 

(35) 


_ m 
D » 


(") 


T 

A 
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Thus,  the  total  amount  of  data  stored  is  in  proportion  to  the  number  of  samples 
(or  degrees  of  freedom)  in  the  OUTPUT  image  field.  This  is  a new  and  very 
interesting  result.  An  Important  consequence  of  Equations  34  and  3S  is  that 
it  does  not  matter  how  the  division  is  made  between  numbers  of  data  pages  (H) 
and  numbers  of  data  points  (P>  per  page. 


The  inverse  square  relationship  with  the  amplitude  signal-to-clutter  ratio 
corresponds  to  an  inverse  relationship  to  the  "energy"  (intensity) 
signal-to-clutter  ratio  familiar  from  coherent  optics.  In  the  special  case 
where  the  output  image  plane  has  the  same  number  of  degrees  of  freedom  as  the 
Fourier  or  hologram  plane,  we  may  write: 


(3b) 


(37) 


An  example  calculation  is  useful  in  establishing  a baseline  for  the  experi- 
ments to  be  described  shortly.  Suppose  we  have  an  image  field  of  409b  points 
and  that  our  data  field  has  many  fewer  points,  so  that 


P « n 

then 

I s \ = m >/409b 

\ N / A >TpH  J~PH” 


M 

. b4 

lNJ 

A -J  PH 

(38) 

(39) 

(40) 


Equation  40  describes  the  experimental  conditions  and  results  which  we  wish 
to  verify. 
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2.3  RESULTS  OF  MULTI PLEXINCI  EXPERIMENTS 


In  order  to  validate  the  theoretical  arguments  presented  in  the  previous 
paragraphs,  a series  of  experiments  was  performed.  The  experiments  involved 
the  formation  and  reconstruction  of  multiplexed  holograms  using  random 
reference  patterns  and  conventionalized  subject  Images. 

The  primary  object  of  these  experiments  was  to  determine  the  dependence  of 
s ignal-to-c lutter  ratio  on  the  amount  of  data  stored.  Also,  we  wished  to 
find  out  whether  or  not  the  data  limit  was  independent  of  page  formatting, 
as  predicted. 

In  a typical  experiment,  several  reference  subject  pairs  were  associated  and 
added  together  to  form  a conventionally  multiplexed  hologram.  Then  each  of 
the  reference  images  was  submitted,  in  turn,  to  the  hologram  and  the  associ- 
ated subject  image  was  reconstructed. 

The  signal-to-clutter  ratio  was  determined  by  averaging  the  amplitudes,  of  the 
data  points  and  comparing  this  value  with  an  average  of  samples  taken  from  the 
data-tree  border  region.  By  varying  the  number  of  data  points  in  the  image, 
the  number  of  samples  In  the  field,  and  the  number  of  subholograms  in  an 
ensemble  a parametric  study  could  be  performed. 

The  apparatus  for  this  experiment  was  a Data  General  Ec’ipse  computer,  com- 
bined with  a Printronix  printer.  The  print  characters  v ore  under  software 
control  and  were  programed  to  produce  Images  with  ten  levels  of  grayscale. 

The  image  format  for  the  experiments  presented  here  consisted  of  a square 
array  of  409b  samples.  The  Fourier  transform  plane  also  contained  409b  samples. 
The  data  arrays  were  also  chosen  as  square  arrays  ranging  from  one  point  to 
sixteen . 

In  Figures  1 and  2,  sample  outputs  from  five  experiments  are  shown.  The  data 
image  for  these  experiments  consisted  of  a square  array  of  four  points.  In 
order  to  show  that  proper  independent  reconstruction  was  taking  place  the  data 
arrays  used  to  form  different  subholograms  were  located  in  different  parts  of 
the  scene. 

A cursory  examination  of  Figures  1 and  2 shows  the  qualitative  feature  that 
the  signal-to-clutter  ratio  (S/N)  decreases  as  the  number  of  images  multiplexed 
in  an  ensemble  increases.  This  is  as  expected.  The  quantitative  results  con- 
formed to  the  theory,  also.  This  is  illustrated  by  the  numbers  of  Table  l, 
and  the  plot  of  these  numbers  given  in  Figure  4.  The  theoretical  values  given 
in  the  second  column  are  computed  using  Equations  34  and  35.  The  standard 
deviation  given  in  the  last  column  is  compiled  from  measurements  of  the 
reconstruction  signal-to-clutter  ratios  of  each  member  of  a given  ensemble. 


A.  Nonmultiplexed,  one  image  stored  B.  Multiplexed,  two  images  stored 
S/N  = 33.56  S/N  = 21.54 


• Sample  field  has  4096  points 

* . * : *•  .*  . . * : 

• Data  Image  has  4 points 

• S/N  ratios  are  for  particular 

• ‘ *' 

images  being  reconstructed 

**•*••.  t • : • • 

• Each  image  represents  a 

■’  •••• 

different  experiment 

j 

C.  Multiplexed,  four  images  stored 
S/N  - 15.71 

93-S-1? 


FIGURE  1.  TYPICAL  RECONSTRUCTIONS  OF  MULTIPLEXED  HOLOGRAM  EXPERIMENTS 


Multiplexed,  eight  images  stored 
S/N  - 10.97 


B.  Multiplexed,  sixteen  images  stored 
S/N  - 8.023 


FIGURE  2.  ADDITIONAL  SAMPLE  MULTIPLEXING  EXPERIMENTS  WITH  FOUR  POINT 
DATA  IMAGES  (SUBJECT) 


Nonmult ip lexed,  one  image 


B.  Multiplexed,  two  imag 


Multiplexed,  four  images 
S/N  = 8.22 


• Sixteen-point  data  image  used 

• Each  image  is  a representative  reconstruct  ion  from  one  experiment 


FIGURE  3.  SECOND  SERIES  OF  MULTIPLEXING  EXPERIMENTS 


I 

l 
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16-POINT  IMAGE 


It  should  also  he  noted  that  the  use  of  a random  reference  results  In  some  net 
clutter  even  when  only  a single  image  is  being  recorded  and  reconstructed . 

This  is  Illustrated  by  Figure  1,  view  A. 


This  series  of  experiments  was  repeated  for  the  case  where  the  data  image  con- 
sisted of  an  array  of  sixteen  points  (Figure  J).  Again,  as  demonstrated  in 
Table  1 and  Figure  4 the  experiments  and  the  theory  are  in  strict  quantitative 
agreement . 

TABLE  1.  SIGNAL- TO-CLUTTER  RATIOS  FOR  MULTI  FLEXED  HOLOGRAMS 

Output  image  domain  has  409b  samples.  Number  of  data  points  in  each  image  is 
as  listed. 


No.  of  Multiplexed 
Holograms  (H) 

S/N  Theoretical 

S/N  Experimental 

0 

n 

4 Data  Point  Image 

1 

32 

33.56 

- 

> 

22. b3 

22.3b 

0.825 

4 

lb 

15.  bb 

0.91b 

8 

11.31 

11.45 

0.31  1 

lb 

8 

8 . bb 

0.  bb2 

lb  Data  Point  Image 

l 

lb 

1 4 . 9(> 

- 

) 

11.31 

12.27 

0.89 

4 

8 

8 . b 7 b 

0 . S 7*JQS> 

8 

S.bS 

5 . b 7 4 

0. 19b 

A final  experiment  in  this  class  was  performed  to  test  the  assertion  that  the 
data  which  could  be  stored  is  primarily  dependent  on  the  number  of  degrees  of 
freedom  of  the  output  image  plane,  in  accord  with  Equation  IS.  In  this  experi- 
ment the  random  reference  employed  was  windowed  so  as  to  be  a minor  portion 
of  the  field  of  view. 

The  consequence  of  this  windowing  is  that  during  reconstruction  the  clutter  is 
weighted  bv  the  autocorrelation  of  the  window  function.  Moreover,  because  the 
clutter  is  confined  in  the  field  of  view,  its  magnitude  is  enhanced.  It  is 
this  enhancing  factor  which  reduces  the  signal-to-clut ter  ratio  and, 
consequently,  the  amount  of  information  which  can  be  stored. 


The  reconstruction  of  this  experiment  is  illustrated  in  Figure  S.  The 
weighting  bv  the  autocorrelation  of  t He  window  function  is  clearly  evident. 
The  enhancement  of  clutter  is  evident  when  this  illustration  is  compared 
with  Figure  3,  view  A. 


k n .*  ■ • 


V • S 9 » 

• •:*  • P| 


• Hologram  is  unmultiplexed 

• Flutter  field  tapered  according  to  auto  correlation  ot  reference  window  function 

• Peak  clutter  very  much  higher  than  if  it  were  more  widely  distributed 
(compare  with  Figure  J,  View  A 

• No  attempt  made  to  measure  signal-to-clut ter  ratio  in  this  experiment 


FllJURF.  S. 


SAMPl.F.  RECONSTRUCT l ON  WITH  ORUUNAI  RFFF.RFNt'F  CONSTRAINED  TO  A FRACTION 


OF  THE  FI  FI  D OF  VIEW 


SECTION  3 


ADAPTIVE  HOLOGRAMS 


3 . 1 THEORY 

3.1.1  INTRODUCTION 

It  was  asserted  by  the  proposal  that  a hologram  could  be  trained  to  distinguish 
between  two  very  similar  reference  patterns. 

The  notion  was  that  the  common  subset  of  the  two  reference  patterns  could  be 
canceled,  leaving  distinguishing  features  which  could  be  used  to  reconstruct 
identifying  subject  images.  The  technique  makes  use  of  the  ability  of  a 
hologram  to  coherently  subtract  by  means  of  a phase  change.  The  isolation  of 
the  characteristic  features  was  considered  to  be  analogous  to  the  Gram-Schmidt 
orthogonalization  technique.  This  concept  is  explored  in  detail  in  Appendix 
A,  which  is  borrowed  from  the  proposal. 

During  a design  of  the  training  algorithm,  which  would  implement  this  concept, 
it  was  realized  that  a negative  feedback  process  would  achieve  the  desired 
goal.  The  idea  was  to  wrap  a feedback  link  around  the  outside  of  the  holo- 
graphic system  so  that  the  reconstructed  image  would  be  subtracted  from  the 
current  subject  input.  The  holographic  recording  mechanism  would  use  this 
difference  image,  combine  it  with  the  current  reference,  and  add  the  result- 
ing subhologram  to  the  preexisting  ensemble  according  to  the  rules  of  conven- 
tional multiplexing.  Figure  6 provides  a schematic  representation  of  this 
process. 

By  supplying  two  different  paired  subject  and  reference  patterns  to  the  system 
inputs  according  to  some  sort  of  alternation  scheme,  it  was  felt  that  the 
common  subsets  of  the  reference  patterns  would  be  subtracted  off.  This  would 
leave  the  reconstructions  to  the  mutually  othorgonal  features  of  the  two 
reference  patterns. 

This  configuration  was  extremely  interesting,  but  it  involved  a great  deal  of 
computation  in  performing  the  forward  and  inverse  Fourier  transformations  for 
each  iteration. 

After  some  study  it  was  realized  that,  because  of  the  linearity  properties  of 
these  transforms,  the  feedback  loop  could  be  brought  inside  the  transforms. 

The  adjustment  process  then  would  be  entirely  in  the  Fourier  domain.  Except 
for  initial  transformations  to  prepare  the  subject  and  reference  and  an 
inversion  to  display  the  ultimate  result,  there  would  be  no  need  for  repetitive 
Fourier  transformations. 
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Tills  change  led  to  an  adaptive  holographic  system  which  is  schematically 
Illustrated  bv  Figure  7.  This  system  is  easy  t o analyze  and  has  a number  of 
exceptionally  interesting  properties.  Implementation  of  the  arrangement  in 
software  showed  that  it  can  be  unstable.  Some  minor  adjustments  in  the 
description  of  the  iteration  process  insured  unconditional  stability  with- 
out violating  the  clarity  of  the  basic  idea.  These  matters  will  be  discussed 
in  detail  in  the  next  section. 


The  possibility  ot  this  simple  teedback  arrangement,  with  its  powerful  properties, 
was  a surprise.  it  is  the  outstanding  result  of  this  study  effort. 

i.i..’  hash:  thkoky 


In  the  adaptive  holographic  scheme  the  set  of  values,  Hn(f.) , of  the  hologram, 
at  any  given  time,  is  derived  from  a summation  of  past  inputs.  The  process  is 
related  to  convent  Iona  1 holographic  multiplexing  because  of  the  formation  of 
an  ensemble  ot  subholograms.  At  each  step  of  the  process  the  current  patterns 
arriving  through  the  subject  and  reference  channels  are  multiplied  together. 
This  product  forms  a now  subhologram  which  is  added  to  the  preexisting  ensemble 
ot  holograms.  I'he  equation  which  describes  this  summation  is 


II  (f.) 
n 


*Bn(° 


* 

K (f.) 
n 


+ H . (O 

n-  l 


(Ml 


In  this  equation  H.._  i (£)  is  the  preexisting  multiplexed  hologram.  is 


the  new  multiplexed  hologram.  Rn(Q  is 
is  the  current  subject  pattern  and  g is 

The  subject  pattern  is  derived  from  two 
version  ot  the  current  subject  input 


the  current  reference  pattern, 
a loop  gain  constant. 

sources.  One  is  a Fourier  transformed 


l;.) 


Subtracted  I rom  this  is  a pattern  A(t.)  . A(f.)  is  the  consequence  ol  an  Instan- 

taneous reruns t rue l ion  process  which  is  on-going.  This  process  is  only 
impll. ii  in  Figures  b and  /.  The  conjugate  reference,  K* (?,),  is  used  to  form 
the  hologram,  as  shown.  I'he  unconjugated  reference  Rf!',)  is  used  t o reconstruct 
tin1  hologram  according  to  the  procedure 


A if.)  - K (O  II  (f.)  (All 

u n n 


An(f.)  is  inverse  transformed,  upon  demand,  and  thereby  produces  an  image, 
an(x),  appropriate  for  examination. 


By  subtracting  A t rom  S we  gel  R . 

n n n 


B (f.)  - S (O  - A (f.) 
n n n 


( U . ) 
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REFERENCE 

r(x> 


The  combination  of  Equat ions  41,  41,  and  44  results  In  a difference  equation. 
For  example,  using  Equations  41  and  44 


II 

n 


(S  - A )K  + 
n n n 


n-1 


Employing  Equation  43  In  45  gives 


II 

n 


E 


S K* 
n n 


(4  4) 


(4b) 


Wo  can  solve  for  II 


Vl  + * S„\ 

' + e|h,|j 


(47) 


This  is  the  proper  Iterative  difference  equation  to  describe  the  operation  of 
the  adaptive  hologram.  It  is  unconditionally  stable.  Notice  that  in  the 
derivation  ot  this  equation  it  was  assumed  that  the  formation  ol  the  increment 
to  the  hologram  ensemble,  g(Sn  - A1a)R*,  used  a reconstruction,  An,  which  itself 
depended  on  the  increment.  Tills  seems  like  a bootstrap  operation  and  is 
counter  intuitive.  The  procedure  is  known  as  backward  differencing  and  is 
essential  for  stability.* 


Suppose  we  were  to  use  a more  intuitive  approach.  Eet  us  make  the  subject 
channel  input  depend  on  the  reconstruction  of  what  has  already  been  recorded. 
Thus 
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Now,  we 

subs  t i t tie 

Equat ion  43 

into 

Equation  48 

to  get 
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Equat ion 

50  is 

t lu' 

forward  diff 

oroneo 

equ.it  ion 

version  ot  the  process. 

It  is 

unstable  for  positive  g when 

g R * R . > 2 , (SI) 

n n- 1 

For  negative  g it  Is  always  unstable,  which  is  also  true  tor  the  backward 
difference  case  of  Equation  47. 


*Oppenhelm,  A.V.,  and  Schafer,  S.W.,  "Digital  Signal  l'rocessing, " 
Englewood  Cliffs,  New  Jersey:  Prentice-Hall  (1975),  pp.  203-206 


In  the  initial  series  of  experiments  on  the  adaptive  hologram  the  forward 
difference  Equation  50  was  employed.  We  found  that  the  process  would  converge 
for  awhile,  but  then  would  begin  to  raptdlv  diverage.  We  were  using  random 
reference  patterns  for  these  experiments.  In  order  to  get  the  convergence 
time  within  reasonable  limits,  the  value  of  g was  made  moderately  large.  It 
would  turn  out  that  some  spatial  frequency  component  of  the  power  spectrum  ol 
the  reference,  |Rn(£)|  » would  be  large  enough  to  satisfy  Equation  51.  This 
spatial  frequency  would  exhibit  positive  feedback  and  the  system  would  blow 
up.  By  switching  to  the  backward  difference  Equation  47,  we  completely  cured 
the  problem. 

3.1.3  CONTINUUM  CASE 

In  order  to  gain  further  insight  into  the  properties  of  the  adaptive  hologram, 
let  us  carry  the  differencing  process  to  the  infinitesimal  limit.  It  the  step 
interval  between  n-1  and  n represents  a time  At,  we  may  rearrange  Equation  4t> 
and  write: 
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Carrying  this  to  the  limit  we  get: 
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Equation  51  is  the  most  general  description  of  the  adaptive  holographic 
process.  Witli  specialization  of  the  terms  we  can  draw  out  some  interesting 
results.  We  choose  as  a special  case  a reference  and  subject  which  both 
switch  on  to  a constant  level. 
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If  11  = 0 for  t < 0,  the  solution  ot  Equation  54  is  straightforward,  and  results 
in  the  equation: 
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3 . L . 4 1 MPORTANT  PROPERT I ES 

There  are  two  things  to  note  about  Equation  55.  First,  the  hologram  converges, 
in  the  limit  of  infinite  time,  to  the  spatial  pattern. 

. 8(0  K*(f.)  s(0 

|R(0  1 2 R(°  (5b) 

We  see  that  tlie  reference  pattern  is  thrown  into  the  denominator. 

Reconstruction  of  the  hologram,  using  this  reference,  is  described  by: 

A(Q  - R(U  H(fJ  - R „ s(f>)  (57) 

| R(£)  | 2 R(f;) 

Thus,  the  reconstruction  is  a perfect  reproduction  of  the  subject,  regardless 
of  the  funct ion  which  is  used  for  the  reference. 

This  is  a substantial  Improvement  over  the  conventional  recording  and  recon- 
struction process.  Even  if  a completely  random  function  Is  used  for  the 
reference,  there  will  be  no  background  hare  with  this  type  of  hologram.  In 
terms  of  the  description  give  by  Equation  11, 

rK(x)  © Tj  (x)  - <$K[  A(x)  + PK(  (x)  , (11) 

the  random  crosstalk  term,  P (x) , is  forced  to  zero. 

Kli 

We  can  regard  this  outcome  as  a sldelobe  cancellation  process  whereby  the 
sidelobes  of  the  point  spread  function,  produced  through  the  autocorrelation 
of  the  reference,  are  suppressed. 

(It  has  not  escaped  notice  that  the  functional  form  of  Equation  5t>  is  essen- 
tially the  same  as  that  employed  in  image  restoration  work.  It  may  very  well 
be  that  a variant  of  this  technique  will  prove  useful  for  image  restoration.) 


The  second  issue  of  consequence  arises  from  consideration  of  the  functional 
ratio  of  Equation  5b.  Selection  of  an  arbitrary  function,  r(x)  , to  serve  as 
a reference  will  usually  result  in  the  transformed  pattern,  R(Q , taking  on 
the  value  zero  for  certain  frequencies,  f,.  We  might  expect  that  the  hologram, 
apparently  having  zeros  in  the  denominator,  will  be  infinite  in  value  at  some 
points.  In  fact,  this  does  not  happen. 

The  resolution  of  this  problem  is  found  by  examination  of  the  time  constant  of 
the  exponent  in  Equation  55: 


The  time  constant  is  a function  of  spatial  frequency  which  is  inversely  pro- 
portional to  the  power  spectral  density.  For  those  spatial  frequencies  at 
which  R(f.)  * 0,  the  time  constant  is  infinite  and  no  change  in  the  hologram 
takes  place  during  the  adaptation. 

There  are  interesting  correlaries  of  this  important  result.  First,  all  por- 
tions of  the  hologram  do  not  adapt  at  the  same  rate.  Thus,  the  hologram  is 
not  formed  as  a unit.  Second,  spatial  frequencies  where  the  power  spectral 
density  is  highest  will  adapt  most  quickly  and  come  to  equilibrium.  For  most 
reference  scenes  these  will  usually  correspond  to  the  dc  background  and  slowly 
varying  terms.  High  spatial  frequency  terms  are  usually  low  in  amplitude  and 
will  adjust  slowly.  This  means  that  finely  detailed  structures  and  features 
will  have  influence  only  after  a substantial  time  has  passed. 

These  results  may  be  summarized  by  the  rule  that  the  basic  form  and  shape  (low 
frequencies)  of  an  object  will  tend  to  impact  the  hologram  before  the  fine 
details  and  texture  (high  frequencies).  This  rule  will  be  violated  only  when 
the  high  spatial  frequencies  dominate  the  power  spectrum  as  when  a scene  is 
composed  mostly  of  contrasty  textures. 

A final  correlary  is  that  spatial  frequencies  of  a subject  image  will  not  be 
recorded  unless  the  corresponding  spatial  frequencies  exist  in  the  reference. 
This  rule  is  true  for  conventional  holography  also. 

These  properties  of  the  adaptive  hologram  will  be  of  great  significance  if  the 
concept  is  employed  for  pattern  recognition. 

3.1.5  DISTINGUISH  INC.  BETWEEN  TWO  PATTERNS 

Given  the  interesting  properties  of  the  adaptive  hologram,  it  becomes  rela- 
tively simple  to  set  up  a global  training  paradigm  to  train  it  to  distinguish 
between  two  similar  patterns.  The  patterns  to  be  recognized  are  fed  in  through 
the  reference  channel.  Associated  with  each  reference  is  a subject  pattern 
which  serves  to  label  its  associate. 

We  did  a number  of  experiments  in  which  these  patterns  were  fed  to  the  system 
in  different  orders  with  variable  amounts  of  adaptation  permitted.  In  a typical 
experiment,  a particular  reference-subject  pair  was  fed  in  and  many  iterative 
cycles  were  computed  so  as  to  permit  the  adaptation  to  go  nearly  to  completion. 
Then  a second  reference  subject  pair  would  be  introduced,  again  with  sub- 
stantial adaptation.  At  the  end  of  a double  pairing  process,  a test  would  be 
made  to  see  whether  or  not  false  reconstructions  occurred.  This  test  con- 
sisted of  applying  each  reference  in  turn  and  seeing  whether  or  not  a trace  of 
the  unassociated  subject  appeared.  We  usually  found  that  it  was  necessary  to 
flip  back  and  forth  between  the  associated  pairs  several  times  before  all 
traces  of  false  reconstruction  were  rendered  invisible  (that  is,  suppressed 
below  the  lowest  level  of  the  print  gray  scale). 
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Wo  also  found  that  the  disappearance  and  form  ot  this  crosstalk  was  path 
dependent.  Thus,  it  made  a difference  how  much  adaptation  was  permitted 
before  the  reference-subject  pairs  were  exchanged.  This  peculiar  tosult  led 
to  some  doubts  about  the  nature  ot  the  pattern  recognition  process.  It  seems 
that  the  tlram-Schmidt  explanation  outlined  in  the  proposal  ( and  Appendix  A) 
may  not  be  correct.  In  order  to  verity  the  validity  ot  the  experiment  an 
extstance  proof  was  constructed.  This  proof,  given  lit  Appendix  H,  shows  that 
the  adaptive  hologram  can  be  trained  to  completely  distinguish  between  two 
partially  correlated  reference  patterns. 

1. 1.5.1  Critical  Issue,  Superpacking.  The  generalization  ol  this  property  ot 
different tal  recognition  to  separation  between  three  or  more  stmtlat  patterns 
has  not  been  demonstrated,  as  vet.  I’root  that  this  Is  possible  for  an 
arbitrary  number  ot  patterns  would  have  verv  important  consequences . 

Perhaps  the  most  significant  consequence  would  be  that  the  amount  ot  informa- 
tion which  could  be  stored  in  a hologram  would  be  vastly  increased.  It  was 
demonstrated  in  an  earlier  section  that  the  storage  capacity  ot  a conven- 
tionally multiplexed  hologram  is  proportional  to  the  number  of  degrees  ot 
freedom  of  the  system  in  which  the  hologram  is  embedded.  The  limitation  was 
established  primarily  by  the  reconst  met  ion  crosstalk  between  various  members 
of  the  multiplexed  ensemble. 

It  the  adaptive  apparatus  were  used  t o term  a multiplexed  hologram  ensemble  it 
mav  be  possible  to  completely  cancel  the  crosstalk  between  ensemble  members. 

It  this  proves  to  bo  the  case,  then  the  storage  capacity  of  the  hologram 
increases  to  the  square  ot  the  number  ot  degrees  ot  treodom. 

The  argument  behind  this  conclusion  runs  as  follows:  The  hologram  represents 
a vector  space  with  n degrees  ot  treodom.  We  mav  thorotoro  embed  n vectors  in 
this  space.  I'.ach  vector  represents  a different  reference-subject  pairing. 
During  t he  reconst  met  ion  the  subject  vector  is  reproduced  along  with  a cross- 
talk term.  It  the  crosstalk  can  he  canceled  awav,  then  the  subject  vector 
will  be  in  pristine  form.  It  lt\e  sub|eet  vector  also  has  n degrees  ot  freedom, 
it  therefore  can  represent  thfs  many  bits  ot  data.  Since,  wiilt  crosstalk 
cancellation,  there  can  be  n ot  tltose  subject  vectors  reconstructed  j'v  t lie 
system,  the  total  number  ot  data  bits  which  can  be  stored  is  n x n • u* . 

This  hypothesis,  ot  course,  has  vet  t o be  validated.  It  it  proves  to  he  cot 
root,  it  will  have  at  least  the  following  Important  consequences . First,  the 
amount  of  effort  to  pack  the  hologram  tull  ot  Information  will  be  verv  great. 
This  is  because  the  crosstalk  cancel lat ton  process  involves  an  iteration  among 
t he  contributing  reference-subject  pairs.  Moreover,  t ho  more  pairs  are  to  be 
stored,  the  greater  the  number  ot  iterations  that  will  probably  be  required  to 
achieve  convergence . 


Second,  addition  of  information  to  an  adapted  hologram  requires  a great  deal 
of  adjustment.  The  new  image  will  generate  cross  terms  with  all  the  other 
images  previously  stored.  These  will  have  to  be  adapted  out  one  by  one. 
Moreover,  this  first  order  adaptation  process  will  in  turn  create  new  cross 
terms  which  will  require  at  least  second  order  correction.  Clearly,  the 
amount  of  iteration  required  to  stuff  a new  pair  into  an  ensemble  will  be  a 
very  strong  function  of  the  number  of  old  pairs  already  making  up  the  ensemble. 

A third  consequence  is  that  the  adaptive  hologram  may  provide  an  immensely 
powerful  mechanism  for  pattern  recognition.  This  is  because,  if  crosstalk 
cancellation  happens,  the  number  of  reference  patterns  which  may  be  recognized 
will  be  equal  to  the  number  of  patterns  which  may  be  embedded.  Which  is,  of 
course,  the  number  of  degrees  of  freedom. 

Consider  a practical  case.  It  is  within  the  state  of  the  art  to  provide 
Fourier  transforms,  over  an  array  of  409b  elements,  many  times  a second.  If 
the  above  argument  is  correct,  a pretrained  system  of  this  size  could  be  made 
to  distinguish  among  as  many  as  409b  similar  and  different  patterns  in  a 
single  pass  (milliseconds). 

The  possibility  that  this  superpacking  of  the  hologram  may  work  sets  a clear 
direction  for  follow-on  research.  We  must  discover  if  full  cross  talk  can- 
cellation is  possible  or  if  the  rules  of  conventional  multiplexing  set  the 
limit.  If  superpacking  proves  possible,  then  an  investigation  should  be  made 
to  find  an  efficient  training  scheme.  Once  this  basic  research  work  is 
successfully  completed  practical  applications  will  rapidly  follow  because  the 
technology  base  is  ripe. 

3.2  EXPERIMENTAL  RESULTS 
3.2.1  TESTS  OF  ADAF’TATION 

The  discovery  of  the  adaptive  hologram  substantially  modified  our  original 
experimental  plan.  We  had  intended  to  move  straight  into  an  exercise  of  train- 
ing a hologram  to  distinguish  between  two  similar  patterns.  Instead,  we 
decided  to  first  explore  some  of  the  convergence  and  clutter  rejection  proper- 
ties of  the  adaptive  system. 

Our  initial  series  of  experiments  employed  a single  subject-reference  pair. 

The  subject  was  a data  array  of  sixteen  uniform  amplitude  points,  and  the 
associated  reference  was  a random  pattern.  The  choice  of  these  images  was 
made  because  we  had  just  concluded  a detailed  properties  study  of  this  class 
of  holograms. 

After  some  wrestling  with  the  software  and  adjusting  the  feedback  loop  gain 
coefficient,  we  found  that  the  adapting  system  gradually  cleaned  up  the  clutter. 
This  was  a notable  and  exciting  event.  Fvirther  iteration  resulted  in  diver- 
gence, however.  The  reconstructed  subject  patterns  showed  characteristic  peri- 
odicities indicating  that  the  divergence  was  occurring  at  localized  points  in 
the  Fourier  plane.  Figure  8 shows  a representative  sequence  of  divergence. 
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Upon  investigation  the  problem  was  traced  to  our  use  of  a forward  difference 
equation  to  describe  the  feedback  process.  In  the  preceding  theoretical  sec- 
tion there  is  a detailed  discussion  of  this  issue.  By  switching  to  a backward 
difference  equation,  the  Instability  was  eliminated  and  no  further  problems 
of  this  kind  were  encountered. 

A typical  early  experiment  with  the  corrected  system  is  illustrated  by  Figures  9 
and  10.  Figure  9 shows  reconstructions  of  the  data  field  subject  with  the  image 
showing  the  characteristic  clutter  caused  by  the  autocorrelation  of  the  random 
reference  pattern.  Note  that  the  initial  signal-to-clutter  ratio  is  close  to 
the  ensemble  average  value  of  sixteen  predicted  by  theory.  The  other  entries 
in  Figure  9 graphically  demonstrate  the  improvement  in  reconstruction  quality 
as  the  system  iteratively  adapts. 

The  quantitative  results  of  this  experiment  are  illustrated  by  Figure  10.  The 
subject  amplitude  was  100.  Initially  the  reconstructed  image  had  a much  lower 
amplitude.  This  was  a consequence  of  the  values  chosen  for  the  magnitudes  of 
the  reference  pattern,  and  the  loop  gain  coefficient,  g.  Upon  iteration  the 
amplitude  of  the  reconstruction  exponentially  converged  to  that  of  the  subject 
in  accord  with  Equations  55  and  57.  The  signal-to-clutter  ratio  also  improved 
dramatically  during  this  convergence  process. 

3.2.2  TESTS  OF  CONVERGENCE 

In  the  initial  series  of  experiments  the  reference  was  a random  pattern.  This 
type  of  reference  had  approximately  a uniform  power  spectral  density.  The 
theory,  as  exemplified  by  Equation  58,  indicated  that,  in  this  case  all  spatial 
frequencies  in  the  reconstruction  would  converge  at  approximately  the  same  rate. 
This  was  apparently  the  case  in  the  initial  set  of  experiments. 


We  decided  to  try  a much  sterner  test  of  the  convergence  of  the  system.  We 
selected  a reference  pattern  almost  completely  dominated  by  low  spatial  fre- 
quencies with  the  high  frequency  terms  being  relatively  very  weak.  This  refer- 
ence was  paired  with  a subject  which  was  small,  complicated  and  composed  pri- 
marily of  very  high  spatial  frequencies.  In  conventional  holography  there 
would  be  no  chance  of  getting  a useful  reconstruction  of  the  subject,  but  the 
new  theory  predicted  chat,  after  adaptation,  a clear  reconstruction  was  possi- 
ble. The  reference  for  this  experiment  was  a two  dimensional  square  wave, 
half  filling  the  field  of  view.  The  Fourier  transform  of  this  function  has 
the  usual  form  of  Sin  £/£  and  thus  assures  there  is  at  least  some  signal  at 
the  highest  spatial  frequencies  of  the  spectrum.  This  function  also  contains 
regions  of  zero  signal,  which  provides  an  interesting  test  of  the  theory. 


The  subject  chosen 
(see  Figure  14B) . 
therefore,  most  of 
end  of  the  spatial 


to  be  associated  was  the  image  of  the 
The  structures  of  this  character  were 
the  signal  energy  was  concentrated  at 
spectrum. 
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Initial  frame  of  sequence. 
Reconstruction  has  S/N  of  15.11, 
compared  to  theoretical  of  16. 


Reconstruction  after  partial 
adaption,  showing  partial  clear) 
of  autocorrelation  produced 
clutter.  S/N  = 21.26 


C.  Continued  cleanup  of  clutter. 
S/N  - 24.30 


D.  l.ate  reconstruct  ion  in  adaptive 
sequence.  S/N  = 32.59 


FICURE  9. 


The  clutter  is  initially  generated 
through  the  autocorrelation  of  a 
random  reference  image.  No  multi- 
plexing occurred  in  this 
experiment . 


SEQUENCE  OF  RECONSTRUCTIONS  SHOWING  CLUTTER  REDUCTION  BY  FEEDBACK 
ADAPTION 


SUBJECT  AMPLITUDE 


RECONSTRUCTION  AMPLITUDE 


SIGNAL  TO  CLUTTER  RATIO 


INITIAL  VALUE 


NOTE  THE  CONVERGENCE  OF  THE  AMPLITUDE 
OF  THE  RECONSTRUCTED  SIGNAL  TO  THE  ORIGINAL 
SUBJECT.  THIS  CONVERGENCE  IS  ACCOMPANIED 
BV  A STEADY  REDUCTION  IN  CLUTTER 


NUMBER  OF  ITERATIONS 


FIGURE  10.  QUANTITATIVE  RESULTS  FOR  EXPERIMENT  ILLUSTRATED  IN  FIGURE  9 


Figure  11  illustrates  the  reference  used,  together  with  reconstructions  after 
various  numbers  of  iterations.  The  early  reconstructions  lack  any  useable 
detail,  as  indicated  by  Figure  11B.  Here  the  energy  of  the  low  spatial  fre- 
quencies in  the  reference  totally  dominate  the  holographic  process.  This 
image  is  representative  of  (although  somewhat  better  than)  conventional  holog- 
raphy. After  many  iterations  the  characteristic  shape  of  the  W becomes  visi- 
ble. Figure  11D  provides  an  example,  late  in  the  iteration  sequence  of  this 
particular  experiment.  There  are  some  peculiar  features  of  the  process  indi- 
cated by  this  reconstruction.  Since  the  effects  of  the  low  spatial  frequencies 
have  been  adapted  away,  the  high  frequency  terms  are  dominant.  Note  the  traces 
of  the  corners  of  the  reference  indicated  by  the  secondary  reconstructions  at 
the  left  and  bottom  edges  of  the  field.  These  reconstructions  illustrate  that 
cross  correlation  between  the  corners  is  very  strong.  Note  also  the  influence 
of  the  edges  of  the  reference  function  which  is  apparent  in  the  streaks  radiat- 
ing from  the  principal  (central)  reconstruction. 

Since  convergence  was  slow,  we  repeated  the  experiment  with  a much  stronger 
feedback  coefficient.  The  results  are  shown  in  Figure  12.  The  first  figure 
in  the  sequence  apparently  picks  up  at  the  end  of  the  previous  experiment. 

As  the  sequence  progresses  the  primary  and  secondary  reconstructions  become 
clear  and  accurate  representations  of  the  original  subject.  The  impact  of  the 
edges  of  the  reference  function  have  nearly  faded  by  the  end  of  the  sequence, 
and  the  reconstruction  process  is  dominated  by  the  corners. 

There  is  some  fluctuation  in  the  secondary  reconstructions  as  the  convergence 
process  progresses.  The  left-hand  reconstructions  start  to  increase  slightly 
in  strength  as  the  bottom  reconstruction  begins  to  fade.  Presumably,  the  left 
reconstructions  would  peak  and  fade  if  the  iteration  process  continued  long 
enough  (about  two  hours  of  CPU  time  was  used  for  this  experiment). 

A third  experiment  in  this  series  validated  the  preliminary  conclusions  drawn 
from  the  earlier  experiments.  For  this  experiment  the  square  wave  was  thinned 
down  to  just  its  edge.  The  reference  thus  took  the  form  of  an  empty  box  (see 
Figure  14D) . This  enhancement  of  the  high  spatial  frequencies  in  the  refer- 
ence greatly  speeded  the  convergence;  Figure  13  evidences  the  result.  In  the 
early  iterations  the  reconstruction  looks  much  like  the  late  forms  of  Figure  12. 
By  the  end  of  the  experiment,  the  principal  reconstruction  is  nearly  perfect, 
whereas  traces  of  the  secondary  reconstructions  are  nearly  eliminated. 

We  conclude  from  these  tests  that  the  major  features  of  the  theory  of  adapta- 
tion are  confirmed.  More  detailed  quantitative  studies  might  show  some  unex- 
pected phenomena.  However,  we  felt  that  there  would  be  a greater  return  for 
our  efforts  by  pressing  on  to  a test  of  differential  pattern  recognition. 

3.2.3  RECOGNITION  OF  DIFFERENCES 

Having  successfully  demonstrated  that  the  adaptive  hologram  follows  the  theo- 
etical  predictions,  the  experimental  effort  returned  to  the  originally  scheduled 
task.  As  described  in  the  foregoing  theoretical  section,  it  was  believed  that 
a hologram  could  be  trained  to  distinguish  between  two  similar  (highly  corre- 
lated) patterns.  With  the  machinery  of  the  adaptive  hologram  available,  a test 
of  this  idea  was  relatively  easy. 
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Reference  pattern.  This  pattern 
consists  mostly  of  low  spatial 
frequencies . 


Reconstruction  after  first  series 
of  iterations.  Note  dominance 
of  low  spatial  frequencies. 


Further  in  the  iteration  sequence. 
High  spatial  frequency  structure 
is  becoming  apparent 

This  experiment  attempts  to  reconstruct 
a small  highly  structured  subject  (the 
image  of  the  letter  W)  with  a much 
larger  nearly  unstructured  reference. 

In  conventional  holography  the 
reconstruction  would  show  virtually 
no  structure. 


Reconstruction  at  end  of  test  run 
The  shape  of  the  subject  is  becom 
ing  clear.  Reconstruction  is  now 
dominated  by  the  reference  edge 
structure.  Repeat  images  are  the 
result  of  high  cross  correlation 
among  the  corners  of  reference 
box  image. 


FIGL'RE  11.  ATTEMPTED  RECONSTRUCTION  OF  SMALL  SUBJECT 


A.  Initial  reconstruction.  The  result 
is  very  similar  to  the  final  result 
of  Figure  11.  This  experiment  is 
properly  regarded  as  a continuation 
of  that  experiment. 


B.  Reconstruction  after  75  iterations. 
The  W structure  is  becoming  clear. 
Repeat  of  the  image  is  caused  by 
corner  cross  correlations. 


C.  Continued  clean  up  of  main  image  D. 

and  further  weakening  of  secondary 
reconstructions . 

Experiment  similar  to  that  of  Figure  11, 
but  with  increased  loop  gain.  The  final 
reconstruction  clearly  shows  the  desired 
image  in  the  center  of  the  frame. 

Secondary  reconstructions  should  fade 
out  after  further  iteration. 


Final  Image  in  experiment.  Basic 
trends  continue  to  prevail  except 
for  a slight  enhancement  of  the 
upper  left  secondary  image. 
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FIGURE  12.  ATTEMPTED  RECONSTRUCTION  OF  SMA1.L  SUBJECT  WITH  INCREASED  LOOP  GAIN 


Reconstruction  at  beginning  of 
iteration  series.  Note  Similarity 
to  12A.  Secondary  reconstructions 
start  out  weaker  in  this  series. 


b.  Continued  clean  up  of 
reconstruction . 


Further  reduction  of  secondary  D.  I- Inal  reconstruction  In  test 

reconstructions  sequence.  Primary  image  is  almost 

completely  clean.  Secondary 
reconstructions  are  nearly 
el iminated. 


FIGURE  13.  EXPERIMENT  SIMILAR  TO  FIGURES  11  AND  12,  EXCEPT  THAT  THE  REFERENCE 
CONSISTED  OF  JUST  THE  EDGE  OF  THE  BOX  (SEE  FIGURE  14D) 


An  experiment  was  designed  in  which  two  different  subject-reference  pairs  were 
presented,  alternately,  to  the  system  for  training.  Upon  presentation  of  one 
such  pair,  the  system  was  allowed  to  adapt  for  several  iterations  and  a test 
was  then  made  of  the  reconstructive  properties  of  the  partially  trained  holo- 
gram. This  consisted  of  separately  presenting  the  two  reference  patterns  to 
see  what  form  the  reconstructions  would  take.  After  this  test  the  alternate 
reference-subject  image  pair  was  injected  and  the  iteration  continued. 

Typical  results  are  illustrated  in  Figures  14  and  15.  In  Figure  14  the  refer- 
ence and  subject  images  are  shown.  The  images  were  chosen  based  on  previous 
experience  with  the  system  convergence  properties,  and  by  the  need  for  simplic- 
ity to  avoid  tedious  hand  entry  of  data. 

As  illustrated  a narrow,  rectangular  box  served  as  the  reference  for  the  sub- 
ject letter  image  N.  A wide,  square  box  was  selected  as  the  reference  to  be 
paired  with  the  letter  image  W.  The  purpose  of  the  experiment  was  to  show 
that  after  suitable  training,  according  to  the  simple  paradigm  just  described, 
the  narrow  box  would  reconstruct  only  the  image  of  the  N,  whereas  the  wide  box 
would  reconstruct  only  the  W. 

The  cross  correlation  coefficient  between  the  two  reference  patterns  was  calcu- 
lated to  be  0.689,  so  a failure  of  the  system  to  distinguish  between  the  two 
reference  patterns  would  show  up  as  strong  secondary  reconstructions.  In  order 
to  make  these  shadow  reconstructions  evident,  the  N and  W images  were  located 
in  slightly  different  positions  within  the  field  of  view. 

Figure  15  demonstrates  the  success  of  this  experiment.  Figures  15A  and  15B 
show  intermediate  results  after  a certain  amount  of  adaptation  has  occurred. 

The  desired  primary  reconstructions  are  becoming  clear.  Substantial  secondary 
reconstructions  of  the  unpaired  subject  shows  the  correlated  subset  of  the  two 
references  has  not  been  entirely  cancelled  out. 


The  last  two  subfigures  picture  the  state  of  the  system  after  further  adapta- 
tion. The  reference  subject  pair  used  immediately  before  this  test  reconstruc- 
tion was  the  W-square  box.  This  shows  in  the  perfectly  reconstructed  image  of 
the  W.  It  also  shows  in  the  slight  imperfections  of  the  N image  when  the 
narrow  box  reference  is  played  through  the  hologram.  Although  this  was  the 
end  of  the  test  sequence,  convergence  was  not  yet  complete.  Evidently  spatial 
frequencies  of  low  amplitude  were  still  impacting  the  reconstruction  process 
and  causing  false  reconstructions. 

Other  experiments  of  a rough  sort  indicated  that  convergence  is  path  dependent. 
That  is,  it  matters  what  the  strength  of  the  feedback  is  or  how  many  iterations 
are  made  before  switching  reference-subject  pairs.  The  degree  of  cross  correla- 
tion between  different  reference  patterns  also  appears  to  have  an,  as  yet, 
undetermined  impact.  Both  of  these  topics  may  be  important  and  should  have 
priority  in  future  investigations. 


A.  Reconstruction  using  the  narrow 
box  image  as  reference.  Note  the 
false  reconstructions.  This  image 
was  produced  after  a sequence  of 
alternating  subject  reference 
pairs  had  been  presented  to 
partially  train  the  system  in 
differential  recognition. 


B.  Reconstruction  with  wide  box  image 
after  partial  training.  Note 
false  reconstructions  of  N near 
the  W.  Note  also  the  false 
reconstructions  due  to  the  other 
corners  of  the  box. 


Reconstruction  with  training 
almost  completed.  The  last 
iteration  sequence  had  used  a 
W-wide  box  combination.  This 
explains  the  residual  cross  talk. 


n. 


Reconstruction  with  wide  box 
reference  and  the  system  in  the 
same  state  of  training  as 
view  C. 

The  system  has  learned  to  distinguish 
between  reference  patterns  that  have 
a O.bfW  correla t ion. 


FIGURE  15.  RESULTS  OF  DIFFERENTIAL  TRAINING  EXPERIMENT 
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The  software  included  in  this  document  is 
configured  for  a Data  General  Eclipse  Computer 
supported  with  a 100  M Byte  disk  memory. 


NOTE 


SECTION  4 


MACHINE  HOLOGRAPHY  SOFTWARE 


4 . 1 SUMMARY 

The  machine  holography  exercise  explored  the  data  storage  properties  and  some 
image  recognition  properties  of  holograms  through  use  of  a computer  simulation. 
The  software  developed  for  the  hologram  simulation  consists  of  five  main 
classes  of  programs: 

(1)  Image  creation  routines 

(2)  Hologram  routines 

(3)  I/O  routines 

(4)  Utility  routines 

(5)  Borrowed  and  miscellaneous  routines 

All  machine  holography  software  resides  on  disk  under  directory  HOLO.DR.  The 
directory  is  saved  on  magnetic  tape.  File  0 of  the  tape  contains  the  holo- 
graphy programs,  and  File  1 contains  the  holography  image  files. 

The  image  creation  routines  create  various  simple  point  images  and  random 
noise  scenes  in  the  integer  disk  file  format  used  by  the  machine  holography 
software.  The  image  creation  routines  have  file  names  of  the  following 
form:  CR — .FR.  All  image  creation  software  is  a variation  to  the  CRSCN.FR 

program  which  creates  a point  image  at  specified  coordinates  in  an  integer 
array  stored  on  disk.  CRNOI.—  program  creates  random  noise  integer  arrays 
in  the  standard  disk  format  used  by  machine  holography  software.  There 
appeared  to  be  some  correlation  between  different  random  scenes,  so  the 
random  scenes  were  input  to  the  program  SCRAM. FR  to  "scramble"  (randomly 
shift)  the  rows  and  columns  of  the  data  file.  The  following  are  scene 
creation  routines: 


1 . 

CRSCN.— 

Creates  single  point  image  file;  ex.  file  P3333. 
point  at  coordinates  X = 33  Y = 33 

f - 

l; 

CRNOI.— 

Creates  a random  noise  scene;  ex.  file  N0ISE1  is 
random  noise  file 

Li 

SCRAM. - 

Used  to  "randomize"  a random  scene.  This  program 
modifies  a file  that  had  been  previously  created. 

L 

CR4SCN .— 

Creates  a four  point  image;  ex.  file  F3232  is 
a file  of  four  points  with  LH  corner  at  X = 32, 

Y = 32 

f • 
!. 

i. 


i 

r 

, 


CRlbSCN.— 


CRLN01 

CRN.— 


CRW.— 


Creates  a sixteen  point  scene;  ex.  S3232  is  .1 
file  of  sixteen  points  with  LH  corner  at 
X-32,  Y-32 

Creates  a random  noise  scene  of  limited 
site  (l.bAxbA) 

Creates  an  linage  "N";  ex.  file  1MAGEN  was 
created  by  CRN.SV 

Creates  an  linage  ot  "W" ; ex.  tile  1MA0EW 
was  created  by  CRW.SV 


CR3. 

CRA. 

CRS. 

CRb. 


CRSQ. 


Creates  an  image  file  of  a 

and  an  integer.  These 
were  used  to  create  the 
Tank  1.0.  files:  T1D1 . TID2, 
TID3,  TIDA , T10S.  and  Tll)6. 
Creates  a square  image 


All  image  tiles  can  be  examined  by  executing  the  program  IMAGER. SV.  IMAGER 
reads  the  image  from  disk  and  outputs  a grayscale  plot  to  the  line  printer. 
The  subroutine  SDEF  (contained  in  disk  file  1REA0.  ) is  used  bv  machine 
holography  software  to  read  the  image  files  from  disk  into  memorv. 

The  actual  holography  routines  consist  ot  tour  programs:  HDRV.  , HSDRV. 
MH0L0.-  EBORVT . , and  CRSUM.  . 


110RV . — 

HSDRV .- 
MH0L0. 


KBORVT. 


Creates  a single  hologram  and  performs 
an  image  reconstruct  ion 

Creates  a multiplexed  hologram 

Reads  a hologram,  performs  an  image  reconstruction, 
and  outputs  a grayscale  plot.  MH01.it.  is  used  in 
conjunction  with  HSDRV. 

Creates  a hologram  feedback  loop 


CRSUM.—  Is  a utility  program  used  to  zero  a complex  array 

file  on  disk.  CRSUM.  is  executed  to  create  a zero 
valued  hologram  file  prior  to  running  HSDRV. - 


I/O  routines  are  RBIX.— , RB1.KC.— , WRI.C.  , and  WB1.KC.  . These  routines  are 
used  to  read  and  write  bA  by  bA  complex  arrays  to  and  from  disk. 

Several  routines  were  borrowed  from  different  directories.  These  are 
FFT2D.RB,  GSPLT .— , PICPLT.RB.  and  Pl.PT.RB. 


r 

l 

0 

t 

f 
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FFT2D.RB  Performs  a two-dimensional  Fourier  transform 

on  a complex  array 


GSPLT 


GSPLTX.— 


PLPT.RB  and 
PICPLT.RB 


Performs  a grayscale  plot  of  a complex  data 
array.  The  program  was  modified  to  calculate 
S/N  of  the  array. 

Is  a modified  version  of  GSPLT.—  which  has  a 
parameter  in  the  call  statement  specifying 
how  many  points  are  used  in  the  S/N  calculation 

Are  used  by  GSPLT  and  GSPLTX 


The  software  operates  on  64  by  64  complex  arrays.  Because  of  limited  memory, 
as  the  programs  execute  they  swap  data  back  and  forth  from  memory  to  disk. 

The  driver  programs  create  files  and  read  existing  files.  The  operator 
specifies  the  required  file  names  when  prompted  by  the  programs.  The  file 
names  are  stored  in  seven  element  integer  arrays,  so  file  names  are  limited 
to  seven  characters.  Any  I/O  error  causes  the  programs  to  stop. 

The  utility  subroutines  and  the  disk  I/O  subroutines  (except  SUBROUTINE  SDEF) 
use  the  format  shown  in  Figure  16  to  read  and  write  data.  The  disk  file  uses 
sixty-five  256-word  blocks. 

NOTE: 

(1)  All  file  names  used  are  seven  characters  or  less. 

(2)  File  names  are  stored  by  the  program  in  integer  arrays. 

(3)  Data  stored  on  disk  files  are  in  standard  format,  i.e., 
either  integer  format  or  complex  format 

(4)  In  the  following  write-up,  upper  case  letters  correspond  to 
program  console  output.  Underlines  lowercase  letters  correspond 
to  operator  console  input. 


(5) 


means  carriage  return. 


HEADER 


64  COMPLEX  VALUES 
266  WORDS 


• DATA  IS  STORED  IN  MEMORY  IN  THE 
COMPLEX  ARRAY  IS  (64.  64) 

/ \ 

COLUMN  ROW 

INDEX  INDEX 

• DATA  IS  STORED  IN  65  BLOCKS  ON  THE 
DISK  FILE.  EACH  BLOCK  CONTAINS  ONE 

ROW  OF  OATA  (64  COMPLEX  VALUES  256  WORDS) 
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FIGURK  16.  DATA  FORMAT  OF  COMPLEX  DISK  FILE  DATA 


4.2  IMAGE  CREATION  ROUTINES 


All  images  and  scenes  used  by  the  machine  holography  software  are  stored  on 
integer  disk  files  according  to  the  format  in  Figure  17.  The  disk  file 
stores  a two-dimensional  array.  Each  row  of  the  array  is  stored  in  one  block 
of  the  file  with  consecutive  rows  stored  in  consecutive  blocks.  Each  block 
is  256  words,  with  each  row  of  the  image  having  256  integer  values.  All 
machine  holography  images  are  64  by  64,  so  only  a quarter  of  each  block  is 
used.  This  is  shown  by  the  shaded  area  in  Figure  17.  The  image  file  can 
have  an  arbitrary  number  of  rows  (blocks)  of  data.  Machine  holography  images 
have  b4  rows  of  data.  All  machine  holography  images  have  array  coordinates 

1.1  at  the  upper  left  hand  corner  of  image  array  file.  All  integer  files  are 
read  using  IMAGER. SV  (for  plots)  or  SUB  SDEF  (contained  in  file  IREAD.-). 

4.2.1  CRSCN.FR 

CRSCN.FR  is  a FORTRAN  program  that  generates  an  integer  image  file.  All 
machine  holography  programs  of  the  form  CR— ,FR  (except  CRSUM)  are  modifications 
of  CRSCN.FR. 

CRSCN.FR  creates  an  integer  image  file  consisting  of  a single  point  of  a 
specified  value,  at  specified  coordinates,  with  the  remaining  points  of  the 
image  set  to  zero.  All  CR— .FR  programs  (except  CRSUM)  ask  for  the  image  array 
size  and  position  within  the  data  file  (refer  to  Figure  17).  All  machine 
holography  images  are  currently  64  by  64  in  size  and  are  positioned  with  the 
upper  left-hand  corner  positioned  at  coordinates  1,1. 

Modifications  to  CRSCN.,  in  order  to  create  images  other  than  a single  point, 
are  made  by  modifying  the  FORTRAN  code  in  the  DO  50  loop.  This  can  be  best 
understood  by  examining  the  listings  of  the  different  CR— .FR  programs.  The 
following  is  the  dialogue  between  the  program  and  operator. 

RUN  n 

COLUMN  integer 
ROW  integer 

(NOTE:  Column  and  row  are  the  coordinates  of  the  image  point  within 
the  image  array. ) 

VALUE  real  number 

(NOTE:  Value  is  the  "intensity"  of  the  point.  VALUE  = 100.  for  most 
machine  holography  images.) 


I 


C CRSCN  FR  7/05/78 

C 

C CRSCN  GENERATES  AN  256  X JRMX  ARRAY  CONTAINING 

C AN  NX  X NY  IMAGE  WITH  UPPER  LH  CORNER  COORDS  (NXC.NYC) 

C THE  PORTIONS  OF  THE  ARRAY  NOT  CONTAINING  THE  IMAGE  ARE 

C FILLED  WITH  ZEROS  AND  THE  ARRAY  IS  STORED  ON  DISK  AS  KNAME 

C 

DIMENSION  I DD  < 256 ) , KNAME ( 7 > 

C 

C KNAME  IN  FILE  #1  IS  USED  TO  STORE  THE  BACKGROUND  IMAGE 

DO  00  100-1. 10 
TYPE  "RUN  ",  100 

ACCEPT  "COLUMN  ",NCOL. "ROW  ", NROW, "VALUE  ".VALUE 
WRITE! 10, 900) 

900  FORMAT ( " NAME  OF  FILE  TO  BE  CREATED  - ",Z> 

READ  <11, 905 ) KNAME ( 1 ) 

905  FORMAT (S7) 

CALL  DF I LW ( KNAME ( 1 ) , I ER ) 

CALL  CF IL W ( KNAME ( 1 ) , 2, IER) 

IF ( IER  NE  1)  STOP  — CFILW  ERROR  (BCD)  — 

CALL  OPEN( 1, KNAME(1 ). 2, IER, 512) 

IF i IER  NE  1)  STOP  — OPEN  ERROR  (BCD)  — 

C 

C INPUT  IMAGE  PARAMETERS 

C 

ACCEPT" IMAGE  SIZE  (NX. NY)  - ",NX,NY 
1 YPE"COORDS  OF  UPPER  LH  CORNER  OF  IMAGE" 

TYPE"E  C (1,1)-  UPPER  LH  CORNER  OF  ENTIRE  ARRAY" 

ACCEPT" (NXC, NYC ) - ", NXC.NYC 

JCMX-NX+NXC-1 

JRMX-NY+NYC-1 

IF  < JCMX  GT.  256)  STOP-- IMAGE  OFF  RIGHT  EDGE  OF  ARRAY- 

PI -3  1415926 
ARG=PI /8 
DO  10  1=1,256 
10  I DD ( I ) =0 

IDD ( 2 > =256 
IDD ( 4 ) -JRMX 

CALL  WRBLK( 1, 0, IDD, 1, IER > 

IF(  IER  NE.  1 ) STOP — WRBLK  ERR— 

IDD( 2 ) =0 
I DD ( 4 ) =0 
NS 1 =NYC- 1 

lF(NSl  Lfc  U)  UO  lO  dO 

DO  20  IR-1.NS1 

20  CALL  WRBLK ( 1 , IR,  IDD,  1,  IER) 

IF  (IER  NE  1)  STOP— WRBLK  ERR  — 

30  CONTINUE 

RMEAN-0 
SICMA-256 
NRAND-235 1 
DO  50  I R=NYC , JRMX 
1 IR-IR  + NYC-1 
DO  5 IC=1 , NX 

I IC  = IC+NXC-1  ) 

IDD ( I IC ) =0 
5 CUNT 1NUE 

IF ( IR  EG  NROW) IDD (NCOL) -VALUE 
CALL  WRBLK ( 1 , IR,  IDD.  1,  IER) 

IF  i IER  NE  1)  STOP  —WRBLK  ERROR  ( BGD ) — 

50  CONTINUE 

CALL  CLOSE! 1, IER) 

00  CONTINUE 


STOP 

CTKjn 


n n non  noon 


M 


CRNOI  FR  7/05/78 

ORSON  GENERATES  AN  256  X JRMX  ARRAY  CONTAINING 

AN  NX  X NY  IMAOE  WITH  UPPER  LH  CORNER  COORDS  (NXC.NYC) 

THE  PORTIONS  OF  THE  ARRAY  NOT  CONTAINING  THE  IMAOE  ARE 
Fil  l ED  WITH  ZEROS  AND  THE  ARRAY  IS  STORED  ON  D.SK  AS  KNAME 

DIMENSION  IDD1256). KNAME (7) 


KNAMF  IN  FILE  #1  IS  USED  TO  STORE  THE  BACKGROUND  IMAOE 

RMEAN-C 
SIGMA*256 
NRAND  2151 

ACCEPT  "RANDOM  SEED" . NRAND 
DO  00  100-1. 10 
TrPE  "RUN  ",  180 
WRITEl 10. TOO) 

TOO  FORMAT!"  NAME  OF  FILE  TO  BE  CREATED  - ".Z> 

READ* 11. 905)  KNAME l 1 ) 

905  FORMA T(S7) 

OAl  L DF 1 1 W ( KNAME ( 1 ) . IER) 

CALI.  CFILWIKNAMEl  1 >.  2.  IER) 

IF t IER  NE  1)  STOP  — CFILW  ERROR  <BOD)  — 

CAM  OPEN! I.  KNAMEl 1 ).  2.  IER.  512) 

11  l(R  Nl  I)  STOP  --OPEN  ERROR  « BOD ) — 

C 

c input  image:  PARAMETERS 

0 

Al  Cl PT" IMAGE  SIZE  iNX.NY)  - “ , NX . NY 
TYPE  COORDS  OF  UPPER  LH  CORNER  OF  IMAOE" 

TYPL'E  G il. 1)-  UPPER  LH  CORNER  OF  ENTIRE  ARRAY" 
ACCEPT" <NXC. NYC)  - ".NXC.NYC 
v)CMX-NX-*NXC  1 
URMX “NY  “NYC  - 1 

IFUCMX  GT  256>  STOP  -IMAGE  OFF  RIGHT  EDGE.  OF  ARRAY  — 

PI-0  1415926 
ARG-PI  '0 
DO  10  1-1.256 
10  lDOiD-0 

i DD l 2 i -256 
1 DD<  4 ) -JRMX 

CAl L NRBLKl 1.0. I DD. I . IER! 

IF  (HR  NE  1 )STOP-~WRBLK  ERR-- 
1 DD i 2 ) -0 
I DD i 4 ) -0 
NS  1 -NYC  1 

IF  l NS  1 LE  0)  00  TO  JO 
DO  20  IR-l  NS1 

20  CAl  L Wf.’BLKl  1.  IR.  I Dl>.  1.  IER) 

IF  l IER  NE  1)  STOP — WRBLK  ERR  — 

30  CONTINUE 

Di  50  IK-NVC. UK MX 
I I R - 1 R * NYC ■ 1 
DO  0 IC-1 . NX 
I IC-1C+NXC-1 

IDDi  I 1 C ) -RGAUSl RMCAN, SIGMA.  NRAND) 

5 CONTINUE 

CALL  WRBLK ( 1 , IK.  IDD  1.  IER) 

IF  < IER  NE  1)  STOP  --WRBLK  ERROR  l BCD ) -- 
50  CONTINUE 

CALL  CLOSl  l 1.  IER) 

00  CONTINUE 
STOP 
END 


II 


i 


i 

« 
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I 
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THIS  ROUT  INC  RANDOM  I ZEB  A RANDOM  ARRAY 
BY  PERFORMINO  RANDOM  CIRCULAR  SHIFTS 
OF  THE  ROWS  AND  COLUMNS 


1 

I 

J 


1 


1 

I 

I 

i 


C 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 


40 


41 

50 


70 


71 

60 


900 

905 


to 

C 

C 

5 


C 


100 


SUBROUTINES  USED  RAND  RB 

IREAD  RB 

DIMENSION  I DD ( 256  > . KNAME ( 7 ) 

COMPLEX  lb ( 64. 64  > . TEMP ( 64 ) 

ACCEPT  " INPUT  SEED  FOR  RANDOM  NUMBER  " . NRND 

IMN-O 

I AT- 1 

h-6 

NR  -64 

NC  -64 

CALL  SDEF i IS. M.  IMN.  IAT. KNAME ) 

SHIFT  AND  INTERCHANGE  COLUMNS 

DO  50  150-1. NR 
CALL  RAND (NRND) 

INDE  X -FLOAT  (NRND)  *64  765536 
l NOE  X- INDEX + 32 

TVPF  INDEX. INDEX. INDEX 

DO  40  I 40- 1 . NC 

MOVE -1 ♦MOD ( 1 40+ INDEX.  64) 

1 FMP ( 1 40 ) - 1 S ( MOVF  150) 

CONTINUE 

DO  4 1 I 4 1 - 1 . NC 

I S < 14  1.  I 50) -TEMP ( 141 ) 

CONTINUE 

CONTINUE 

SHIM  AND  INI ERC MANGE  ROWS 

DO  60  160-1. NC 
(.ALL  RAND  (NRND) 

INDEX -MOAT  (NRND)  *64  / 65516 
1 NDE  X - I NDF  X ♦ □S 

TYPE  INDEX. INDEX. INDEX 

DO  70  I 70- 1 . NR 

MOVE- 1 ♦MOD ( I 70^ INDEX,  64) 

TEMP ( I70)-IS( 160.  MOVE) 

CONTINUE 
DO  71  l 7 1 - 1 . NR 

160-  I'M  ' - T E MP  ( 1M  ' 

CONI  I NVE 
LON) 1 NUE 


WRITE  FILE  TO  DISK 
WRI TE ( 10  900) 

FORMAT (M  NAME  OF  FILE  TO  BE  CREATED  - Z > 
READ ( 1 1 905 ) KNAME ( l ) 

L DRMAT(S?> 

CALL  DF  tl.W(  KNAME  ( t ).  IER) 

CALL  CFILW( KNAME ( 1 ). 3. IER) 

IF  (UR  NE  DSTQP  --CREATE  FILE  ERROR- 
CALL  OPEN ( 1 . KNAME ( 1 ) , 3.  IER.  513) 

IF (IER  NF  t>  STOP  - OPEN  ERROR  - 

NX  -64 

NY  -64 

NXC-1 

NYC-1 

JCMX-NX ►NXC-1 
JRMX-NN  * NYC  - l 

DO  10  II 0-1. 356 
l DD ( I 10 ) -0 
CONI INUk 
I DD ( 3 > -356 
I DD( 4 • -JRMX 

CALL  WRBLM  l.  0.  IDD.  1 . IER  ) 

IF (IER  NE  t)STOP  WRBLK  ERROR  - 
I DD<  3) -0 
I DD ( 4 ' -0 

DO  100  1R -NYC. JRMX 
I I R • l R *NYC  1 
DO  5 IC-l . NX 
IlC- IC^NXC  1 

I DD • 11C ) -RE Al  ( 18< 1C.  1R ) ) 

CONI INUF 

CALL  WRBl K( t IR  IDD.  1 . IFR) 

U ( IER  NE  1 ) ST  UP  WRULK  l RROR  - 
CONTINUE 

CALL  CLOSE  (1.  IER  ) 

STOP 

END 


> 

r 


I ! 

I1 


1 1 


r 

' 

! i 


ft  n o r>  n n 
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CRSCN  GENERATES  AN  256  X JRMX  ARRAY  CONTAINING 
AN  NX  X NY  IMAGE  WITH  UPPER  LH  CORNER  COORDS  (NXC.NYC) 

THE  PORTIONS  OF  THE  ARRAY  NOT  CONTAINING  THE  IMAGE  ARE 
FILLED  WITH  ZEROS  AND  THE  ARRAY  IS  STOREO  ON  DISK.  AS  KNAME 

OtMENilON  1 DO v 256 > . KNAME ( 7 > 

C 

C KNAME  IN  FILE  #1  IS  USED  TO  STORE  THE  BACKGROUND  IMAGE 

DO  80  130-1 . 16 
1 YF  t ‘ RuN  " . ItiO 

ACCCf'T  "COLUMN  ",NCOL,"ROW  “ , NROW,  "VALUE  “.VALUE 
WRITEUO.  900) 

900  FORMAT  ( " NAME  OF  FILE  TO  BE  CREATED  - Z > 

REAL i 1 1 , 905)  KNAME ( 1 ) 

905  FORMAT ( S7 ) 

■•At  L OF  U.W  ( KNAME  l 1 ) . IER  ) 

< At  I CFII  W l KNAME  ( 1 ),  2.  IER) 

IF  > IER  NE  t)  STOP  — CFILW  ERROR  (BCD)  — 

(ALL  OPEN(  1.  KNAME(  1 ).  2,  IER.  512) 

I F t IER  NE  1)  STOP  --OPEN  ERROR  (BCD)  -- 
C 

C 1NPU1  IMAGE  PARAMETERS 

C 

ACCEPT" IMAGE  SIZE  (NX, NY)  - " , NX . NY 
TypfcnoRDS  OF  UPPER  LH  CORNER  OF  IMAGE" 

TY PE ' E G (1.1)-  UPPER  LH  CORNER  OF  ENTIRE  ARRAY" 

ACCEPT  ■ (NXC. NYC ) - ".NXC.NYC 
,tCMX-NX*NXC  1 

tPM  » “N  t-NYT-  1 

if<  •(.("«  CT  256)  STOP — IMAGE  OFF  RIGHT  EDGE  OF  ARRAY  — 

F!“3 

APG-P  I/O 
bu  1 0 i ~ i ■ c.  Jo 
10  I DC ( l ) -0 

I uD ( 2 ) - ^56 
I DD  ( 4 ) -JRMX 

Call  HRBlal  1.  U.  IDD.  1.  IER) 

l F ( IER  NE  1 )STOP--WRBLK  ERR -- 

I DD ( 2 / -0 

IDD( 4 > -0 

NS  I -NYC-1 

IFlNSI  LE  0)  GO  TO  30 

00  20  IR-l.NSl 

20  (.ALL  WRBlK(1,  IR.  I DD  1,  IER) 

1 F t IER  NE  1)  STOP--WRBLK  ERR— 

30  CONTINUE 

RMtAN-0 
SIGMA-.  56 
NRaND-2351 
DO  50  l W -NYC .JRMX 
I IR- IR  «NYC- 1 
DO  5 IC-1 ■ NX 
IIC-IC+NXC-1 
IbD( I IC )=0 
CONTINUE 

I F ( I R FG  NROW) I OD(NCOL) -VALUE 
!F( IR  F Q NROWt IDD (NCOL^Tt -VALUE 

1 *•  t I R tO  NR'  IW+ J ) I t)D  ( NCOL  ) -VALUE 
IFilH  tO  NROW (3)1 DD ( NCOL  »3 ) “VALUE 
CALL  WRBLK( I. IR  IDD  1, IER) 

Irilt-R  Nfc  1'  STOP  — WRBLK  ERROR  ( BCD ) — 

50  CON  I I HUE 

CALL  CLOSE  < 1 IER  > 

Bu  luniinue 

,!Tor' 

END 
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C CR16SCN  FR  11/09/70 

C 

C CRSCN  GENERATES  AN  256  X JRMX  ARRAY  CONTAINING 

C AN  NX  A NY  IMAGE  WITH  UPPER  LH  CORNER  COORDS  (NXC.NYC) 

C THE  PORTIONS  OF  THE  ARRAY  NOT  CONTAINING  THE  IMAGE  ARE 

C FILLED  WITH  ZEROS  AND  THE  ARRAY  IS  STORED  ON  DISK  AS  KNAME 

C 

DIMENSION  I DD  * 256 ) • KNAME ( 7 ) 

C 

C KNAME  IN  FILE  #1  IS  USED  TO  STORE  THE  BACKGROUND  IMAGE 

DO  00  ISOM.  16 
TYPE  “RUN  100 

ACCEPT  'COLUMN  “.NCOL.  "ROW  ",  NROW.  "VALUE  ".VALUE 
WRITE* 10, 900) 

900  FORMAT ( " name  OF  FILE  TO  BE  CREATED  - " . Z ) 

READ* 1 i. 905/  KNAME ( 1 ) 

905  FORMAT <S7) 

(AIL  DF  I LW  ( KNAME  * t ) . I ER  ' 

CALL  CK ILW* KNAME* 1 ). 2. 1ER> 

I T < ! CR  NF  1)  STOP  — CFILW  ERROR  (BCD)  — 

CALL  OPENd.  KNAME*  1 ).  2.  IER.  512) 

IF  * ILR  NE  1)  STOP  --OPEN  ERROR  (BCD)  — 

C 

C INPUT  IMAGE  PARAMETERS 

C 

ACCEPT'  IMAGE  SIZE  *NX.NY)  - ",  NX,  NY 
TYRF "COORDS  OF  UPPER  LH  CORNER  OF  IMAGE" 

TYPE'E  G *1.1)-  UPPER  LH  CORNER  OF  ENTIRE  ARRAY" 

ACCEPT" (NXC. NYC ) - ".NXC.NYC 

jCMX»NX*NXC-l 

JRMX«NY*NvC-i 

1F*JCMX  OT  ?56 > STOP-- IMAGE  OFF  RIGHT  EDGE  OF  ARRAY -- 


P I *3  1415926 
ARG-PI/0 
DO  10  1-1.256 
10  I DD  * I ) -0 

IDO*  2 ) -256 
I DD  * 4 ) * JRMX 

CA|  t.  WRBI  K*  1. 0.  I DD,  1.  IER) 

IhdfcN  NE  1 JSTOP--WRBLK  ERR-- 
I DD  * 2 ) -0 
T DD  * 4 ) -0 
NSl-NYC-1 

IF  * NS  1 LE  0)  GO  TO  30 
DO  20  IR-1. NS 1 

20  CALL  WRBLK(1.  IR,  I DD.  1.  IER) 

TF«  IER  NE  1)  STOP- - WRBLK  ERR  — 
3u  CONTINUE 

RMEAN-0 
SIGMA-256 
NRAND-2351 
DO  50  IR-NYC , JRMX 
I IR«IR*NYC- 1 
DO  5 IC-1 . NX 
1 I C - IC ♦NXC- 1 
mo*  I IC  ) -0 
5 * ONTINUE 


IF*  1R 

EG 

NROW.'  I DD  * NCOL  > * VALUE 

IF  i IR 

EG 

NROW) TDD»NC0L*4) -VALUE 

i . 

IF'  IR 

EG 

NROW ) i DO  * NCOL  *8 ) -VALUE 

m IR 

EG 

NROW) I DD  * NCOL ♦ 12) -VALUE 

IF*  i R 

t*j 

NROW* 4 ) 1 DD  * NCOL > -VALUE 

IF*  IR 

EG 

NROW-* 4 > I DD*  NCOL  *4) -VALUE 

IF*  IR 

EG 

NROW-* 4 ) IDD*NCUL-*S) -VALUE 

If  * IR 

EG 

NROW* 4) I DD* NCOL* 12) -VALUE 

IF*  IR 

EG 

NROW  * 0 ) I DD  * NCOL ) -VALUE 

IE*  IR 

EG 

NROW*0 ) I DD  * NCOL *4 ) -VALUE 

IF  ( IR 

EG 

NROW*0) I DD* NCOL *0) -VALUE 

IF  ( IR 

EG 

NROW*  fl ) IDD* NCOL* 12) -VALUE 

IF*  IR 

EG 

NROW* 12) IDD* NCOL ) -VALUE 

IF*  IR 

EG 

NROW* 12) I DD  * NCOL ♦ 4 ) -VALUE 

IF(  IR 

EG 

NROW* 12) IDD (NCOL *0) -VALUE 

I] 

IF < IR 

EG 

NROW* 12) IDD (NCOL* 12) -VALUE 

CALL  WRBLK ( 1 - I R . I DD.  1,  IER) 

IF ( IER 

NE 

1)  STOP  --WRBLK  ERROR  (BGD) 

50  CONTINUE 


CALL  CLOSF ( 1, ICR ) 
90  CONTINUE 
STOP 
END 
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c 

c 

c 

c 

c 

c 

c 


c 

c 


CRLNOI.  FR 


11/10/78 


900 

905 


C 

C 

C 


10 


20 

30 


334 

333 


30 


6 
336 


80 


CRSCN  GENERATES  AN  236  X JRMX  ARRAY  CONTAINING 

AN  NX  X NY  IMAGE  WITH  UPPER  LH  CORNER  COORDS  (NXC.NYC). 

THE  PORTIONS  OF  THE  ARRAY  NOT  CONTAINING  THE  IMAGE  ARE 
FILLED  WITH  ZEROS  AND  THE  ARRAY  IS  STORED  ON  DISK  AS  KNAME 

DIMENSION  I DD ( 236 ) » KNAME ( 7 ) 

KNAME  IN  FILE  *1  IS  USED  TO  STORE  THE  BACKGROUND  IMAGE 

RMEAN-O. 

SIOMA-236. 

NRAND-2331 

ACCEPT  "RANDOM  SEED".  NR AND 
DO  80  180-1.  10 
TYPE  "RUN  ".180 

ACCEPT  "REFERENCE  SIZE  - ".NSIZE 
ACCEPT  "REFERENCE  POSITION  - «,NPOS 
WRITE ( 10. 900) 

FORMAT ("  NAME  OF  FILE  TO  BE  CREATED  - ".Z> 

READ(1 1.903)  KNAME ( l > 

FORMAT <S7> 

CALL  DF I LW ( KNAME ( 1 ) . I ER ) 

CALL  CFILW(KNAME ( 1 ) . 2.  IER ) 

IFdERNE.  1)  STOP  — CFILW  ERROR  (BCD)  — 

CALL  OPENd.  KNAMEd  ).  2.  IER.  312) 

IFdER  NE  1)  STOP  — OPEN  ERROR  (BOD)  — 

INPUT  IMAGE  PARAMETERS 

ACCEPT" IMAGE  SIZE  (NX, NY)  - ".NX. NY 
TYPE“COORDS  OF  UPPER  LH  CORNER  OF  IMAGE" 

TYPE"E. 0.  (1.1)-  UPPER  LH  CORNER  OF  ENTIRE  ARRAY" 

ACCEPT" (NXC,  NYC)  - ".NXC.NYC 

JCMX-NX+NXC-1 

JRMXaNY+NYC- 1 

IF( JCMX  . OT  236)  STOP — IMAGE  OFF  RIGHT  EDGE  OF  ARRAY — 

PI-3.  1413926 
ARO-PI/8 
DO  10  1-1.236 
IDD( I >-0 
IDD(2)— 236 
IDD(4)-JRMX 

CALL  WRBLKd.O.  IDD.  1.  IER) 

IF( IER  NE  1)ST0P — WRBLK  ERR — 

I DO ( 2 ) -O 
I DD ( 4 ) -O 
NS1— NYC-1 
IF(NS1  . LE 
DO  20  IR-1 


O) 
NS  1 


00  TO  30 


CALL  WRBLK ( 1.  IR.  IDD.  1.  ICR) 

IF( IER  . NE  1)  STOP — WRBLK  ERR — 
CONTINUE 

DO  30  IR— NYC, JRMX 
IIR— IR+NYC-1 

_ __  ^ 


IIR  QE  34)00  TO  333 


30 

31 

32 


OR 

AND 

AND 

AND 


IIC 
HC 
I IC 
IIC 


EQ. 

EQ 

EQ 

EQ 


IF ( HR  . LT.  28 
DO  3 IC— 1 , NX 
IIC-IC+NXC-1 
IFdIC  EQ  29 
IF( I IR  EQ 
IF( IIR  EQ 
IF( I IR  EQ 
CONTINUE 
CONTINUE 
CONTINUE 

CALL  WRBLK ( 1 . IR.  IDD.  1.  IER) 

IF(IER  NE  1)  STOP  — WRBLK  ERROR 

DO  6 16-1, NX 

IDD( 16) — O 

CONTINUE 

CONTINUE 

CONTINUE 

CALL  CLOSE (1,  IER) 

CONTINUE 

STOP 

END 


33 ) I DD ( 1 1 C ) - 1 00 

30)  IDDdIC  >-100 

31) 1 DD ( 1 1 C ) - 1 00 

32)  IDDdIC  >-100 


(BOD)— 
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C CPW  FR  12/20/70 

c 

l CRSCN  GENERATES  AN  256  1 JRMX  ARRAY  CONTAINING 

AN  N<  X NY  IMAGE.  WITH  UPPfcR  LH  CORNER  COORDS  (NXC.  NYC) 

C THE  PORTIONS  OF  THE  ARRAY  NOT  CONTAINING  THE  IMAGE  ARE 

C FILLED  WITH  ZEROS  AND  THE  ARRAY  IS  STORED  ON  DISK  AS  KNAME 

C 

DIMENSION  1DD'256>  KNAMEt7> 

C 

C KNAME  IN  FILE  •!  IS  OSED  TO  STORE  THE  BACKGROUND  IMAGE 

RMEAN-0 

SIGMA-256 
NR AND "2 J 5 1 

ACCEPT  -RANDOM  SEED"  NR AND 
DO  80  I0O-1. 10 

type  -run  ".ioo 

accept  -REFERENCE  SIZE  • "NSIZE 
accept  -REFERENCE  POSITION  - " . NPOS 
WRITE' 10  POO) 

900  FORMAT!"  NAME  OF  FILE  TO  BE  CREATED  - ",Z> 

R(  AD<  1 1 90S)  KNAME ( 1 ) 

905  FORMAT (S7> 

CALL  DF  ILW(  KNAME  ( 1 ) ■ IER  ) 

CALL  CF  ILW(  KNAME  < l ).  2.  IER) 

IF. IER  NE  1)  STOP  — CFIlW  ERROR  (BCD)  — 

CALL  OPEN' 1 . KNAME  l ).  2.  IER.  512) 

IF 'IER  NE  1)  STOP  --OPEN  ERROR  (BOD)  — 

L 

C INPUT  IMAGE  PARAMETERS 

C 

ACCEPT- IMAGE  SIZE  'NX. NY)  - "NX. NY 
TYPE  COORDS  OF  UPPER  LH  CORNER  OF  IMAGE" 

T /PE "E  G (11>-  UPPER  LH  CORNER  OF  ENTIRE  ARRAY" 

ACCEPT" (NXC. NYC ) - ".NXC.NYC 

jrMX-NX*NXC-l 

JRM*  -N i ♦NYC  - 1 

IMJ'N'  GT  256)  STOP-' IMAGE  OFF  RIGHT  EDGE  OF  ARRAY  — 


P I -3  1415926 
ARC-PI /8 
DO  10  1-1.256 

io  ioodf-0 

IDD'2>-256 
I DD  > 4 > -JRMX 

Call  wRBLKll.O.  I DD-  1.  IER) 

IF ' IER  NE  1 >STOP--WRBLK  ERR  — 

I DD<  2 > m0 
IDD<4>-0 
NSl-NYC-1 

IF(NS1  LE  0)  GO  TO  30 
DO  20  IR-1. NS  1 

20  CALL  WRBLKd,  IR.  I DD.  1.  IER  > 

IF  (IER  NE  1)  STOP--WRBLK  ERR  — 
30  CONTINUE 

DO  50  IR-NYC . JRMX 
I I R - IR+NYC- 1 


IF ' I IP 

LT 

29 

OR 

I IR 

GT 

35)00  TO  555 

DO  5 IC 

-1  ■ NX 

I IC-IO 

NXC-  1 

i 

I DD ( 32 ) 

-100 

IF'  I IR 

EO 

32 

AND 

IIC 

EO 

30) IDD< I IC )-lOO 

IF'  I IR 

EO 

32 

AND 

lie 

EO 

31  ) IDD(  1 10-100 

IF<  I IR 

EQ 

32 

AND 

IIC 

EO 

33)  IDD(  1 10-100 

IF*  I IR 

EO 

32 

AND 

IIC 

EQ 

34)  IDD(  110-100 

if  ( hr 

EQ 

32 

AND 

I IC 

EQ 

29)  I DD  ( 1 10-100 

IF ( I IR 

EO 

32 

AND 

IIC 

EQ 

35)  IDD(  1 10-100 

IF<  I IR 

EO 

29 

AND 

IIC 

EO 

49) IDD( I IC )-l00 

IF'  I IR 

EG 

29 

AND 

IIC 

EQ 

50)  IDD(  1 10-100 

IF'  I IR 

EO 

29 

AND 

IIC 

EQ 

S1>IDD(IIC)-100 

IF(  I IR 

EQ 

30 

AND 

IIC 

EQ 

51  ) IDD(  1 10-100 

IF(  I IR 

EO 

31 

AND 

IIC 

EO 

51  ) I DD ( 1 10-100 

IF ( I IR 

EO 

32 

AND 

IIC 

EQ 

49  ) I DD ( 110-100 

IF  t I IR 

EO 

32 

AND 

IIC 

EO 

50)  I DD(  1 10-100 

IF'  iir 

EO 

32 

AND 

IIC 

EQ 

51  ) IDD(  110-100 

I F ( I IR 

EO 

33 

AND 

IIC 

EO 

51  ) IDD(  1 10-100 

IF'  I IR 

EO 

34 

AND 

IIC 

EQ 

51  ) IDD(  I 10-100 

IF'  I IR 

EO 

35 

AND 

IIC 

EQ 

49  ) I DD(  110-100 

IF  « IIR 

EO 

35 

AND 

lie 

EQ 

50)  IDD(  I 10-100 

IF . [ IR 

EO 

35 

AND 

ire 

EO 

51 ) IDD< IIC  >«100 

C 

554  CONTINUE 

5 CONTINUE 

555  CONTINUE 

CALL  WRBLKd.  IR  I DD-  1 IER) 

IF ( IER  NE  1)  STOP  --WRBLK  ERROR  (BCD) — 
DO  6 1 6-1. NX 
I DDi  1 6 > -0 
6 CONTINUE 

556  CONTINUE 

50  CONTINUE 

CALL  CLOSE ( 1. IER) 

00  CONTINUE 
STOP 
END 
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C CR3  12/20/78 

C 

C CRSCN  GENERATES  AN  236  X JRMX  ARRAY  CONTAINING 

C AN  NX  X NY  IMAGE  WITH  UPPER  LH  CORNER  COORDS  (NXC.NYC) 

C THE  PORTIONS  OF  THE  ARRAY  NOT  CONTAINING  THE  IMAGE  ARE 

C FILLED  WITH  ZEROS  AND  THE  ARRAY  IS  STORED  ON  DISK  AS  KNAME 

C 

DIMENSION  I DD ( 236  > ‘ KNAME  < 7 ) 

C 

C KNAME  IN  FILE  »]  IS  USED  TO  STORE  THE  BACKGROUND  IMAGE 

RMEAN-0 

SIGMA-236 

NRAND-2331 

ACCEPT  ''RANDOM  SEED",  NR  AND 
DO  80  180-1. 10 
TYPE  "RUN  ",  ISO 

ACCEPT  "REFERENCE  SIZE  - ".NSIZE 
ACCEPT  "REFERENCE  POSITION  - ",NPOS 
WRITE! 10, 900) 

900  FORMAT!"  NAME  OF  FILE  TO  BE  CREATED  - ",Z> 

READ! 11. 903)  KNAME  < 1 ) 

903  FORMAT ( S7 > 

CALL  DF ILW( KNAME ( 1 ),  IER) 

CALL  CF I LW ( KNAME ( 1 ) , 2.  IER) 

IF! IER  NE  1)  STOP  — CFILW  ERROR  (BCD)  — 

CALL  OPEN!  1.  KNAME!  1 ),  2.  IER,  312) 

IF! IER  NE  1)  STOP  — OPEN  ERROR  ( BCD ) — 

C 

C INPUT  IMACE  PARAMETERS 

C 

ACCEPT" IMAGE  SIZE  (NX, NY)  - ".NX, NY 
TYPE"CDORDS  OF  UPPER  LH  CORNER  OF  IMAGE" 

TYPE"E  0 (1,1)-  UPPER  LH  CORNER  OF  ENTIRE  ARRAY" 

ACCEPT" (NXC.  NYC ) - ".NXC.NYC 

JCMX-NX+NXC-1 

JRMX-NY+NYC-l 

IF1JCMX  OT  236)  STOP  — IMAGE  OFF  RIGHT  EDGE  OF  ARRAY  — 

PI -3  1413926 
ARG-PI/8 
DO  10  1-1,236 
10  I DD ! I ) -0 

I DD  < 2 ) “236 
I DD ( 4 ) » JRMX 

CALL  WRBLK!  1, 0.  IDD,  1.  IER) 

IF!  IER  NE  DSTOP— WRBLK  ERR  — 

I DD ( 2 ) -0 
IDD ( 4 ) -0 
NS1-NYC-1 

IF! NS  1 LE  O)  GO  TO  30 
DO  20  IR-1, NS1 

20  CALL  WRBLK  ( 1 , IR.  IDD,  1.  IER) 

IF! IER  NE  1)  STOP — WRBLK  ERR — 

30  CONTINUE 

DO  30  IR-NYC, JRMX 
I IR  = IR+NYC- 1 

IF ( I IR  LT  29  OR  IIR  GT  33)G0  TO  333 
DO  3 IC-1,  NX 
1 1C- IC  + NXC-1 
IDD! 32>-100 


IF! I IR 

EQ 

32 

AND 

I IC 

EQ 

30) IDD! I IC > = 100 

IF! I IR 

EQ 

32 

AND 

I IC 

EQ 

31 ) IDD! I IC >=100 

IF! I IR 

EQ 

32 

AND 

IIC 

EQ 

33) IDD! I IC  > = 100 

IF! I IR 

EQ 

32 

AND 

I IC 

EQ 

34)  IDD!  I 10  = 100 

IF! I IR 

EQ 

32 

AND 

I IC 

EQ 

29)  IDD!  1 10  = 100 

IF!  I IR 

EQ 

32 

AND 

IIC 

EQ 

33) IDD! I IC  >-100 

IF!  IIR 

EQ 

29 

AND 

I IC 

EQ 

49) IDD! I IC  > = 100 

IF ( I IR 

EQ 

29 

AND 

IIC 

EQ 

30)  IDD!  I 10  = 100 

IF!  IIR 

EQ 

29 

AND 

I IC 

EQ 

31  ) IDD!  I 10-100 

IF! I IR 

EQ 

30 

AND 

I IC 

EQ 

31 ) IDD! IIC  >-100 

IF! IIR 

EQ 

31 

AND 

IIC 

EQ 

31  HDD!  I IC  > = 100 

IF! I IR 

EQ 

32 

AND 

I IC 

EQ 

49)  IDD!  1 10-100 

IF! I IR 

EQ 

32 

AND 

I IC 

EQ 

30) IDD! 1 IC  > = 100 

IF! I IR 

EQ 

32 

AND 

I IC 

EQ 

31 ) IDD! 1 IC  >-100 

IF!  I IR 

EQ 

33 

AND 

I IC 

EQ 

31 > IDD! I IC  > = 100 

IF!  I IR 

EO 

34 

AND 

IIC 

EQ 

31  ) IDD!  I 10  = 100 

IF! I IR 

EQ 

3* 

AND 

1 IC 

EQ 

49) IDD! I IC  > = 100 

IF! I IR 

EQ 

3 3 

AND 

I IC 

EQ 

30) IDD!  I IC  > = 100 

IF! I IR 

EQ 

33 

AND 

I IC 

EQ 

31 ) IDD! 1 1C  1 = 100 

C 

334  CONTINUE 
3 CONTINUE 

333  CONTINUE 

CALL  WRBLK!  1.  IR,  IDD.  1.  IER) 

IF! IER  NE  1)  STOP  — WRBLK  ERROR  (BCD) — 
DO  6 16-1, NX 
IDD! 16 ) -O 
6 CONTINUE 
336  CONTINUE 
30  CONTINUE 

CALL  CLOSE! 1 . IER ) 

80  CONTINUE 
STOP 


C CR4  12/20/78 

C 

C CRSCN  GENERATES  AN  2S6  X JRMX  ARRAY  CONTAINING 

C AN  NX  X NY  IMAGE  WITH  UPPER  LH  CORNER  COORDS  (NXC.NYC) 

C THE  PORTIONS  OF  THE  ARRAY  NOT  CONTAINING  THE  IMAGE  ARE 

C FILLED  WITH  ZEROS  AND  THE  ARRAY  IS  STORED  ON  DISK  AS  KNAME 

C 

DIMENSION  100(236 >. KNAME< 7) 


C 

C KNAME  IN  FILE  »1  IS  USED  TO  STORE  THE  BACKGROUND  I MACE 

RMEAN-0 

S10MA=236 

NRAND-2331 

ACCEPT  "RANDOM  SEED" , NR AND 
DO  80  100*1. 10 
TYPE  "RUN  ISO 

ACCEPT  "REFERENCE  SIZE  * ".NSIZE 
ACCEPT  "REFERENCE  POSITION  - ".NPOS 
WR I TE ( 10, 900) 

900  FORMAT ( " NAME  OF  FILE  TO  BE  CREATED  = ",Z> 

READ  <11, 903 ) KNAME ( 1 ) 

903  FORMAT <S7> 

CALL  DF I LW ( KNAME ( 1 ) , I ER ) 

CALL  C F I LW ( KNAME ( 1 ) , 2 , I ER ) 

I F ( I ER  NE  1)  STOP  — CFILW  ERROR  ( BCD ) — 

CALL  OPEN < 1 , KNAME  < 1 ) , 2,  I ER,  312) 

IF ( IER  NE  1)  STOP  — OPEN  ERROR  (BCD)  — 

C 

C INPUT  IMAGE  PARAMETERS 

C 

ACCEPT"IMACE  SIZE  (NX, NY)  » ",NX,NY 
TYPE "COORDS  OF  UPPER  LH  CORNER  OF  IMAGE" 

TYPE " E G (1,1)*  UPPER  LH  CORNER  OF  ENTIRE  ARRAY" 

ACCEPT" (NXC, NYC ) * ",NXC,NYC 

JCMX-NX+NXC-1 

JRMX*NY+NYC-1 

IF ( JCMX  GT  236)  STOP — IMACE  OFF  RIGHT  EDGE  OF  ARRAY — 

PI *3  1415926 
ARG=P 1/8 
DO  10  1*1 , 236 
10  IDD( I ) *0 

I DD ( 2 ) =256 
I DD ( 4 ) * JRMX 

CALL  WRBLK (1,0, I DD, 1, IER) 

IF ( IER  NE  1 >STOP--WRBLK  ERR — 

IDD( 2) *0 
I DD ( 4 ) *0 
NS1*NYC-1 

IF ( NS 1 LE  0)  GO  TO  30 
DO  20  IR*1, NS 1 

20  CALL  WRBLK ( 1 , IR, I DD, I, IER) 

IF ( IER  NE  D STOP--WRBLK  ERR  — 

30  CONTINUE 

DO  30  IR*NYC, JRMX 
I I R = I R +NYC  - 1 


I F ( I I R 

LT 

29 

OR 

I IR 

GT 

33 ) GO  TO  333 

DO  3 IC 

*1 . NX 

I IC • IC +NXC- 

1 

IDD (32) 

*100 

IF(  I IR 

EQ 

32 

AND 

I IC 

EQ 

30) IDD( I IC  >*100 

IF  < I IR 

EQ 

32 

AND 

I IC 

EQ 

31  ) I DD ( IIC  > = 100 

I F ( I I R 

EQ 

32 

AND 

I IC 

EQ 

33) I DD ( I IC ) = 100 

IF<  I IR 

EQ 

32 

AND 

I IC 

EQ 

34  ) IDD( IIC  > = 100 

IF(  I IR 

EQ 

32 

AND 

I IC 

EQ 

29) IDD < IIC ) *100 

IF(  I IR 

EQ 

32 

AND 

I IC 

EQ 

33) IDD( I IC  > = 100 

IF(  I IR 

EQ 

29 

AND 

IIC 

EQ 

49) IDD< I IC ) =100 

IF < I IR 

EQ 

30 

AND 

IIC 

EQ 

49) IDD( I IC  > = 100 

IF(  I IR 

EQ 

31 

AND 

l IC 

EQ 

49) IDD ( I 1C ) = 100 

IF< I IR 

EQ 

32 

AND 

I IC 

EQ 

49 ) IDD( I IC ) =100 

IF ( I IR 

EQ 

32 

AND 

I IC 

EQ 

50)  I DD < I IC  > = 100 

IDD<  51  >*100 

C 

334  CONTINUE 
3 CONTINUE 

333  CONTINUE 

CALL  WRBLK  < 1 , IR  IDD.  1.  IER) 

IF< IER  NE  1)  STOP  --WRBLK  ERROR  ( BGD ' — 
DO  6 16*1 , NX 
IDD( 16 ) *0 
6 CONTINUE 
336  CONTINUE 
30  CONTINUE 

CALL  CLOSE ( 1, IER ) 

00  CONTINUE 
STOP 
END 
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C CR5  12/20/78 

C 

C CRSCN  GENERATES  AN  236  X JRMX  ARRAY  CONTAINING 

C AN  NX  X NY  IMAGE  WITH  UPPER  LH  CORNER  COORDS  (NXC.NYO 

C THE  PORTIONS  OP  THE  ARRAY  NOT  CONTAINING  THE  IMAGE  ARE 

C FILLED  WITH  ZEROS  AND  THE  ARRAY  IS  STORED  ON  DISK  AS  KNAME 

C 

DIMENSION  I DD ( 236 ) . KNAME ( 7 ) 

C 

C A NAME  IN  FILE  *1  IS  USED  TO  STORE  THE  BACKGROUND  IMAGE 

RMEAN*0 

SIGMA«236 

NRAN0»2331 

ACCEPT  "RANDOM  SEED". NR AND 
DO  80  180-1. 10 
TYPE  "RUN  ",  180 

ACCEPT  "REFERENCE  SIZE  = ".NSIZE 
ACCEPT  "REFERENCE  POSITION  - ".NP03 
WRITE! 10. 900) 

900  FORMAT ( " NAME  OF  FILE  TO  BE  CREATED  = " . Z ) 

READ  <11. 903 ) KNAME  < 1 ) 

903  FORMAT (S7> 

CALL  DF ILW ( KNAME ( 1 ) . IER ) 

CALL  CF I LW ( KNAME ( 1 ) . 2 IER ) 

IF (IER  NE  1)  STOP  — CFILW  ERROR  (BCD)  — 

CALL  OPENl 1, KNAMEC 1 >. 2. IER. 312) 

IF ( I ER  NE  1)  STOP  — OPEN  ERROR  (BCD)  — 

C 

C INPUT  IMAGE  PARAMETERS 

C 

ACCEPT"IMACE  SIZE  (NX. NY)  - ".NX. NY 
TYPE  COORDS  OF  UPPER  LH  CORNER  OF  IMAGE" 

TYPE'E  G (1.1)=  UPPER  LH  CORNER  OF  ENTIRE  ARRAY" 

ACCEPT" (NXC. NYC > » "NXC.NYC 

JCMX-NX+NXC-1 

JRMX=NY*NYC“t 

I F ( JCMX  CT  236)  STOP — I MAGE  OFF  RIGHT  EDGE  OF  ARRAY-- 

PI  =3  1415926 
ARG»P I / 8 
DO  10  1=1,236 
10  I DD ( I ) =0 

I DD  < 2 ) =236 
I DD ( 4 ) * JRMX 

CALL  WRBLK1 1. 0, IDD. 1. IER ) 

IF (IER  NE  1 >ST0P--WR8LK  ERR — 

I DD ( 2 ) >0 
IDD( 4 ) =0 
NS1-NYC-1 

IF (NS1  LE  0)  CO  TO  30 
DO  20  IR-I. NS 1 

20  CALL  WRBLK ( 1 . IR,  IDD.  1.  IER ) 

I F ( I ER  NE  1)  STOP — WRBLK  ERR — 

30  CONTINUE 

DO  30  IR-NYC.JRMX 
I IR=IR*NYC-I 


IF( I IR 

LT 

29 

OR 

I IR 

GT 

33 ) GO  TO  333 

DO  3 IC 

■1.  NX 

I IC-IONXC- 

1 

IDD ( 32 ) * 100 

IF( I IR 

EQ 

32 

AND 

I IC 

EQ 

30)  I DD ( I 10  = 100 

IF(  I IR 

EQ 

32 

AND 

I IC 

EQ 

31 > IDD< I IC )=100 

IF ( I IR 

EQ 

32 

AND 

I tc 

EQ 

33) IDD( I IC )=100 

IF<  I IR 

EQ 

32 

AND 

I IC 

EQ 

34 ) IDD( IIC )=100 

I F ( I I R 

EQ 

32 

AND 

I IC 

EQ 

29) IDD<  I IC  > = 100 

IF(  I IR 

EO 

32 

AND 

I IC 

EQ 

33) I DD ( I IC  > = 100 

IF(  l IR 

EQ 

29 

AND 

1 1C 

EQ 

49) IDD<  I IC  >*100 

IF(  I IR 

EQ 

29 

AND 

I IC 

EQ 

30) IDD< I IC  > = 100 

IF(  I IR 

EQ 

29 

AND 

I IC 

EQ 

31 > IDO< I IC )=100 

IF(  I IR 

EQ 

30 

AND 

1 IC 

EQ 

49) IDD( I IC )=100 

IF(  I IR 

EQ 

31 

AND 

1 1C 

EQ 

49)  I DD ( I 10  = 100 

IF( I IR 

EQ 

32 

AND 

I IC 

EQ 

49) I DD ( I IC 1=100 

IF(  I IR 

EQ 

32 

AND 

I IC 

EQ 

SO) IDD(  I IC  > = 100 

IF(  I IR 

EQ 

32 

AND 

I IC 

EQ 

31  ) I DD  ( 110  = 100 

IF(  I IR 

EQ 

33 

AND 

I IC 

EQ 

31  ) I DD  ( 1 10=100 

IF< I IR 

EQ 

34 

AND 

I IC 

EQ 

31  ) I DD ( I IC  > = 100 

IF(  I IR 

EQ 

33 

AND 

IlC 

EQ 

49)  IDD ( 1 10  = 100 

if t iir 

EQ 

33 

AND 

IIC 

EQ 

50  > I DD  ( 110  = 100 

IF(  I IR 

EQ 

33 

AND 

IIC 

EO 

31  ) IDD ( 110=100 

C 

334  CONTINUE 
3 CONTINUE 

333  CONTINUE 

CALL  WRBLK ( 1 . IR. IDD  1. IER) 

IF ( IER  NE  1)  STOP  — WRBLK  ERROR  ( BGD  > -- 
DO  6 16*1 . NX 
I DD ( 16 ) *0 
6 CONTINUE 
336  CONTINUE 
30  CONTINUE 

CALL  CLOSE ( 1 . IER ) 

80  CONTINUE 
STOP 
END 
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I 

I 

1 

I 

[ 


I 
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C CRb  12/20/78 

C 

C CRSCN  GENERATES  AN  236  X JRMX  ARRAY  CONTAINING 

C AN  NX  X NY  IMAGE  WITH  UPPER  LH  CORNER  COORDS  (NXC.NYC) 

C THE  PORTIONS  OF  THE  ARRAY  NOT  CONTAINING  THE  IMAGE  ARE 

C FILLED  WITH  ZEROS  AND  THE  ARRAY  IS  STORED  ON  DISK  AS  KNAME 

C 

DIMENSION  I DD < 236 ) . KNAME ( 7 ) 


C 

C KNAME  IN  FILE  #1  IS  USED  TO  STORE  THE  BACKGROUND  IMAGE 

RMEAN-0 

SIGMA-236 

NRAND-2351 

ACCEPT  "RANDOM  SEED" . NR AND 
DO  80  180-1. 10 
TYPE  “RUN  ". 180 

ACCEPT  "REFERENCE  SIZE  - ".NSIZE 
ACCEPT  "REFERENCE  POSITION  » ",NPOS 
WRITE! 10,  900) 

900  FORMAT ( " NAME  OF  FILE  TO  BE  CREATED  - ",Z> 

REA0(1 1.903)  KNAME! 1) 

903  FORMAT !S7> 

CALL  DFILWIKNAME! 1 ) , IER ) 

CALL  CF ILW 1 KNAME ( 1 ) . 2.  IER) 

IF! IER  NE  1)  STOP  — CF I LW  ERROR  IBGD)  — 

CALL  OPEN! 1, KNAME! 1 ).  2.  IER,  312) 

IF! IER  NE  1)  STOP  — OPEN  ERROR  (BGD)  — 

C 

C INPUT  IMAGE  PARAMETERS 

C 

AC C EP T " I MAGE  SIZE  1NX.NY)  - ".NX, NY 
TYPE"COORDS  OF  UPPER  LH  CORNER  OF  IMAGE" 

TYPE "E  C (1,1)-  UPPER  LH  CORNER  OF  ENTIRE  ARRAY" 
ACCEPT" (NXC, NYC ) - ".NXC.NYC 
JCMX-NX+NXC-t 
JRMX— NY+NYC— 1 

IFfJCMX  OT  236)  STOP — IMAGE  OFF  RIGHT  EDGE  OF  ARRAY — 

PI-3  1415926 
ARG-PI/8 
DO  10  I - 1 . 236 
10  IDD( I ) -0 

I DD! 2 ) -236 
I DD ( 4 ) - JRMX 

CALL  WRBLK ( 1 . 0.  I DD,  1.  IER) 

IF!  IER  NE  1 (STOP — WRBLK  ERR— 

I DD! 2 ) -O 
I DD ( 4 ) -0 
NS 1 -NYC- 1 

IF ! NS  1 LE  0)  GO  TO  30 
DO  20  IR-1 . NS1 

20  CALL  WRBLK! 1,  IR,  IDD,  1.  IER) 

IF!  IER  NE  1)  STOP  — WRBLK  ERR  — 

30  CONTINUE 

DO  30  IR— NYC , JRMX 
I IR-IR+NYC-1 

IF!  I IR  LT  29  OR  HR  GT  33)G0  TO  333 
DO  3 IC-1 , NX 
I IC-IC+NXC-1 
I DD ! 32 ) - 1 00 


IF! I IR 

EO 

32 

AND 

I 1C 

EQ 

30) IDD! I IC ) -100 

IF!  I IR 

EO 

32 

AND 

I IC 

EQ 

31 > IDD! I IC >-100 

IF! I IR 

EG 

32 

AND 

I 1C 

EQ 

33) IDD! I IC  >-100 

IF!  I IR 

EG 

32 

AND 

I IC 

EQ 

34) IDD! I IC ) -100 

IF!  I IR 

EQ 

32 

AND 

IIC 

EQ 

29) IDD! I IC ) -100 

IF! I IR 

EO 

32 

AND 

I IC 

EQ 

33)  IDD!  I 10-100 

IF! I IR 

EQ 

29 

AND 

I IC 

EQ 

49) IDD! I IC )-lOO 

IF! I IR 

EQ 

29 

AND 

IIC 

EQ 

30) IDD! I IC  >-100 

IF! I IR 

EQ 

29 

AND 

IIC 

EQ 

31 > IDD! 1 1C  >-100 

IF!  I IR 

EQ 

30 

AND 

I IC 

EQ 

49) IDD! 1 1C ) -100 

IF!  I IR 

EQ 

31 

AND 

I IC 

EQ 

49) IDD! I IC ) -100 

IF!  I IR 

EQ 

32 

AND 

I IC 

EQ 

49) IDD! I IC )-100 

IF!  I IR 

EQ 

32 

AND 

I IC 

EQ 

30) IDD! I IC ) -100 

IF!  I IR 

EQ 

32 

AND 

1 IC 

EQ 

31 ) IDD! 1 IC  >-100 

IF ( I IR 

EQ 

33 

AND 

1 IC 

EQ 

31 ) IDD! I IC  >-100 

IF!  I IR 

EQ 

33 

AND 

I IC 

EQ 

49) IDD! I IC  >-100 

IF!  I IR 

EQ 

34 

AND 

IIC 

EQ 

31 > IDD! I IC >-100 

IF!  I IR 

EQ 

34 

AND 

I IC 

EQ 

49) IDD! I IC  >-100 

IF! I IR 

EQ 

33 

AND 

I IC 

EQ 

49) IDD!  1 IC  >-100 

IF!  I IR 

EQ 

33 

AND 

I IC 

EQ 

30) IDD! I IC  > = 100 

IF! I IR 

EQ 

33 

AND 

1 1C 

EQ 

31 ) IDD! I IC )-100 

C 

334  CONTINUE 
3 CONTINUE 

333  CONTINUE 

CALL  WRBLK ! 1 , IR. IDO. 1. IER) 

IF! IER  NE  1)  STOP  — WRBLK  ERROR  (3GD) — 
DO  6 16-1 . NX 
IDD! 16 ) -0 
6 CONTINUE 
336  CONTINUE 
30  CONTINUE 

CALL  CLOSE (1, IER) 

80  CONTINUE 
STOP 
END 
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U/10/7B 


i 


C CR6Q  FR 

C 

C CRbCN  GENERATES  AN  256  X URMX  ARRAV  CONTAINING 

C AN  NX  X NV  IMAGE  WITH  UPPER  LH  CORNER  COORDS  (NXC.NVC) 

C THE  PORTIONS  OF  THE  ARRAV  NOT  CONTAINING  THE  IMAGE  ARE 

C FILLED  WITH  ZEROS  AND  THE  ARRAV  IS  STORED  ON  DISK  AS  KNAME 

C 

DIMENSION  IDD!256>.  KNAMEI7) 

C 

C KNAME  IN  FILE  #1  IS  USED  TO  STORE  THE  BACKGROUND  IMAGE 


RMEAN-0 

SICMA-256 

NRAND-2351 

ACCEPT  "RANOOM  SEED".  NR  AND 
DO  80  180-1. 10 
TvPF  "RUN  ISO 

ACCEPT  "REFERENCE  SIZE  - ".  NSI ZE,  ISI ZE 
ACCEPT  "REFERENCE  POSITION  - " , NPOS 
WRITE!  10.  900) 

■700  FORMAT!"  NAME  OF  FILE  TO  BE  CREATED  - ",  Z > 

READ! I 1 . 905)  KNAME! 1 ) 

905  FORMAT (S7> 

CALL  DFILWIKNAME! 1 ).  IER) 

CALL  CFILWIKNAME!  1),  2.  IER) 

IF ! IER  NE  l)  STOP  — CFILW  ERROR  IBCD)  — 

CALL  OPEN! 1, KNAME! 1 ).  2.  IER.  512) 

IF!  IER  NE  I)  STOP  —OPEN  ERROR  ( BGD ) — 

C 

C INPUT  IMAGE  PARAMETERS 

C 

ACCEPT"IMAGE  SIZE  INX,NV)  - " , NX . NY 
TVPE"COQRDS  OF  UPPER  LH  CORNER  OF  IMAGE" 

TVPE"E  G <1.  1)-  UPPER  LH  CORNER  OF  ENTIRE  ARRAV" 
ACCEPT'MNXC.NYC)  - ".NXC.NVC 
JCMX— NX+NXC— 1 
JRMX-NY+NYC- 1 

IF  ( JCMX  GT  256)  STOP  — IMAGE  OFF  RIGHT  EDGE  OF  ARRAY- 

PI -3  1415926 
ARG-PI/8 
DO  10  1-1.256 
10  IDD< I >-0 

IDD<2)-256 
IDD ! 4 ) -URMX 

CALL  WRBLK!  1,0.  IDD.  1.  IER) 

IF!  IER  NE  DSTOP  — WRBLK  ERR  — 

IDD 12 ) -0 
I DDl 4 ) -0 
NSI-NVC-I 

IF1NS1  LE  0)  GO  TO  30 
DO  20  IR-1.  NSI 

20  CALL  WRBLK!  I.  IR.  IDD.  1.  IER  ) 

IF!  IER  NE  1>  STOP — WRBLK  ERR  — 

30  CONTINUE 

DO  50  IR-NVC.JRMX 
I IR-IR+NYC-1 


IF! IIR  LT  NPOS 

OR 

IIR 

GE 

NPOS*NS I ZE ) GO 

TO 

555 

DO  5 IC-1,  NX 

I I C-IC+NXC- 1 

300  IF! 1IC  LT  NPOS 

OR 

IIC 

CE 

NPOS* ISI ZE ) GO 

TO 

554 

IDD!  1 10-100 

554  CONTINUE 

5 CONTINUE 

555  CONTINUE 

CALL  WRBLK!  1.  IR.  IDD.  1.  IER) 

IF! IER  NE  1)  STOP  --WRBLK  ERROR  ! BGD ) — 

IF! I IR  LT  NPOS  OR  IIR  OE  NP0S-NS1 ZE >C0  TO  556 
00  6 1 6- 1 . NX 
IDD! 161*0 
6 CONTINUE 
55e  CONTINUE 
50  CONTINUE 

CALL  ClOSE! t IER ) 

80  CONTINUE 
STOP 
END 


n n n n n 


IMAGER  FR  11/3/78 

THIS  ROUTINE  READS  AN  INTEGER  IMAGE 
AND  OUTPUTS  A GRAYSCALE  PLOT 

COMPLEX  IS  < 64. 64  > 

DIMENSION  KNAME ( 7 ) 

NR-64 

NC-64 

M=6 

IAT-1 

IMN-1 

CALL  SDEF< IS.  M,  IMN.  IAT . KNAME ( 1 ) ) 
CALL  QSPLT< IS,  NC.  NR, 3) 

STOP 

END 
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IREAD.FK  09/01  / 78 
C 

SUBROUTINE  SD£F < IS. N,  INN,  IAT.  KNAME ) 

C 

C SDEF  READS  A PORTION  OF  A REAL  IMAOE  STORED  ON  OISK 

C AND  STORES  IT  AS  THE  REAL  COMPONENT  OF  COMPLEX  ARRAY  IS 

C 

DIMENSION  KII256.  5>.KNAME(7> 

COMPLEX  IS(64.  64).  CZERO 

C 

IMN-IMN*1 

C 

NC-2«»M 

NR-NC 

HM-N 

NCC-2«»MM 

NRR-NCC 

C 

WRITE) 10. 900) 

900  FORMAT  < " INPUT  FILENAME  • ",Z> 

READ) 1 1 . 905 ) KNAME ( 1 ) 

905  FORMAT (S7> 

CALL  OPEN)  1.  KNAME)  1 ).  0.  IER.  512) 

IFHER  NE  1 (STOP--OPEN  ERR  — 

910  FORMAT)"  INPUT  FILENAME  “,S10> 

CZER0-CMPLX10  . 0 ) 

DO  17  15-1.  NR 

DO  17  16-1,  NC 

17  IS) 16.  1 5 ) -C  ZERO 

CALL  RD8LK!  1. 0.  HI.  1.  IER) 

IF) IER  NE  1 (STOP — RDBL  ERR-- 
NCI-KI12. 1) 

NR  I -K  1)4,  1 ) 

TYPE  "NCI.  NR  I".  NCI.  NR  I 

TYPE" INPUT  ! X . Y > COORDS  OF  UPPER  LH  CORNER  OF  IMAOE" 
TYPE“E  0 t 1 , 1 ) - UPPER  LH  CORNER  OF  IMAOE  ON  DISK" 
ACCEPT"  - ".NXC.NYC 

920  FORMAT)"  UPPER  LH  COORDS  IX,  Y)  -".215) 

JCMX-NXC*NCC-1 

JRMX-NYC*NRR-1 

IF ) JCMX  OT  NCI)  STOP--NXC  TOO  LARGE 

IF ) URHX  OT  NR  I ) STOP--NYC  TOO  LARGE 

1 1 -1 

IR  — 4 

NBLOCK-S 

LOOP-NRR / NBLOCK 

lob-nrr->loop»nblock> 

IF ( LOOP  EQ  0)00  TO  1001 

DO  40  140-1.  LOOP 

IR-) 140-1 >«5*NYC 

CALL  RDBLK11.  IR.  KI.  NBLOCK,  IER) 

IF)  IER  NE  DSTOP--RDBL  ERR- - 
1002  CONTINUE 
JJ-0 

DO  42  142-1. NCC 

IC»I42*NXC-1 

JJ-JJ-M 

DO  43  143-1. NBLOCK 
1 143-143-1 

SS-FL0AT1KI 1 IC. 143) ) 

IF ) SS  OT  1023)SS-1023 
IS1UJ. II*II43)-CMPLX(SS.O  > 

43  CONTINUE 
42  CONTINUE 
11-11*5 

1001  CONTINUE 
40  CONTINUE 

IF1L0B  EO  0)00  TO  50 

NBLQCK-LOB 

LOB-O 

IR-IR*5 

CALL  RDBLKl  1,  IR.  KI.  NBLOCK,  IER  ) 

IF) IER  NE  1 ) STOP  --RDBL  ERR- 
OO  TO  1002 
50  CONTINUE 

TYPE  IS)  1.  1 ).  IS (64.  64).  IS)  10.  2) . IS (2.  10) 

CALL  CLOSE) 1,  IER) 

RETURN 

END 


i 

i ; 

. n 


P 


\i 


' 
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4.2.2  HDRV.FR 


IHDRV.FR  is  a FORTRAN  program  which  creates  a single  hologram  and  performs 
an  image  reconstruction.  The  holograph  is  created  by  multiplying  the  two- 
dimensional  Fourier  transform  of  an  image  by  the  complex  conjugate  of  the 
Fourier  transform  of  a reference  image.  The  multiplication  is  done  on  a 
point-to-point  basis.  The  reconstruction  is  accomplished  by  multiplying  the 
hologram  by  the  Fourier  transform  of  the  reference  on  a point-by-point  basis, 
and  taking  the  inverse  Fourier  transform.  The  reconstructed  image  is  written 
to  a disk  file,  in  complex  format. 

I « 

1 | * 

1 1 * 


10/ SO/ 78 


REV  A 


10/31/70 


C HPRV  FR 
C 

C THIS  IS  A DRIVER  PROGRAM  TO  GENERATE  HOLOCRAM 

C DATA  FILES 

c 

c 

Dimension  kname* 7>.  jrf.f*  7), ksum*  7> , kfname* 7 >, jfref < 7 > 
COMPLEX  16*64, 64),  TEMP  1 64) 

IMN-l 

M-6 

1AT-1 
NR»2*»M 
NC  -NR 

TVPE  "SUDJECT  INPUT- 

FALL  SDEF  ( IS,  M,  1MN,  I AT . KNAMF  l 1 ) > 

CALL  FFT2D*  IS,  NC , NR,  M,  1 ) 

CAM  WBLKC ' IS, KFNAME ( 1 ) ) 

TVPE  "REFERENCE  INPUT- 

CALL  SDEF  i IS,  M,  IMN,  I AT , JREF  * 1 ) > 

CALL  FFT2D*  IS.  NC,  NR,  M.  1) 

CALL  WBLKC* IS, JFREF* D > 

C 

C OPEN  FILE  OF  COMPLEX  IMAGE  FREQUENCIES  TO  BE  READ 

C 

900  FORMAT* “ NAME  OF  FILE  TO  BE  READ  - ",Z> 

905  FORMAT  *S7 ) 

CALL  OPEN*  1,  KFNAME*  1 ).  2,  IER,  512) 

IF  HER  NE  DSTOP  --OPEN  ERROR-- 

C 

C CONJUGATE  REFERENCE  AND  PERFORM  POINT  TO  POINT  MULTIPLY 

C 

NBLOCK- 1 

DO  40  140-1,64 

IR-I40 

CALL  RU3l*'  l,  IR.  TEMF  . NBLOCK,  IER) 

IF  * IER  NE  DSTOP  - -RDBLK  ERROR- 
999  FOPMAT* IX, 8E14  4) 

DO  40  141-1,64 

IS*  14  I . I 10 ) “CONJG * IS*  141 , 140) >»TEMP* 141 ) 

40  CUN/IMUt 

C 

C CLOSE  IMAGE  FILE 

C 

CALI  Cl  USE  l 1,  IER) 

IF* IER  HE  DSTOP  --CLOSE  ERROR- • 

C 

C OPEN  REFFPENCC  FILE 

C 

CALL  OPEN* 1, JFREF* D, 2,  IER.  512) 

IF*  IER  NE  DSTOP  --OPEN  ERROR  — 

C 

C PERFORM  RECONSTRUCTION 

C 

DO  50  l it)*  l . 64 
IP-150 

CALL  RDBLK*  l,  IR,  TEMP,  NBLOCK,  IER) 

DO  30  1 5 1 » 1 . 64 

IS*  151  *50) -IS* 151,  1 50 ) »TEMP * 151) 

50  CON  i IN  XL 

C 

C PER) CRM  INVERSE  FFT 

C 

CALL  FFT2D*  IS.  NC- NR.  M.  2) 

C 

C STORE  PfC JNSTRUCTED  I MAGE 

C 

CALL  CLOSE* 1, IER) 

IF  t IER  Nfc  DSTOP  --CLOSE  ERROR- 
TYPE  "HOLOGRAM  OUTPUT" 

CALL  WBLKC I IS. KSUM* D ) 

STOP 

END 
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C HDRV  FR 

C 

C THIS  18  a DRIVER  PROGRAM  to  generate  MCH  OGR AM 

c data  files 

c 
c 

DIMENSION  KNAME  ( 7) 

COMPLEX  C VALUE 

COMPLEX  IS(6464)  TEMP (64) 

IMN-l 
M-6 
I AT- 1 

NO -NR 

ACCEPT  "VALUE  - ".VALUE 
CVALUE-CMPLX (VALUE.  0 ) 

CALL  SDEF( IS. M. IMN  IAT) 

CALL  WBLKC(IS) 

CALL  FFT 2D < IS. NC . NR  . M 1 > 

DO  22  1-164 
DO  22  J-l  64 
I S d • J ) -C  VALUE*  I S ( I.  J) 

22  CONTINUE 

CALL  FFT20( 16  NC. NR  H 2> 

CALL  WBLKC(IS) 

STOP 

CALL  SDEF * IS  M.  IMN-  IAT  > 

CALL  FFT2D( IS  NC  NR  M 1 ) 

CALL  WBLKC(IS) 

C 

C OPEN  FILE  OF  COMPLEX  IMAGE  FREQUENCIES  TO  BE  READ 

C 

WRI  TE<  10  <*00! 

900  F ORMAT < •'  NAME  OF  FILE  TO  BF  READ  - ".7> 

READ  < 11  90S) KNAME ^ l ) 

905  FORMAT (S7) 

CALL  OPEN< 1. KNAME< 1) . 2.  IER.  512) 

IF  * If R NE  t'STOP  --OPEN  ERROR 


C CONJUGATE  REFERENCE  AND  PERFORM  POINT  TO  POINT  MULTIPLY 

C 

NSLOCK-l 
DO  40  140-1. 64 
IB* I 40 

CALL  RDBLKd.IR  TEMP  . NBLOCK , IER  > 

IF  « IFR  NE  USTOP  RPHl  K f RROR 
999  FORMAT  < 1 X. 0E 1 4 4) 

DO  40  141-1  64 

lSd4l  140)-C0NJG(  IS(  141  140'  ' • T E HP  * 141  > 

40  CONTINUE 

C 

C CLOSE  IMAGE  FILE 

C 

CALL  CLOSE ( 1 IFR) 

IF  * IER  NE  DSTOP  --CLOSE  ERROR 
C 

C OPEN  REFERENCE  FILE 

C 

WRITE  « 10  900' 

HEAD (11  905'KNAHF \ \ ' 

CALL  OPEN d • KNAME d ' 2 IER  512) 

I F * I ER  NF  1)ST0P  OPFN  r RROR 

C 

C perform  reconstruction 

C 

DO  50  150-1  64 
IR-I50 

CALL  RDPlKd.IR  TEMP  NBLOCK.  IER  > 

DO  50  151-1  64 

IS* 151  I 30> -16(151  I 50 '• TEMP* 151 ) 

50  CONTI NUE 

C 

C PERFORM  INVERSE  FFT 

C 

CALI  FF T 2D < IS  NC  NR  M 2 ' 

C 

C STORE  RE  CONSTRUCTED  IMAGE 

C 

CAU  CLOSE <1  IER) 

IF*  IER  NE  l)6T0P  -Cl  08r  ERROR 
CALI  WBi  KC  < t S ' 

STOP 

END 


* 


4.2.3  HSDRV.FR 


HSDRV.FR  is  a FORTRAN  program  used  to  create  a multiplexed  hologram.  Tin- 
program  accepts  image  and  reference  integer  files,  computes  the  hologram,  and 
adds  the  result  to  a specified  hologram  file.  The  program  requires  a hologram 
file;  therefore,  CRSUM.SV  is  usually  run  prior  to  running  HSDKV.SV  to  create 
a file  in  the  proper  format  filled  with  zeros.  1 he  images  are  reconstructed 
using  the  MHOLO.SV  program. 

The  operator  and  program  interface  via  the  console  is  similar  in  terms  of 
file  inputs  to  the  FBDRVT.FR  program.  Refer  to  FBDRVT.FR  operator /program 
dialogue  (Paragraph  4.2.5)  for  file  formats. 

4.2.4  MHOLO.FR 

MHOLO.FR  is  a FORTRAN  program  which  reads  a specified  hologram  file  and 
performs  an  image  reconstruction  and  gray  scale  plot  for  a specified  reference 
file.  The  reference  file  must  contain  the  Fourier  Transform  of  the  reference 
image . 


The  operator  and  program  dialogue  uses  the  standard  machine  holograph  nomen- 
clature detail  in  the  FBDRVT.FR  description  (Paragraph  4.2.5). 
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C HSDRV  FR 

C 
C 
C 
C 
C 


THIS 

DATA 


IS  A DRIVER  PROGRAM  TO  GENERATE  HOLOORAN 
FILES 


DIMENSION  KNAME ( 7 ) . JREF ( 7 > , KSUM( 7 ) , KFNAME ( 7 > . JFREF < 7 ) 
COMPLEX  IS(64.  64).  TEMP ( 64 ) . CVALUE 
CVALUE-CMPLXd  ,0  > 

ACCEPT  " WEIGHTING  COEFFICIENT  ",  WC 
CVALUE-CMPLX(WC,  0 > 

WRITE! 10.  900) 

READ! 11. 9Q5)HSUM< 1 ) 

IMN*1 


M-6 


C 

C 


IAT*1 
NR-2**M 
NC  -NR 

ACCEPT  "NUMBER  OF  RUNS  “ , NRUNS 
DO  00  180-1.  NRUNS 
TYPE  "RUN  ",I0O 
TYPE  "SUBJECT  INPUT" 

CALL  SDEF ( IS.  M.  IMN.  IAT , KNAME ( 1 ) J 
CALL  FFT2D(  IS.  NC,  NR.  M,  1) 

CALL  WBLKC  < IS.  KFNAME  ( 1 > ) 

TYPE  "RUN  ". 180 
TYPE  "REFERNCE  INPUT  " 

CALL  SDEF(  IS.  M.  IMN,  IAT,  JREF ( 1 ) ) 

CALL  FFT2D< IS.  NC, NR, M,  1 ) 

CALL  WBLKC (IS,  JFREF ( 1)  > 

OPEN  FILE  OF  COMPLEX  IMAGE  FREQUENCIES  TO  BE  READ 


900  FORMAT ( " NAME  OF  FILE  TO  BE  READ  - ",Z> 
905  FORMAT ( S7 ) 

CALL  OPEN<  1 , KFNAME ( 1) . 2,  IER.  512 > 

I F ( I ER  NE  1IST0P  —OPEN  ERROR— 


(_  CONJUGATE  REFERENCE  AND  PERFORM  POINT  TO  POINT  MULTIPLY 

C 

NBLOCK-1 

DO  40  1 40-1 , 64 

IR-IaO 

CALL  RDBLKd.  I R.  TEMP.  NBLOCK.  IER) 

IF  ( TER  NE  DSTOP  — RDBLK  ERROR- 
999  f-URhATc  IX,  BE  14  4) 

DO  4A  141-1, 64 

ibl 141. 1 40 ) -CONJO ( IS ( 141, 140) >*TEMP< 141 ) 

40  CONTINUE 

C 

C CLOSE  IMAGE  FILE 

C 

CALL  CLOSE  < (.IER) 

IT (IER  NE  DSTOP  --CLOSE  ERROR  — 

C 

C OPEN  HOLOCRAM  SUM  FILE 

C 

C All  OPEN!  1,  K3UM<  1),  2,  IER.  512) 

IF  (IER  NE  DSTOP  —OPEN  ERROR  — 

C 

C PERFORM  SUMMATION 

C 

DO  50  150-1, 64 

io-r*>c 

CALL  RDBLKd.  IR,  TEMP,  NBLOCK.  IER) 

DO  50  151-1, 64 

IS (I  51,  1 50) -CVALUE* IS ( 151,  1 50) -TEMP (151 ) 

50  CONTINUE 

C 

C STORE  HOLOGRAM  SUM 

C 

CALL  CLOSE <1, IER ) 

IF  ( IER  NE  DSTOP  —CLOSE  ERROR- 
CALL  WBLKC ( IS,  KSOM( 1 ) ) 

80  CONTINUE 
bVQP 
END 
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MHOLO.FR  11/3/78 

THIS  PROGRAM  READS  A MULTIPLEXED  HOLOGRAM 
FROM  DISK. PERFORMS  IMAGE  RECONSTRUCTION 
AND  OUTPUTS  A GRAYSCALE  PLOT 


DIMENSION  KNAME ! 7 > , KSUM ( 7 ) 

COMPLEX  I S ! 64, 64 ) . TEMP ( 64 ) 

M-6 

NC-64 

NR  =6  4 


READ  FILE  CONTAINING  MULTIPLEXED  HOLOGRAM 


NBLOCK-1 

ACCEPT  "NUMBER  OF  RUNS  " , NRUNS 
DO  80  180*1 . NRUNS 
TYPE  "RUN  ". I 80 
TYPE  "HOLOGRAM  FILE" 

CALL  RBLKC(IS) 

INPUT  REFERENCE  IMAGE  TRANSFORM  (F-DOMAIN) 
WRITE! 10. 900) 

900  FORMAT < " NAME  OF  REFERENCE  IS  ",Z> 

READ ( 11, 905 ) KNAME  < 1 ) 

905  FORMAT (S7) 

CALL  OPEN! 1, KNAME ( 1 >, 2. IER, 512) 

IF  HER  NE.  1JST0P  — OPEN  ERROR — 

AND  DO  RECONSTRUCTION 

DO  50  150*1,64 
IR= 150 

CALL  HDBLK  < 1 . IR,  TEMP, NBLOCK,  IER) 

DO  50  151-1, 64 

IS ( 151,  1 50) -IS! 151,  I 50 ) *TEMP (151 ) 

50  CONTINUE 

PERFORM  INVERSE  TRANSFORM 

CALL  FFT2D! IS, NC. NR, M, 2) 

CALL  CLOSE! 1. IER) 

IF!  IER  NE  1)ST0P  —CLOSE  ERROR- 
CALL  GSPLT! IS, NC, NR, 3) 

WRITE! 12. 907 ) KNAME  < 1 ) 

907  FORMAT! 1 HO,  S7) 

80  CONTINUE 
STOP 
END 


4 - 2 (> 


4.2.5  FBDRVT 


l. 


t 


l 


FBDRVT.FR  is  a FORTRAN  driver  program  for  hologram  feedback  systems.  FBDRVT 
requires  the  following  relocatable  binary  files: 

WBLC.RB 

1READ.RB 

FFT2D.RB 

FNORM.RB 

CONS.RB 

RBLC.RB 

HCRE.RB 

HRCON.RB 

HADD.RB 

GSPLT.RB 

P1CPLT.RB 

PLPT.RB 

FORT. LB 


The  proper  load  statement  is  contained  in  file  FBDRVT. MC 
4.2.6  FBDRVT.FR 


FBDRVT.FR  is  a FORTRAN  driver  program  used  to  examine  multiplexed  holograms 
in  feedback  systems  and  is  illustrated  in  Figures  18  and  19.  The  program 
performs  the  equation 


H , + c R*S 

m = n-1 a 

n (1  + c RR*) 

where 

« hologram  in  Fourier  domain  at  iteration  n 

R - reference  image  in  Fourier  domain 

R*  - complex  conjugate  of  R 

S ■ Image  scene  in  Fourier  domain 
n 


(59) 


» ' 


■ 


j 


c - convergence  constant 

for  operator  specified  scene  and  reference  image  pairs. 
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FIGURE  19.  PSEUDO  FLOWCHART  FOR  FBDRVT 


Reconstruction  of  the  image  is  accomplished  by  using  the  formula 


I . - K-1  (R  -H)  (hf 

A A 

whore 

I,  « the  reconstructed  image 
A 

R.  » the  reference  in  the  Fourier  domain  paired  with  the  image  1. 

A A 

H - the  hologram  in  the  Fourier  domain 

and 


F~ ^ « the  inverse  Fourier  transform  operator 

The  driver  program  is  functionally  divided  into  two  divisions.  The  first 
division  sets  up  the  scene  and  reference  data  files  through  console  inter- 
face with  the  operator.  The  second  division  performs  the  hologram  feedback 
equation  (Equation  5d)  and  the  reconst ruct ion  equation  (Equation  bO)  for 
operator  specified  image  and  scene  pairs. 

In  the  first  division  the  program  asks  the  operator  for  files  containing 
integer  reference  images  and  integer  scene  images.  For  each  reference  and 
scene  pair  the  program  will  ask  for  four  file  names  for  files  the  program 
creates.  The  second  division  the  program  performs  the  feedback  equation  and 
image  reconstruction.  The  output  is  a grayscale  plot  for  each  operator 
spec  if ied  reconst ruct ion. 

The  dialogue  between  the  operator  and  the  program  is  as  follows: 

NUMBER  OF  PROGRAM  RUNS  - integer 

NUMBER  OF  POINTS  TO  BE  AVERAGED  IN  GSP1.T  - integer 

(NOTE:  The  grayscale  plot  finds  the  spec  it  led  number  ot  peak  (highest) 
values  and  averages  them  t o calculate  S/N  (see  GSPl.T. 
descr ipt ion) . ) 


CONVERGENCE  CONSTANT  « real number 

NAME  OF  HOLOGRAM  FILE  TO  BE  CREATED  - f i lename 

NUMBER  OF  REFERENCES  AND  SCENE  PAIRS  » integer 

(NOTE:  The  program  can  perform  the  feedback  equation  on  a maximum  of 

six  reference  and  scene  pairs,  and  the  hologram  reconstruction 
equation  on  a maximum  of  nine  references.  The  maximum  integer 
response  to  the  message  above  is  nine.) 

The  following  set  of  responses  is  repeated  for  the  specified  number  of 
reference  and  scene  pairs. 

REFERENCE  INPUT 

INPUT  FILENAME  = f i lename 

(NOTE:  The  specified  file  must  be  an  integer  image  in  the  standard 
format  (see  CRSCN  descr ipt ion) . ) 

NCI , NR  I b4,b4 

(NOTE:  The  number  of  columns  and  the  number  of  rows  of  all  hologram 

images  are  64  by  64.) 

INPUT  (X, Y)  COORDS  OF  UPPER  LH  CORNER  OF  IMAGE 

E.G.  (1,1)  = UPPER  LH  CORNER  OF  IMAGE  ON  DISK  - 1 , 1 

(NOTE:  All  hologram  images  have  coordinates  1,1  for  the  upper  left 
hand  corner  of  the  image.) 

IN  EFT 2D 

OUT  EFT 2D 

indicates  that  Fourier  transform  has  been  performed. 

NAME  OF  REF  FILE  IN  F DOMAIN  = filename 

(NOTE:  The  program  deletes  the  specified  file  and  creates  a new  file 

containing  the  b4  by  64  complex  array  of  the  reference  Fourier 
transform.  In  the  past  it  has  been  the  practice  to  start 
file  names  of  transformed  images  with  the  letter  F to  indicate 
the  data  is  in  the  Fourier  domain.) 

SCENE  INPUT 

INPUT  FILENAME  » f 1 lename 

(NOTE:  The  specified  file  must  be  an  integer  scene  in  the  standards 

format  (see  CRSCN  description).) 
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NCl.NRl  64,64 

INPUT  ( X , Y ) COORDS  OF  UPPER  LH  CORNER  OF  IMAGE 

E.C..  (1,1)  - UPPER  LH  CORNER  OF  IMAGE  ON  DISK  = 1,1 

IN  FFT2D 

OUT  EFT 2D 

NAME  OF  IMAGE  FILE  IN  F DOMAIN  - filename 

(NOTE:  The  program  deletes  the  specified  file  and  creates  a new  file 

containing  the  ti4  by  64  complex  array  of  the  scene  Fourier 
transform.  In  the  past  it  has  been  the  practice  to  start  file 
names  of  Fourier  transformed  files  with  the  letter  F to  indicate 
the  data  is  in  the  F domain.) 

NAME  OF  WORKING  FILE  = f i lename 

(NOTE:  The  program  creates  a file  containing  */(l  + c RR*)  (see 

Equation  59)  for  each  reference  and  scene  pair.  This  is  called 
the  working  file.  The  file  name  must  be  different  for  each 
reference  and  scene  pair.  If  a reference  is  to  be  used  for 
reconstruction  only,  a file  name  for  the  working  file  still 
must  be  specified  for  each  reference  and  scene  pair.) 

NAME  OF  CONSTANT  FILE  ' filename 


(NOTE:  The  program  creates  a file  containing 

c R*S 
n 

(1  + i RK* 

for  each  reference  and  scene  pair  (see  Equation  59).  This  is 
called  the  constant  file.  A different  file  name  must  be 
specified  for  each  reference  and  scene  pair.  It  a reference 
is  used  reconstruct  ion  onlv,  a file  name  for  the  constant  file 
still  must  be  specified  for  each  reference  and  scene  pair.) 

This  ends  the  file  setup  division  of  the  program. 

The  following  is  the  dialogue  between  the  operatin'  and  t tie  program  for  t lie 
feedback  and  image  reconstruct  ion  section  of  the  program. 

NUMBER  OF  ITERATIONS  = integer 

(NOTE:  The  program  will  perform  Equation  5^  the  specified  number 

of  times  and  then  stop,  or  return  to  the  beginning  of  t lie 
program. ) 

NUMBER  OF  ITERATIONS  FIR  PRINTOUT  - integer 


J 


h,  - 


(NOTE:  The  first  time  through  the  loop  the  program  computes 
Equation  59  only  once  prior  to  a grayscale  plot.) 

The  following  message  is  repeated  every  iteration  until  a printout. 

WHICH  REFERENCE  AND  SCENE  - integer 

(NOTE:  The  integer  response  indicates  which  reference  and  scene  pair 

is  used  in  the  hologram  feedback  equation  (Equation  59),  i.e.,  1 
means  the  first  pair,  2 the  second  pair,  etc.) 

HOW  MANY  PRINTOUTS?  integer 


The  program  will  perform  the  specified  number  of  reconstructions  (Equation  60). 
The  program  can  handle  up  to  nine  image  reconstructions,  provided  nine  refer- 
ence and  scene  pairs  were  input  in  the  first  section  of  the  program. 

The  following  is  repeated  for  the  specified  number  of  printouts: 

WHICH  REFERENCE  FOR  RECONSTRUCTION?  integer 

A 1 means  first  reference,  a 2 means  second  reference,  etc. 

IN  FFT2D 
OUT  FFT2D 

indicates  inverse  Fourier  transform  has  been  performed. 

ONE 

IN  GSPLT 

TWO 

THREE 

FOUR 

FIVE 

SIX 

SEVEN 

EIGHT 

After  the  specified  number  of  printouts  the  program  will  repeat  the  hologram 
iteration  routine.  After  the  specified  number  of  iterations  and  program  runs 
the  program  stops. 
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FBDRVT  FR  12/21/78 


THIS  ROUTINE  IS  A DRIVER  TO  CREATE  A HOLOGRAM 
FEEDBACK  LOOP 

INTEGER  SNAME ( 7 ) . RNAME < 7 ) , FRNAME 1(7), HNAME ( 7 ) . BNAME ( 7 ) , ANAME ( 7 ) 
INTEGER  FSNAME < 7 ) , WC0NS1 (7),  WCONS2C7),  WNAME1 (7) , WNAME2C7) 
INTEGER  FRNAME2 ( 7 ) , FRNAME3<  7 ) , FRNAME4 ( 7 > , FRNAME 5 ( 7 ) , FRNAMEe>  ( 7 ) 
INTEGER  WCGNS3 C 7 ) . WCONS4C7). WCONS5<7>,  WC0NS6 ( 7 > 

INTEGER  WNAME3 ( 7 > , WNAME4C7). WNAME5C7).  WNAME6C7) 

I NTEGER  FRNAME 7 ( 7 ) . FRNAME8 ( 7 ) , FRNAME9 < 7 ) 

COMPLEX  ISC64, 64) . CZERO.  CONE 


CONE-CMPLXC 1 ,0  ) 

M“6 
NC  -64 
NR-64 
IMN*  1 
I AT’  1 

C Z£RO=CHPLX ( 0 .0  ) 

ACCEPT  NUMBER  OF  PROGRAM  RUNS  » ", NPRUNS 
DO  303  1303*1 , NPRUNS 

ACCEPT  "NUMBER  OF  POINTS  TO  BE  AVERAGED  IN  CSPLT  - " , NP 
IFLAG*NP  + .3 


ACCEPT  "CONVERGENCE  CONSTANT  * '• , C 
NR  I TE ( 10, 900) 

900  FORMAT ( " NAME  OF  HOLOGRAM  FILE  TO  BE  CREATED  - ",Z> 
READ ill, 905 ) HNAME ( 1 ) 

905  FORMAT (S7) 

DO  10  110*1, 64 
DO  10  J10-1 , 64 
10  ISC  I 10, J10)*CZER0 

CALL  WBLC< IS, HNAME) 

READ  REFERENCE  IMACE  AND  STORE  TRANSFORM  ON  DISK 

ACCEPT  "NUMBER  OF  REFERENCES  AND  SCENE  PAIRS  - " ■ NREFS 
DO  30  130*1, NREFS 
TYPE  "REFERENCE  INPUT" 

CALL  SDEF( IS, M, IMN, IAT, RNAME ( 1 ) ) 

CALL  FFT2D( IS, NC. NR, M, 1 ) 

CALL  FNORM( IS, XNORM) 

XNORM* 1 /XNORM 
CALL  CONSC IS. XNORM) 

NR  I TEC  10, 999) 

999  FORMAT ( " NAME  OF  REF  FILE  IN  F DOMAIN  * " , Z ) 

IF ( 130  EG  1 ) READC 1 1 . 905 ) FRNAME 1(1) 

IFCI30  EG  2 ) READ  <11, 905 IFRNAME2 ( 1 ) 

IF ( 1 30  EG  3 ) READC 1 1 , 905 )FRNAME3( 1 > 

IF  C 130  EG  4 )READ  <11, 905 ) FRNAME4  < 1 > 

IF  C 130  EG  5 ) READ (11, 905 1FRNAME5 ( 1 ) 

IF ( 130  EG  6>READ< 1 1 , 905)FRNAME6( 1 ) 

IF (130  EG  7 ) READC 1 1 , 905 ) FRNAME 7 ( 1 ) 

IF  t 130  EG  8)READ( 1 1 , 905)FRNAMES( 1 ) 

IF i 130  EG  9 ) READ (11, 905 ) FRNAME9 ( 1 ) 

IF (130  EG  l ) C ALL  WBLC ( IS, FRNAME1 ) 

IF ( 130  EQ  21CALL  WBLC < IS, FRNAME2 ) 

I F < 1 30  EQ  3 ) CALL  WBLC < IS, FRNAME3 ) 

I F ( 130  EQ  4>CALL  WBLC l IS. FRNAME4 ) 

I F < 130  EQ  5JCALL  WBLC ( IS, FRNAME5 ) 

IFCI30  EQ  6 ) CALL  WBLC ( I S- FRNAME6 > 

IF (130  EG  7 (CALL  WBLC < IS, FRNAME7 ) 

IF ( 130  EG  8) CALL  WBLC < I S. FRNAME8 > 

IF <130  EQ  9 ) C ALL  WBLC < I S. FRNAME9 ) 

READ  IMAGE  AND  STORE  TRANSFORM  ON  DISK 
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TYPE  "SCENE  INPUT" 

CALL  SDEF ( IS. M. IMN, I AT , SNAME ( 1 ) > 

CALL  FFT2D! IS, NC. NR. M. 1 ) 

WRITE! 10, 998) 

998  FORMAT!"  NAME  OF  IMAGE  FILE  IN  F DOMAIN  - “,Z> 
READ ( 11 . 905  > FSNAME ! 1 ) 

CALL  WBLC! IS. FSNAME) 

READ  NAME  OF  HOLOGRAM  FILE  TO  BE  CREATED 


INPUT  NAMES  OF  WORKING  DISK  FILES 


1000 


1001 


WRITE! 
FORMAT 
IF( 130 
IF! 130 
IF! 130 
IF! 130 
IF! 130 
IF! 130 
WRITE! 
FORMAT 
IF! 130 
IF! 130 
IF! 130 
IF! 130 
IF ! 130 
I F < 130 


10, 1000) 

I “ NAME  OF  WORKING  FILE  = ",Z> 
) WNAME 1 ! 1 > 

) WNAME2 ( 1 ) 
)WNAME3( 1 > 

) WNAME4 ! 1 ) 

) WNAME 5 ! 1 ) 

) WNAME 6 ( 1 ) 


EG 

EQ 

EQ 

EO 

EQ 

EQ 


1 (READ! 11, 905) 
2) READ! 11. 905) 
3 (READ! 1 1 , 90S ) 
4 (READ! 1 1 , 905) 
5 ) READ! 11. 905) 
MREAD!  11, 905) 


10, 1001 ) 

! " NAME  OF  CONSTANT  FILE  » ",Z> 
(WC0NS1 ( 1 ) 

) WC0NS2  < 1 ) 

) WC0NS3 ( 1 ) 

) WC 0NS4 ! 1 > 

) WC0NS5 ! 1 > 

) WC0NS6 ! 1 ) 


EQ 

EQ 

EQ 

EQ 

EQ 

EQ 


1 )READ! 1 1 , 905) 

2 >READ! 11, 905) 

3 (READ! 11, 905) 
4) READ! 1 1. 905) 
5 >READ (11,  905) 
t> ) READ  (11, 905) 


PERFORM  R*R  CONJUGATE 


IF! 130 

EQ 

l (CALL 

IF( 130 

EQ 

2 (CALL 

IF l 130 

EO 

3 (CALL 

IF! 130 

EQ 

4 (CALL 

IF ( 130 

EQ 

5 (CALL 

IF  1130 

EQ 

6 (CALL 

IF! 130 

EQ 

1 (CALL 

IF i 130 

EQ 

2 ) CALL 

IF! 130 

EQ 

3 (CALL 

IF( 130 

EQ 

4 (CALL 

IF! 130 

EQ 

5 (CALL 

IF! 130 

EQ 

6 > CALL 

CALL  CONS! IS, C ) 


RBLC ! IS, FRNAME1 ) 
RBLC ( IS, FRNAME2 ) 
RBLC! IS. FRNAME3 ) 
RBLC ( IS, FRNAME4 ) 
RBLC!  IS, FRNAME5) 
RBLC! IS, FRNAME6 ) 
HCRE ( IS, FRNAME1 ) 
HCRE ( IS, FRNAME2 ) 
HCRE! IS, FRNAME3  > 
HCRE! IS. FRNAME4 ) 
HCRE! IS, FRNAME5 ) 
HCRE! IS, FRNAME6) 


FIND  1 / ! 1 + C * R«R  CONJ ) 


DO  20  120=1. NC 
DO  20  J2Q=1 , NR 

20  IS'  120, J20)=C0NE/!C0NE*IS( 120,  J20 ) ) 


IF! 130 

EQ 

IF! 130 

EO 

IF! 130 

EO 

IF! 130 

EQ 

IF! 130 

EO 

IF! 130 

EQ 

1 )  CALL  WBLC 

2)  CALL  WBLC 

3)  CALL  WBLC 
4 (CALL  WBLC 
5 ) C ALL  WBLC 
61CALL  WBLC 


! IS.  WC0NS1 ) 
( IS.  WC0NS2) 
( IS.  WC0NS3 ) 
! IS,  WCQNS4 ) 
( IS.  WC0NS5  > 
( IS. WC0NS6 ) 


FIND  C*RCONJ*S/(  ) 


C 


CALL  RBLC! IS. FSNAME) 


IF! 130 

EQ 

1 (CALL 

IF! 130 

EQ 

2 (CALL 

IF! 130 

EQ 

3 (CALL 

IF! 130 

EQ 

4 (CALL 

IF! 130 

EQ 

5 (CALL 

IF! 130 

EQ 

6 (CALL 

CALL  CONS! I 

S,  C ) 

IF! 130 

EQ 

1 (CALL 

IF! 130 

EQ 

2 (CALL 

IF! 130 

EQ 

3 (CALL 

IF! 130 

EQ 

4 (CALL 

IF! 130 

EQ 

5 (CALL 

HCRE! IS, FRNAME1 ) 
HCRE! IS, FRNAME2  > 
HCRE! IS, FRNAME3 ) 
HCRE! IS. FRNAME4 ) 
HCRE! IS, FRNAME5 ) 
HCRE! IS. FRNAME6) 


HR CON ( IS. WC0NS1 ) 
HR CON ( IS. WC0NS2 ) 
HRCON! IS. WC0NS3 ) 
HRCON ( IS. WC0NS4 ) 
HRCON! IS. WC0NS5 ) 
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A a a a ft  ftftft  ftftft  ftftft  ftftft  ftftft 


IF  * I 30 

EC 

61 CALL 

HRCOW*  IS  hCONSa  » 

IFiDO 

EC 

1 *CALL 

WOCC  * IS  AWE1I 

IF  « 130 

EC 

3*C  ALL 

MSEC* IS  HWAfttC 

IF* 130 

EC 

3 » CALL 

NSL.C*  IS  . HMAHt  j 1 

IF  * 130 

EC 

4 fCAUL 

AOILC*  IS  WHA-k  4 ; 

IF* 130 

EC 

3*  CAUL 

KLCUS-iAMEI' 

IF* 130 

EC 

A I CALL 

ICLCI  IS.  IAWEAI 

COMTIMlC 

ACCEPT 

~ «Mn*EP  OF 

HERAT  I CMS  - * - HR  .MS 

ACCEPT 

•'MtfW 

•ER  OF 

ITERATIONS  PEP  PPIMTOUT 

IPFEAO-O 

DO  300  I -JO- 1 St=  ."VS 

ACCEPT  "MHtlCM  REFERENCE  ANT  SCEWC  - * A* 


ej*C  «uvn 


CAUL  RILCdS  HWWEI 

NUL1IP1T  1%  IIUU  * c«as<i)*«2> 


IF*  ML 

EC 

1 ICALIL 

MRCQWlIS  hC  CMS  I 1 

I»  *!ML 

EC 

2* CALX 

HR  C 3W t IS  WCCIMS3 1 

IF*  ML 

EC 

3*CAUL 

HR  COW*  IS  WDCM53I 

IF  *ML 

EC 

4»CAUL 

HR  COW*  IS  HCCIHS4* 

IF  INI 

EC 

SLCAUL 

HR  COW*  IS  nCCINSM 

IF  (N» 

EC 

AICALX 

HR  CC1W*  IS  WCCMSaI 

Ad"  weighted  SCENE  hOlOORAW 

IF  (HI 

EC 

ItCALX 

HA DC*  IS  HHAHE  1 > 

1F«M» 

EC 

2ICALX 

HAOOilS  -ma-ex 

IFnlWL 

EC 

3 "CALL 

hAXC  . IS  HHAHF  3 * 

IF.  M» 

EC 

41  CALX 

IMOOtlS  HMAHE4I 

IF  *M% 

EC 

5ICALL 

HMQO*  IStMNNESt 

IFlMK 

EC 

A*  CALL 

HAD®*  IS  LEAK  & * 

HR  J IE 

M»M» 

TO  OISE 

CAUL  «®l  C«  IS-  *A*t  > 

IfilftlM  ME  0»«0  TO  Ml 
IP* LAC- 1 PF 

ftWtIP  IN®.  OBHAW  MECOMSIPUCT  lOM 


ACCEPT  •'MOW  IMMY  PPlPIfOtlfS  * * PPPIM 
OO  330  I JtJO  -I  NPPIM 

accept  *uhich  reference  top  reconstruction  y 


CALL  RRLCt 
IF* ML 4 EC 
IFuW**  FC 
IF* ML*  to 
IF«M*»  EC 
IF*M*4  EC 
ITAII  EC 
IFAlt  EC 
IF«M»4  EC 
IFAtl  EC 


IS.  HWAHE 
I "CALL 
2>CMX 
JICALL 
* I CALL 
VCAL 

6 » CALL 

7 "CALL 
9 "CALL 
<»»CAUL 


HR  CCAM 
HRCONU 
HRCCM* 
HRCOWL 
HRCOWL 
HRCONL 
HR  COWL 
HRC'OW  l 
HR  COW  I 


IS  F«MAMES» 
is.  FRMAAE.2* 
IS- FHMAME3* 
IS.  FRNAAL4* 
IS  r«WEJI 
IS-  frnauca* 
is  f«wwri 
is  rnwncai 
IS-  FRMAHER* 


CALL  FT!M»IS.WC  WB  « Jl 
CAt  II  C^fttTdlS  MC  HR  IFtLAC* 

in  COM!  1MUE 


LOOP  TO  DIFFERENCE  FOR  NUT  ITERATION 

.VI  JP*L  AC-IPFLAC-  1 
AM  CONI  IML4E 
M I COAT  IML« 

CALL  RESET 

STOP 

I* 


IPF 


*.  Kit 


4-  Tft 


nnnnnnn 


I 

I 

I 


r 


CRSUM  FR  11/3/78 

THIS  PROGRAMS  WRITES  A BLOCK  OF  COMPLEX  ZEROS  ON 
A DISK  FILE 

SUBROUTINES  USED:  WBLKC.RB 

COMPLEX  IS(64. 64). CZERO 
CZERO-CMPLX<0.  . O ) 

DO  10  1-1.64  r 

DO  10  J-l. 64 
IS( I.  J) -CZERO 
10  CONTINUE 

CALL  WBLKCtlS. NAME) 

STOP 
END 

I 


I 

I 
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A. 2. 7 RBLC.FR 


RBLC.FR  is  a FORTRAN  subroutine  that  reads  a 64  by  64  complex  array  from  disk 
as  illustrated  In  Figure  20.  The  routine  uses  channel  1 for  the  disk  transfer. 

Two  arguments  are  passed  in  the  call  statement.  The  first  is  the  name  of 
the  complex  array  in  memory  which  will  receive  the  disk  data.  The  second  is 
an  integer  array  that  contains  the  file  name  of  the  data  on  disk.  The  CAM. 
statement  is  of  the  form 

CALL  RBLC  (arrayname,  filename  array) 


The  subroutine  assumes  the  data  on  disk  is  a 64  by  64  complex  array  with  the 
first  row  of  data  in  BLOCK  1 of  the  file. 

4.2.8  RDBLC.FR 

RDBLC.FR  is  a FORTRAN  subroutine  that  reads  a file  name  from  the  console,  and 
then  reads  a 64  by  64  complex  array  from  the  specified  file  (see  Figure  21). 

The  subroutine  outputs  to  the  console  the  message: 

"NAME  OF  FILE  TO  BE  READ  = " 

to  which  the  operator  responds  with  a file  name  (seven  characters  maximum). 

A 64  by  64  complex  array  is  then  read  from  the  specified  disk  tile.  The 
subroutine  assumes  the  data  on  disk  is  a 64  by  64  complex  array  with  the 
first  row  of  data  in  BLOCK  1 of  the  file. 

The  complex  array  is  returned  to  the  calling  program.  The  CALL  statement  is 
of  t he  f o rm 

CALL  RBLKC  (arrayname) 

4.2.9  WBLC.FR 

WBLC.FR  is  the  name  of  a FORTRAN  subroutine  that  is  used  to  write  a 64  by  64 
complex  array  out  to  disk  (see  Figure  22).  The  subroutine  uses  channel  1 for 
the  disk  transfer. 

Two  arguments  are  passed  to  the  subroutine.  The  first  is  the  array  name  of 
the  complex  data  which  is  to  be  written  out  to  disk.  The  second  is  an  integer 
array  which  contains  the  file  name  of  the  disk  file  to  be  created.  The  CALL 
statement  is  of  the  form: 

CALL  WBLC  (arrayname,  filename  array) 


4-  18 


ci  n ci 


64  COMPLEX  VALUES 
266  WORDS 


BLOCK  1 


BLOCK  64 


• KNAME  | ) IS  READ  FROM  CONSOLE 

• FILE  KNAME  ( ) IS  READ  FROM 

DISK  INTO  ARRAY  IS  ( . I 

• IS  I . I IS  PASSED  TO  CALLING 
PROGRAM 


C RBLKC  FR  10/23/78 

C 

SUBROUTINE  RBLKC lIS) 

THIS  SUBROUTINE  READS  A 64  X 64  COMPLEX  ARRAY 

COMPLEX  IS(64, 64) 

DIMENSION  KNAME (7) 

C 

WRITE (10. 900) 

900  FORMAT ( " NAME  OF  FILE  TO  BE  READ  - ".Z) 

READ (11. 905  > KNAME ( 1 ) 

905  FORMAT ( S7 ) 

CALL  OPEN ( 1 , KNAME ( 1 ) . 2.  IER.512) 

IFUER  NE  1)  STOP  —OPEN  ERROR  — 

NBLQCK“64 
IR  = 1 

CALL  RDBLK ( 1 . IR,  IS. NBLOCK,  IER ) 

IFUER  NE  1 > STOP  — RRBLK  ERR— 

CALL  CLOSE (1.  IER) 

IFUER  NE  1 >STOP  —CLOSE  ERR  — 

RETURN 

END 


FIGURE  20.  RBLKC.FR 
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64  COMPLEX  VALUES 
266  WORDS 


BLOCK  J 


BLOCK  64 


• KNAME  ( i IS  PASSED  TO 
SUBROUTINE 

• IS  I I IS  READ  E ROM  DISK  PILE  K NAME  ( ) 
AND  PASSED  TO  CALLING  PROGRAM 

• KNAME  IS  RETURNEO  TO  CALLING 
PROGRAM  UNCHANGED 


RBLO  Fk  12/6/79 

SUDROUT iNE  RBLC ( 15, KNAME > 

C THIS  SUBROUTINE  READS  A 64  X 64  COMPLEX  ARRAY 

C 

COMPLEX  I S ' 64 , 64  ) 

DIMENSION  KNAME ( 7 ) 

C 


CALL  OPENM  KNAME!  1 > . 2,  IER,  512) 
IFiIER  NE  i)  STOP  — OPEN  ERROR  — 
NBlOCK=64 
IR-l 

CALL  RDBLK ( 1 , IR, IS, NBLQCH, IER) 

IF  HER  NE  1 > STOP  --RRBLK  ERR-- 
CALL  CLOSE! 1 , IER ) 

IF  ( IER  NE  DSTOP  --CLOSE  ERR-- 

RETURN 

END 
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64  COMPLEX  VALUES 
266  WORDS 
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• IS  164  641  PASSED  TO  SUBROUTINE 

• KNAME  I7|  PASSED  TO  CALLING  PROGRAM 
IS  <64  641  RETURNED  TO  CALLING 
PROGRAM  UNCHANGED 


C WBLKC  PR  10/23/70 

C 

SUBROUTINE.  WBlKC  i I S KNAME  ' 

C 

C THIS  SUBROUTINE  WRITES  A £,4  X 64  COUPLET  ARRAY 

C 

COMPLEX  1 S ! 64  64  1 
DIMENSION  I DD ( 256  > . KNAME I 7 ) 

C 

WRITE! 10. 900) 

900  FORMAT!"  NAME  OF  FILE  TO  BE  CREATED  - ",  Z > 

READ! 1 1 903>  KNAME 1 1 ) 

903  FORMAT <S7) 

CALL  DF ILW ( KNAMF  < 1 > , IER  > 

CALL  CF  ILW  ( KNAME  ( 1 I 2,  IER) 

IF! IER  NE  I)  STOP  --  CFILW  ERROR  IBCD)  -- 
CALL  OPEN! 1. KNAME! 1 ) 2 IER. 312) 

IF! IER  NE  l)  STOP  --OPEN  ERROR  (BOD)  -- 
IDD(2) *64 
I DD( 3) • 1 
IDDi 4 ) -64 
I DD i 5 > » 1 
C 

C WRITE  HEADER  INFORMATION 

C 

IR-0 

CAlL  WR3LK!!.  IR.  IDD  |.  IER) 

IFiIER  NE  DSTOP  --WRBLH  ERR-- 

NBl0CK'64 

IR*  I 

call  WRBLKIl.  IR.  I S.  NBLOCR  . IER) 

IFltER  NE  DSTOP  --WRBLK  ERR- - 

call  close  < i . ier  > 

IF ! IER  NE  DSTOP  --CLOSE  ERR- 
REToRN 

END  * 


FIGURE  22.  WBLKC.FR 
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The  subroutine  writes  some  header  information  (number  of  rows  and  columns)  into 
BLOCK  0 of  the  file,  then  writes  the  array  out  to  disk  with  the  first  row  of 
the  array  in  BLOCK  1 of  the  file.  Subsequent  rows  are  written  into  sub- 
sequent blocks. 

4.2.10  WBLKC.FR 

WBLKC.FR  is  a FORTRAN  subroutine  that  accepts  a file  name  from  the  console 
and  writes  a b4  by  b4  complex  array  out  to  the  specified  disk  file  (see 
Figure  23).  The  subroutine  outputs  to  the  console  the  message: 

"NAME  OF  FILE  TO  BE  CREATED  - " 

to  which  the  operator  responds  with  a file  name  and  carriage  return  (seven 
characters  maximum) 

The  subroutine  writes  some  header  information  (number  of  rows  and  columns)  into 
BLOCK  0 of  the  array,  then  writes  the  64  by  64  complex  array  (passed  from  the 
calling  program)  out  to  disk  with  the  first  row  of  data  in  BLOCK  1 of  the 
specified  file.  Subsequent  rows  are  written  into  subsequent  blocks. 

The  subroutine  returns  an  integer  array  to  the  calling  program  containing 
the  specified  disk  file  name.  The  CALL  statement  is  of  the  form: 

CALL  WBLKC  (arrayname,  filename  array) 
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r 
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• AUGMENTS  IS  AND  KNAME  PASSED 
TO  SUBROUTINE 


• IS  AND  KNAME  RETURNED  TO 
CALLING  PROGRAM  UNCHANGED 


C WBLC  FR  12/6/70 

C 

SUBROUTINE  WBLC ( IS. KNAME) 

C 

C THIS  SUBROUTINE  WRITES  A 64  X 64  COMPLEX  ARRAY 
C 

COMPLEX  I S < 64. 64 ) 

DIMENSION  I DD ( 256 ) . KNAME ( 7 ) IJ 

C 


C 

C 

C 


CALL  DF I LW ( KNAME ( 1 ) . I ER ) 

CALL  CF ILW( KNAME ( 1 > , 2.  IER) 

IF ( IER  NE  1)  STOP  — CFILW  ERROR  (BCD)  — 
CALL  OPENi 1.  KNAME ( 1) , 2,  IER.  512) 

IF(  IER  NE  1)  STOP  --OPEN  ERROR  (BCD)  — 

IDD<2>«64 

IDD(3)-1 

I DD ( 4 ) “64 

IDDi 5)-l 

WRITE  HEADER  INFORMATION 
1R-0 

CALL  WRBLK ( 1 , IR.  I DD.  1,  IER) 

IF < IER  NE  1 (STOP  — WRBLK  ERR  — 

NBLOCK-64 

IR-1 

CALL  WRBLM1,  IR,  IS,  NBLOCK,  IER) 

IF  (IER  NE  1 ) STOP  —WRBLK  ERR  — 

CALL  CLOSE ( 1 . IER ) 

IF (IER  NE  1 (STOP  —CLOSE  ERR- 

RETURN 

END 

FIGURE  23.  WBLC .FR 
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4 . 3 UT 1 L 1 TY  SUBROUT 1 NES 


There  are  six  utility  subroutines,  each  of  which  performs  one  simple  operation 
on  a complex  64  by  64  array:  (1)  DIFF.-,  (2)  HADD.-,  (3)  HCRE , (41  HRCON.-, 
(5)  FNORM.— , and  (6)  CONS.—,  Each  of  these  routines  use  channel  1 if  data 
needs  to  be  transferred  to  or  from  disk.  The  subroutines  open  and  close 
channel  1,  and  will  cause  a program  ABORT  if  any  1/0  errors  occur  on  channel  1. 

DIFF.—  computes  the  difference  between  two  complex  arrays.  HADD.—  adds  two 
complex  arrays.  CONS.—  multiplies  a complex  array  by  a real  number.  FNORM.— 
computes  the  square  root  of  the  sum  of  the  squares  of  a complex  array. 

HRCON.—  multiplies  two  complex  arrays  on  an  element  by  element  basis.  HCRE.— 
multiplies  a complex  array  by  the  complex  conjugate  of  a different  complex 
array  on  an  element  by  element  basis. 

4.3.1  DIFF. FR 

DIFF.FR  is  a FORTRAN  subroutine  named  DIFF.  The  subroutine  takes  the 
difference  of  two  complex  64  by  64  arrays  element  by  element,  one  array  being 
in  memory  and  the  other  array  on  disk. 

Two  parameters  are  passed  to  the  subroutine.  The  first  is  the  b4  by  b4 
complex  array  (in  memory).  The  second  is  an  integer  array  containing  the  file 
name  of  the  complex  array  residing  on  disk.  The  call  statement  is  of  the  form 

CALL  DIFF  (arravname,  filename  array) 

The  subroutine  assumes  the  disk  data  is  a 64  by  64  complex  array  with  the 
1st  row  of  data  in  BLOCK  1 of  the  disk  file. 

The  data  in  memory  is  subtracted  from  the  data  on  disk,  witli  the  result 
stored  in  the  64  by  64  complex  array  and  returned  to  the  calling  program. 

4.3.2  HADD . FR 

HADD.FR  is  a FORTRAN  subroutine  named  HADD.  The  subroutine  adds,  element  by 
element,  a 64  by  64  complex  array  in  memory,  with  a 64  by  64  complex  array  on 
disk.  Channel  l is  used  by  the  subroutine. 

Two  parameters  are  passed  to  the  subroutine.  The  first  is  the  64  by  64 
complex  array  in  memory.  The  second  is  an  integer  array  containing  the  tile 
name  of  the  complex  array  residing  on  disk.  The  call  statement  is  ot  the 
form: 

CALL  HADD  (arrayname,  filename  array) 

The  subroutine  assumes  the  disk  data  is  a 64  by  64  complex  array  with  the 
first  row  of  data  in  BLOCK  l of  the  disk  file.  Subsequent  rows  are  in  sub- 
sequent blocks. 

The  two  arrays  are  added,  with  the  resulting  64  by  64  complex  array  returned 
to  the  calling  program. 
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4.3.3  FNORM.FR 

FNORM.FR  is  a FORTRAN  subroutine  named  FNORM.  The  subroutine  calculates  the 
square  root  of  the  sum  of  the  squares  of  the  elements  in  a 64  by  64  complex 
array.  The  result  can  be  used  to  normalize  the  complex  array. 

Two  arguments  are  passed  in  the  CALL  statement.  The  first  is  the  64  by  64 
complex  array.  The  second  argument  is  the  square  root  of  the  sum  of  the 
squares  returned  by  the  subroutine  to  the  calling  program.  The  CALL  state- 
ment is  of  the  form: 

CALL  FNORM  (arrayname,  realvariable) 

The  64  by  64  complex  array  is  returned  to  the  calling  program  unchanged. 

4.3.4  CONS. FR 

CONS.FR  is  a FORTRAN  subroutine  named  CONS.  The  subroutine  multiplies  each 
element  of  a 64  by  64  complex  array  by  a real  variable. 

Two  arguments  are  passed  in  the  CALL  statement.  The  first  is  the  64  by  b4 
complex  array.  The  second  is  the  real  variable.  The  CALL  statement  is  of  the 
form 


CALL  CONS  (arrayname,  realvariable) 

The  resulting  64  by  64  complex  array  is  returned  to  the  calling  program. 

4.3.5  HRCON.FR 

HRCON.FR  is  a FORTRAN  subroutine  named  HRCON.  The  subroutine  multiplies  a 
64  by  64  complex  array  in  memory  with  one  on  disk,  on  an  element  by  element 
basis. 

Two  arguments  are  passed  from  the  calling  program  to  the  subroutine.  The 
first  is  the  64  by  64  complex  array  in  memory.  The  second  is  an  integer  array 
containing  the  file  name  of  the  disk  file.  The  CALL  statement  is  of  the  form 

CALL  HRCON  (arrayname,  filename  array) 

The  subroutine  assumes  the  disk  file  is  a 64  by  b4  complex  array  with  the  first 
row  of  data  in  BLOCK  1 of  the  disk  file,  with  sequential  rows  in  sequential 
file  blocks. 

The  resulting  complex  array  is  returned  to  the  calling  program. 
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A. 3.6  HCRE.FR 


HCRE.FR  is  a FORTRAN  subroutine  named  NCRE.  The  subroutine  multiplies  a 
6A  by  6A  complex  array  in  memory  with  the  complex  conjugate  of  a 6A  by  6A 
complex  array  on  disk,  on  an  element  by  element  basis. 

Two  arguments  are  passed  from  the  calling  program  to  the  subroutine.  The 
first  is  the  6A  by  64  complex  array  in  memory.  The  second  is  an  integer  array 
containing  the  file  name  of  the  disk  file.  The  CALL  statement  is  of  the  form: 

CALL  HCRE  (arrayname,  filename  array) 

The  subroutine  assumes  the  disk  file  is  a 64  by  64  complex  array  with  the  first 
row  of  data  in  BLOCK  1 of  the  file. 


The  subroutine  returns  the  resulting  64  by  64  complex  array  to  the  calling 
program. 


nfinn  r> 


C DIFF.  FR  12/5/78 


THIS  SUBROUTINE  TAKES  THE  DIFFERNCE  OF  TWO  ARRAYS 
ONE  PASSED  IN  THE  ARGUMENT  LIST  AND  THE  OTHER  ON 
DISK 


SUBROUTINE  DIFF( IS. CSNAME ) 
COMPLEX  IS ( 64. 64 ) . CTEMP ( 64 > 
INTEGER  CSNAME ( 7 ) 


CALL  OPEN ( 1. CSNAME. 2. IER. 512) 

IF<  IER  .NE.  1 ) STOP  — S<X>  OPEN  ERROR— 

NBLOCK-1 

DO  10  110-1.64 
IR-I10 

CALL  RDBLK(1,  IR.  CTEMP. NBLOCK,  IER) 

IF(  IER  . NE  1 ) STOP  — RDBLK  ERROR- 
DO  20  J20-1.64 

IS( J20.  1 10 ) -CTEMP ( J20 ) - IS ( J20.  110) 

20  CONTINUE 
10  CONTINUE 

CALL  CLOSE ( 1 > IER) 

IF(  IER  NE  DSTOP  — CLOSE  ERROR — 

RETURN 

END 


I 


n n n n n 


HADD.  FR  12/5/78 

THIS  SUBROUTINE  ADDS  TWO  ARRAYS,  ONE  IN  MEMORY 
WITH  ONE  ON  DISK 

SUBROUTINE  HADD( IS, CSNAME ) 

COMPLEX  IS (64, 64 > , CTEMP ( 64) 

INTEGER  CSNAME ( 7 ) 

C 

CALL  OPENd. CSNAME. 2,  IER, 512) 

IF ( IER  .NE.  l)STOP  ~S(X>  OPEN  ERROR— 

C 

NBLOCK-1 

C 

DO  10  1 10“1 1 64 
IR-I10 

CALL  RDBLK(1,  IR,  CTEMP,  NBLOCK,  IER) 

IF ( IER  . NE.  DSTOP  — RDBLK  ERROR — 

DO  20  J20-1.64 

IS ( J20,  I10)-IS(J20,  1 1 0 ) +C TEMP ( J20 ) 

20  CONTINUE 
10  CONTINUE 

CALL  CLOSE (1,  IER) 

IF(  IER  .NE.  1)  STOP  —CLOSE  ERROR- 

RETURN 

END 


4-48 


C FNORM.  FR  12/7/78 

C 

C THIS  SUBROUTINE  FINDS  THE  SORT  OF  THE 

C SUM  OF  THE  SQUARES  OF  ELEMENTS  IN 

C A COMPLEX  ARRAY 

C 

SUBROUTINE  FNORM ( IS, XNORM > 

COMPLEX  IS(64i 64) , CSUM 
C 

CSUM«CMPLX(0  . 0.  ) 

DO  10  110*1. 64 
DO  10  J10-1. 64 

csum«csum+is<  no,  jio)*is<  no,  uio> 

10  CONTINUE 

XNORM-SQRT< CABS < CSUM > ) 

RETURN 

END 


4-49 


AD-A068  638 


UNCLASSIFIED 


FORD  AEROSPACE  AND  COMMUNICATIONS  CORP  NEWPORT  BEACH  —ETC  F/6 
MACHINE  HOLOGRAPHY. (U) 

MAR  79  C L RICHARDS  N00019-78-C-0659 

U-6516  NL 


F/6  19/5 


n 


non  o 


CONS.  FR  12/5/78 

THIS  SUBROUTINE  MULTIPLIES  A COMPLEX  ARRAY  BY  A CONSTANT 

SUBROUTINE  CONS<IS.  C> 

COMPLEX  IS <64.  64) » COMC 

COMC=CMPLX<C,  0.  > 

DO  10  110=1.64 
DO  20  J20=l.  64 
IS ( J20. I 10)=IS( J20. I 10)*C0MC 
20  CONTINUE 
10  CONTINUE 
RETURN 
END 


4-50 


O U (J  o 


HRCO N.  F R 12/5/78 

THIS  SUBROUTINE  MULTIPLIES  AN  ARRAY  IN  MEMORY 
BY  ONE  ON  DISK  .ELEMENT  BY  ELEMENT. 

SUBROUTINE  HRCON( IS,  CSNAME)  ' 

COMPLEX  IS (64.  64) . CTEMP (64) 

INTEGER  CSNAME(7) 

C 

CALL  OPEN(l,  CSNAME, 2,  IER,  512) 

IF ( IER  . NE.  DSTOP  --S(X)  OPEN  ERROR— 

C 

NBLOCK-1 

C 

DO  10  110*1.64 
IR-I10 

CALL  RDBLK(1.  IR, CTEMP,  NBLOCK,  IER) 

IF ( IER  . NE.  DSTOP  — RDBLK  ERROR- 
DO  20  J20-1.64 

IS( J20, I 10 ) *IS( J20, 110) *CTEMP ( J20 ) 

20  CONTINUE 
10  CONTINUE 

CALL  CLOSE d,  IER) 

IF ( IER  . NE  DSTOP  — CLOSE  ERROR — 

RETURN 

END 


4-51 


~ i “ 
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HCRE  FR  12/5/78 

THIS  SUBROUTINE  MULTIPLIES  AN  ARRAY  IN  MEMORY 
BY  THE  CONJUGATE  OF  AN  ARRAY  ON  DISK  TO 
PRODUCE  A HOLOGRAM 

SUBROUTINE  HCRE< IS. CSNAME ) 

COMPLEX  IS(64.  64). CTEMP(64) 

INTEGER  CSNAME ( 7 ) 

C 

CALL  OPEN  < 1.  CSNAME.  2. IER. 512) 

IF < IER  NE  1 ) STOP  — S(X)  OPEN  ERROR — 

C 

NBLQCK-1 

C 

DO  10  1 10-1 . 64 
IP-110 

CALL  RDBLK < 1 . IR, CTEMP, NBLOCK.  IER) 

IF  ( IER  NE  llSTOP  —RDBLK  ERROR- 
DO  20  J20-1.64 

IS  < J20.  1 10) -IS ( J20.  110) *CONJO (CTEMP ( J20 ) > 

20  CONTINUE 
10  CONTINUE 

CALL  CLOSE ( 1. IER) 

IF  ( IER  NE  1 ) STOP  —CLOSE  ERROR- 

RETURN 

END 


CFFT2D  FR  04/11/78 

SUBROUTINE  FFT2D< A, NC. NR. M. INV) 
COMPLEX  A<64,  64).  W<6,  2).  T.  U 
TYPE  "IN  FFT2D" 

N-2**M 

FN-< 1 /N)»*2 

NV2-N/2 

NM1-N-1 

DO  90  L-l.  M 

AN0--2  »3  141596/2*»L 

W<L.  1 ) “CMHLX ( COS < ANC ) . SIN*  ANG  > > 

50  W<L. 2»-CMPLX<C0S<-ANC>.  SIN(-ANC) ) 

DO  26  J2- 1 . 2 
J-l 

DO  7 I « 1 , NM1 
IF  < I GE  J)  CO  TO  5 
DO  4 KK- 1 . N 
T-A<  J.  KK) 

A< J. KK)-A< I. KK) 

4 A < 1 . KK  > ■ T 

5 K-NV2 

6 IF (K  GE  J)  00  TO  7 

J-J-K 

K-K/2 
GO  TO  6 

7 J«J*K 
LE-1 

DO  20  l «■! . M 
LE1-L.E 
LE-LEHE 
U-  < 1 .0  > 

DO  20  J> 1 . LE 1 
DO  10  1 » J.  N. LE 
IP-IH  E 1 
DO  10  KK  1 . N 
T-A< IP. KK)»U 
A ( IP.  KK)-A< I. KK ) -T 
10  A<I,KK)*A(I.KK)*T 

20  U-U*WU,INV> 

DO  25  1*1.  N 
DO  25  KKr- 1 . N 
T-A< 1. KK) 

A< I . KK)-A(KK.  1 ) 

25  A<KK.  I >*  T 

26  CONTINUE 

IF< INV  EO  1 ) CO  TO  55 

IF ( IFT  FO  2)  FN«FN/<N**?> 

DO  30  1 * 1 . N 
DO  30  KK*  1 . N 

30  A< I . KK ) «A< I • KK ) *FN 

55  TYPE  "OUT  f FT2D’* 

RETURN 

END 


4-S ) 


5/10/78 


REV  A 11/1/78 
REV  B 11/2/78 
REV  C 11/28/78 


C GSPL  T FR 
C 
c 
c 

SUBROUTINE  CSPLT  < IS.  NC . NR.  Of-  LG  > 

C 

C CSPLT  MAKES  A OREY  SCALE  PLOT  OF  THE  MAGNITUDE 

C OF  THE  NR  X NC  COMPLEX  ARRAY  IS 

C GFLC-  0 NEW  LINEAR  ORAY  SCALE 

C 1 PREVIOUS  CRAY  SCALE 

C 2 LOG  ORAY  SCALE 

C 3 NORMALIZE  PEAK  VALUE  TO  99 

C 4 NORMALIZE  PEAK  VALUE  TO  1000 

CQNMGN/CG/VAL  < 10  > . CRAT(30> 

COMMON/ JJ/LBUFI 132.  3) 

DIMENSION  KBUF ( 260 ) . TAB ( 1024) . AVE <64 ) 

integer  tab. val. orat. gflc 

COMPLEX  ISC  64, 64) , DD<64>. CPKMO. CPNORM.  CSUM 

C 

DATA  LBUF) 131,  1 > /S.  12K/ 

DATA  LBUF( 131. 2)/3.  12K/ 

DATA  LBUF( 131. 3)/3.  12M/ 

TYPE  "ONE" 

C 

C 

C INPUT  NUMBER  OF  POINTS  TO  BE  AVERAGED 

C ZERO  AVE  ARRAY 

C 

ACCEPT  "NUMBER  OF  POINTS  AVERAGED" , NAVE 
DO  202  1202-1.64 
AVE  < 1 202 ) *0 
202  CONTINUE 
C 

c 

TYPE" IN  GSPLT" 

C 

DATA  0RAT/100K.  100K,  122K.  11  IK.  122K.  155R. 4*177K. 

* 100K.  122K.  100K,  144K,  155K,  3*122K.  155K,  177K. 

* 100K,  100K.  122K.  1 1 IK,  122K.  155K, 4*177K/ 

C 

CALL  OPEN l 2. "SCRATCH" , 0. IER, 512) 

C 

DO  1 11-1,260 
1 KBUF (II) -0 

C 

DO  23  I 1-1. NR 
DO  21  12-1 . NC 
21  DD( 1 2 > — I S ( 12.  11 ) 

CALL  WRBLM2.  X 1-1 . DD.  1,  IER) 

IF<  IER  NE  1)  STOP  — URBLK  ERR  — 

23  CONTINUE 

30  CONTINUE 

TYPE  "TWO" 

PKMG-0 
DO  5 J- 1 , NR 
DO  5 I-l.NC 
SS-CABS) IS( I. U) ) 

IF ( GFLO  EO  2)  SS-10  *AL0G10(SS*1  ) 

IF ( SS  OT  PKMG ) PKMO-SS 
5 CONTINUE 

C 

TYPE  " THREE" 

C NORMALIZE  INPUT  ARRAY 

C 

IFcjFLQ  LT  3)00  TO  07 


on  noon  n o ooooooo  ooo 


00  68  188-1.  NC 
DO  88  U88-».NR 

!S< 188. J88)-IS< 188.  JB8)*CPKHC 

FIND  THE  "NAVE"  NUMBER  OF  PEAK  POINTS 

SS-CAB6< I S < 188. J88> ) 

IF ( SS  LE  AVE ( NAVE ) ) 00  TO  204 
AVE ( NAVE ) -SS 
NAVEL-NAVE- 1 
00  206  1206-1. NAVEL 

IF ( AVE < NAVE* 1-1206)  LE  AVE < NAVE- 1206 > >00  TO  204 
TEMP -AVE  t NAVE- 1 206  > 

A VE  < NAVE - 1 206  > -AVE  < NAVE ♦ 1 - I 206 ) 

AVE < NAVE* 1 - 1 206  > -TEMP 
206  CONTINUE 
204  CONTINUE 


88  CONTINUE 

TYPE  "FOUR" 


CALCULATE  STATISTICS  OF  BACKGROUND 

87  CONTINUE 

CSUM-CMPUUO  ,0  ) 

BSUH-O 

8SUM2-0 

DO  89  1 89-1 , NC 
A-CABS l IS i 189. 1 > > 

B-CABS  < IS< 189. NR  > ) 

CSUM-CSUM*IS< 189,  l > ♦ I S « 189. NR) 

I S < 189,  1 * -CPNORM 
ISt 189. NR  > -CPNORM 
BSUM2«BSUM2*A*A*B*B 
89  CONTINUE 

TYPE  "FIVE" 

NROWS-NR  1 

DO  91  191-2. NROWS 
A-CABS1  ISU  191  ) ) 

B-CABSl ISiNR. 191 ) ) 

C SUM-C  SUM+ 1 S < 1,  191 > ♦ IS  < NC . 191 ) 

ISt 1. 191 ) -CPNORM 
IStNC, 191 > -CPNORM 
BSUM2-BSUM2*A*A*B*B 
91  CONTINUE 

TYPE  "SIX" 


CALCULATE  AVERAGE  OF  THE  PEAK  VALUES 
PAVE-0 

DO  208  1208-1.  NAVE 
PAVE-PAVE+AVE l 1208) 

208  CONTINUE 

PAVE-PAVE/FLOAT ( NAVE ) 


BSUM-CABS(CSUM) 

BMEAN-BSUM/ ( 2*NC*2*NR-4 ) 

BVAR I ANCE-( BSUM2/FL0AT ( 2*NC*2*NR-4 ) ) -BMEAN-BMEAN 
BSTDDEV-SGRT  < BVAR I ANCE ) 

BS I GNO I -PAVE / BSTDDEV 


4-5  5 


cc 

c 

c 

c 

c PLOT  MAGNITUDE 

C 

DO  92  I 92- 1 , NC 
DO  92  J92-1.  NR 
SS-CABS< IS< 192.  J92) ) 

IS( 192, J92I-CMPLX1SS. 0 > 

92  CONTINUE 
IFLG-O 
SUN 1-0 
SUM2-0 
NN-0 

TYPE  "SEVEN" 

DO  8 16-1. NR 
DO  8 15-1 , NC 
SS-REAL ( IS< 15.  16) ) 

IFlSS  LE  0 ) SS-0 
IS  < 15.  I61-CMPLX (SS, 0 ) 

C DISCOUNT  0 BACKGROUND  IN  GRAYSCALE  DEFINITION 

IF(SS  GT  0 ) NN-NN+ 1 
SUMl-SUMl+SS 
8 SUM2-SUM2*SS«*2 

C 


C 


■ 

35 

i C 

I C 

f C 


10 


12 

C 

C 


11 

3 

C 


TYPE  "EIGHT" 

RMEAN-SUM1/NN 

RSTD-SQRT  < ( SUM2-SUM1 **2/NN ) / (NN-  1 ) ) 
DELTA-RSTD/2 

CHECK  FOR  MAGNITUDE  TOO  HIGH  FOR  LINEAR  PLOT 

UP-RMEAN+SDEV 

IP ! UP  LT  1023)  GO  TO  35 

I F i CFLG  EG  2)  GO  TO  35 

TYPE"SWITCH  TO  LOO  PLOT" 

GFLC-2 
GO  TO  30 
CONTINUE 

MAKE  A NEW  GRAY  SCALE 

IP (CFLG  EG  1 ) GO  TO  3 
DO  10  1-1, 10 

KTFMP-IFIX( RMEAN+FLOAT l 1-5  ) -DELTA ) 

IF (KTEMP  LT  1)  KTEMP-1 

IF<KTEMP  GT  1023)  KTEMP-1023 

UAL ( I ) -KTEMP 

DO  12  1-1,  10 

UAL ! I >-100  *FLOAT( I > 

IFiGFLO  EG  3>VAL( I > -FLOAT ( 1-1 )*10  +9 
CONTINUE 


J-l 

DO  11  1-1,  1024 
IF(VALIJ)  LT  I)  J-Jfl 
IF(J  GT  10)  J-10 
TAB( I )-J 
CONTINUE 


NNC-NC 
NNR-NR 
IF (GFLG  EG 
FORMAT  ! / / / " 
IF (CFLG  EG 
FORMAT  < ///" 


WRITE! 12. 915)  PKMG, RMEAN. RSTD 


0 OR  CFLG 
LINEAR  PLOT") 

2)  WRITE! 12. 910) 
LOO  PLOT" ) 


EG  1 (WRITE! 12.  905) 


t 

| 


J 


I 


4-56 


913 

900 

l 

C 


20 

c 


13 

C 


23 

27 


441 


FORMAT  ( " MAX--.F6  0."  MEAN-". F6  0.  "SDCV-" . F6  0> 

MR  l TE  i 12. 900>NNH,  NNC  . VAL 
FORMAT (IX. "NROWS-". 14. 

" NCOLS-". I4/1X, "ORAV  SCALE-" . 1013/ / > 

DO  IS  1 1-1.  NR 
DO  20  12-1 . NC 
KTEMP-REALi 1S( 13. II ) ) 

IF (KTEMP  OT  1023)  K TEMP -1023 

IF ( KTEMP  CT  VAL ( 10) ) KTEMP-VAL( 10) 

KBDF (12*10) -KTEMP 

CALL  P ICPLT (KBUF . 260. VAL.  ORAT,  LBUF.  TAB) 

CALL  PLPT (LBUF ( 1 . 1 ) ) 

CALL  PLPT ( LBUF ( 1. 2) ) 

CALL  PLPT (LBUF ( 1 . 3) ) 

CONTINUE 

DO  27  11-1. NR 

CALL  RDBLK (2.  1 1 - 1 . DD  1.1ER' 

IF  ( IER  N£  l)  STOP  —RDBLK  ERR  — 

DO  23  12-1. NC 
IS( 12. II >-DD( 12) 

CONTINUE 

CALL  CLOSE (2.  IER) 

MR  I TE ( 12  441 >PKMO. BMEAN. BSTDDEV. BSIONOI 

FORMAT ( 1H0. 3HPEAK-. 2X. EI4  4/ 1 X.  SHMEAN-, 2X.  E 14  4. 2X.  7MSTDDEV-.  2X 
. E 1 4 4/ 1 X . 4HS  'N-, 2X . E 14  4) 

RETURN 

END 


4-S7 


S/10/78 


r 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

1 

c 


GSPLTX  FR 


REV  A 11/1/78 
REV  B 11/2/78 
REV  C 11/28/78 
REV  0 12/6/78 


SUBROUTINE  OSPLT < IS. NC.  NR.  OFLO > 

OSPLT  MAKES  A OREV  SCALE  PLOT  OF  THE  MAQNITUDE 
OF  THE  NR  X NC  COMPLEX  ARRAY  IS 
QFLO-  0 NEW  LINEAR  ORAY  SCALE 

1 PREVIOUS  ORAY  SCALE 

2 LOO  ORAY  SCALE 

3 NORMALIZE  PEAK  VALUE  TO  99 

4 NORMALIZE  PEAK  VALUE  TO  1000 
COMMON/OO/VAL ( 10  > . ORAT t 30 > 

COMMON/ JJ/LBUF< 132,  3> 

DIMENSION  KBUF t 260) . T AB (1024 ) , AVE ( 64 > 

INTEOER  TAB. VAL.  ORAT,  OFLO 

COMPLEX  IS(64. 64). DD ( 64 ) , CPKMO.  CPNORM.  CSUM 


DATA  LBUF ( 131 . 1)/S. 12K/ 
DATA  LBUF( 131, 2)/S, 12K/ 
DATA  LBUFt 131, 3)/S. 12K/ 
TYPE  "ONE" 


INPUT  NUMBER  OF  POINTS  TO  BE  AVERAOED 
ZERO  AVE  ARRAY 

IF (OFLO  OT  3INAVE-0FL0-3 
00  TO  333 

ACCEPT  "NUMBER  OF  POINTS  AVERAOED" . NAVE 
333  CONTINUE 

DO  202  1202-1.64 
AVE  t 1 202 ) -0 
202  CONTINUE 


TYPE" IN  OSPLT" 

DAI  A ORAT/ iOOK.  100K.  122K.  111K.  122K.  13SK.  4*177K, 

• IOOK.  122K,  IOOK,  144K.  1SSK.3*122K.  135K.  177K. 

• IOOK.  IOOK.  122K.  11  IK,  122K,  1S3K, 4*177K/ 

CALL  0PFN<2, "SCRATCH". 0. IER. S12) 

DO  1 I 1 - 1 . 260 
KBUF t I 1 > -0 


DO  23  I 1 - 1 . NR 
DO  21  12-1.  NC 
21  DD i 12  > -I  St  12.  II > 

CALL  WRBLK  < 2.  1 1- 1 . DD,  1 . IER  > 

IF  t IER  NE  1)  STOP  —WRBLK  ERR  — 

23  CONTINUE 

30  CONTINUE 

TYPE  "TWO" 

PKMG-0 
DO  S J-l. NR 
DO  S I — 1 . NC 
SS-CABS ( IS( I. J) > 

IF ( OFLO  EQ  2)  SS-10  *AL0010(SS-H  ) 
IFISS  OT  PKMO ) PKMO-SS 


j 

1 


i 


f 


4-S8 


I 


rtrtnn  on  non  r>  on  o o> 


CONT INUE 

type  "three- 

NORMALIZE  INPUT  ARRAY 

IP  < GFLO  LT  3)00  TO  87 
PNORM-tOOO 

IFICFLC  QE  3)PN0RM-99 
CPNORM-CMPLX(PNORM.  0 ) 

CPKMC-CMPLX<PNORH/PKNO,  0 > 

DO  88  188-1.  NC 
DO  88  J88-1.  NR 

IS) 188. J88)-IS( 188.  J88)*CPKMC 

FIND  THE  "NAVE-  NUMBER  OF  PEAK  POINTS 

SS-CABS ( IS< 188. J88) ) 

IF ( SS  LE  AVE < NAVE ) ) 00  TO  204 
AVE l NAVE  > -SS 
NAVEL-NAVE- 1 
DO  206  1206-1. NAVEL 

IF (AVE l NAVE* 1-1206)  LE  AVE (NAVE- I 206 )) 00  TO  204 
TEMP -AVE  < NAVE- 1 206  > 

A VE ( NA VE - 1 206 ) -AVE  < NAVE ♦ 1 - 1 206 ) 

AVE ( NAVE* 1 - I 206 ) -TEMP 
206  CONTINUE 
204  CONTINUE 


88  CONTINUE 

TYPE  "FOUR" 


CALCULATE  STATISTICS  OF  BACKGROUND 

87  CONTINUE 

CSUM-CMPLX (0  .0  ) 

B SUM-0 
8SUM2-0 

DO  89  189-1.  NC 
A-CABSl IS( 189, 1 > ) 

B-CABS ( IS ( 189. NR  ) ) 

CSUM-CSUM* IS ( 189, 1 )*IS( 189, NR) 

IS ( 189, 1 )-CPNORM 
IS< 189. NR)-CPNORM 
8SUM2-BSUM2+A*A-»8*8 
99  CONTINUE 

TYPE  "FIVE" 

NROWS-NR- 1 

DO  91  191-2. NROWS 
A-CABS( IS( 1. 191 > ) 

B-CABS) l S i NR,  191  ) ) 

C SUM-C  SUM-*  I S ( 1 . 191  )MS)NC,  191) 

IS) 1.  1 9 1 ) -CPNORM 
IStNC,  191 ) -CPNORM 
B3UM2-BSUM2-*-A*A+B*B 
91  CONTINUE 

TYPE  "SIX" 


4-59 


i PLPT.  SR 


04/23/78 


. TITLE  PLPT 
. ZREL 
. TXTM  1 
. ENT  PLPT 

. EXTD  . CPYL,  FRET,  . STOP 
NREL 

j SUBROUTINE  PLPT (LBUF( 1, J) ) 

i LBUF  IS  A 132.  WORD  BUFFER 

LBUF— 167 
1 

PLPT:  JSR  t.  CPYL 

LDA  3. LBUF, 3 
STA  3,  TEMP 
LDA  0, PASS 
MOV  0.  0,  SZR 
JMP  PASS2 
ISZ  PASS 
SUBC  0. 0 
STA  O,  N 
JMP  . +2 

LP:  ISZ  N 

LDA  0, ABUF 

LDA  2,  N ; CHANNEL  # 

SYSTM 
. CHSTS  77 
JMP  . +1 
LDA  0,  BUF 
LDA  l.LN 
SUB  0,  1, SNR 
JMP  PASS2 
DSZ  NN 
JMP  LP 

i — ERROR,  THE  *LPT  HAS  NOT  YET  BEEN  OPENED 

; (IT  MUST  BE  WRITTEN  TO  BEFORE  CALLINC  PLPT) 


PASS2: 


K132D: 
TEMP: 
CNT : 
PASS: 
N: 

NN: 

ABUF: 

LN: 

APBUF: 

BUF: 

PBUF: 


CHANNEL  # 


- 4LPT  NOT  YET  OPENED/ 


i OLD  BUFFER  POINTER 
; ADDRESS  OF  PLOT  BUFFER 


JSR  « STOP 
. TXT  /(PLPT)  - 

LDA  1 , K132D 
MOVZR  1, 1 
STA  l.CNT 
LDA  3, TEMP 
LDA  2, APBUF 
LDA  O,  0,  3 
MO VS  O, 0 
INC  3,  3 
LDA  1.0,3 
ADD  1 , O 
STA  0,  0,  2 
INC  2, 2 
INC  3, 3 
DSZ  CNT 
JMP  NXT 

LDA  0, APBUF 
MOVZL  0,  O 
LDA  1.K132D 
LDA  2,  N 
8YSTM 
. WRS  77 
JMP  . +1 
JSR  «. FRET 


132 

0 

0 

0 

0 

62. 

BUF 

“ 4*400+ "L 

PBUF  I ADDRESS  OF  PLOT  BUFFER 
. BLK  22 

. BLK  67.  i -132/2  + 1 EXTRA 

. END 


I PLOT  BUFFER  ADDRESS 
i BYTE  POINTER 

,'  GET  CHANNEL  NUMBER 


os/oa 


MASTER  FILE 

DIRECTORY 

HOLO 

03/27/79  12 

•LPT. 

0 

RAP  03/27/79 

12 

01 

03/27/79 

10000003 

• TTI 

0 

APM  03/27/79 

12 

01 

03/27/79 
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APPENDIX  A 

FEATURE  EXTRACTION  BY  HOLOGRAPHY 


We  are  presented  with  the  problem  where  we  wish  to  use  a hologram  to  distin- 
guish between  two  similar  reference  patterns,  Pa(x)  and  Pb(x).  Because  the 
patterns  are  similar,  we  may  say  they  are  partially  correlated.  We  may  also 
express  the  similarity  by  saying  the  two  patterns  share  a common  subset. 

The  measure  of  similarity  of  the  two  patterns  is  given  by  the  correlation 
coefficient: 


"ab 


,P  (x)p.  (x)  dx 

v a b 

y lpa‘(x)dxiPb‘(y')dy 


(A-  1 ) 


In  terms  of  the  Fourier  transformed  patterns  we  may  write 


~ab 


R '(?)  R.(?)  d£ 

♦ a d 

V'!|r  |'d?  * R,l 'd? 
V . a • b 


(A- 2) 


We  recognize  that  R l‘  and  1 R.  I"  are  constants,  by  condition  of  the 
experiment . 


The  numerator  of  Equation  A-2  can  also  be  recognized  as  the  integral  over  a 
constant.  The  argument  for  this  is  as  follows:  Reconstruction  of  a hologram 
made  by  reference  p.t(x>  depends  on  the  similarity  of  reconstruction  reference 
pb( x'l  . The  similarity  nay  be  denoted  by  the  correlation  between  these  two 
functions.  If  we  treat  the  functions  as  vectors  in  a space  spanned  by  the 
independent  variable,  x.  we  mav  represent  the  process  as  in  Figure  A-l: 


FIGURE  A-l.  VECTOR  REPRESENTATION  OF  CORRELATION 


A-l 


We  see  that  P^,  may  he  resolved  Into  two  vectors.  One  vector  lies  along  the 
original  reference  and  is  the  correlated  part  of  p^.  We  can  arbitrarily  define 
this  correlated  vector  as  the  subset  of  which  is  common  to  both  Pa  and  P^. 
The  other  vector  which  is  a component  of  p^  is  orthogonal  to  pa-  Thus,  in 
the  reconstruction  process  we  have  a term  like: 


- ^coitmon^a  * °or thogona l^a^ 
■ (constant)  x 7 + (noise 


(A-3) 


We  get  the  reconstruction  because  oconrnort^)a  must  meet  the  reconst ruct ion 
condition  of  being  delta  like.  This  means  that  the  Fourier  transform  of  this 
term  must  be  a constant.  We  can  therefore  write: 

*b^a  * (^om^  + Perth  ® P J RCom  Ra  + Rorth  Ra  (A-4) 


Thus,  ■ -e  may  write 

Vr  Vr 

R R,  " R * constant.  (A-S) 

abba 


Since  we  have  integrals  over  constants,  the  area  ot  the  recording  cancels  out 
from  the  numerator  and  denominator  of  Kquation  A-2.  Thus,  Kqu.U  ion  A-2  is 

equr.ui.eiu  to 


C 


ab 


R 


R 

a 


■it 

R.  R 
b a 

R R, 

a b 


correlation  coefficient  between  a + b 


(A-h) 


The  correlation  between  the  two  similar  patterns  f and  ~ may,  by  an  appro- 
priate twist  of  Equation  A-h,  be  put  to  use. 


Suppose  we  wish  to  form  a multiplexed  hologram  such  that  advent  of  the  refer- 
ence pattern  c a will  cause  a unique  reconstruction  of  the  associated  object 
pattern  ja.  Similarly,  will  reconstruct  only 


An  ordinary  multiplexed  hologram  will  lead  to  a false  reconstruction.  For 
example,  the  hologram 


H 


R 

a 


S 

a 


+ 


(A- 7) 


will  produce  the  output 


R R,  S. 
a b b 


R H 
a 


(A-S) 


when  Che  reference  Pa  is  incident.  Equation  A-8  may  be  interpreted,  with  the 
aid  of  Equation  A-b  as: 


R H 
a 


IrJ  IrJc 


ab 


(A-9) 


We  see  that  the  unwanted  output  is  reconstructed  in  proportion  to  the  degree 

of  correlation  between  P and  p.  . 

a b 

If,  however,  we  make  a hologram  which  includes  subtractive  terms  it  becomes 
possible  to  cancel  the  false  reconstruction. 


Such  a hologram  would  take  the  following  form: 


* * * 

R S + R.  S.  - k.R,  S 
a a d b 1 b a 


- k.R  *S. 
2 a b 


(A- 10) 


In  Equation  A-10  and  k^  are  constants.  We  evaluate  these  constants  by  per- 
forming a reconstruction.  Thus,  with  P^  incident  we  get  the  output 

- (IrJ2  - V.Os.  + <»A*-  k2lR.|2>Sb  <A-11) 

If  ki  is  adjusted  so  that  there  is  no  output  of  the  false  reconstruction,  S^, 
it  will  take  the  value: 


R K*  IR.I 

k ■ — " ■ ■ C , , bv  Equation  (A-6). 

- t R I2  !r  1 "b 


(A— 12) 


a 


a 


By  a symmetrical  experiment  we  find  that 
a 


W*  IrJ 

k,  - — - — — C . 

1 ! ^ i 2 i _ i ab 


R, 


V l“b! 


(A- l 3) 


gives  no  output  of  when  reference  is  incident. 


A- 3 


By  proper  adjustment  of  kj  and  k>  the  hologram  may  be  tuned  to  give  unique 
outputa  when  either  of  the  two  reference  patterns  are  Introduced.  There  is, 
however  a penalty:  The  desired  reconstruction  is  somewhat  reduced.  To  see 
the  magnitude  of  this  effect  we  introduce  Equations  A-12  and  A-l.l  into  the 
reconstruction  Equation  A— 11.  The  result  is 


(A-U) 


R H 

a 


Similarly, 

V - V’1' 


fA-lM 


iA-loi 


In  the  treatment  of  the  correlation  as  a geometrical  process,  the  correlation 
coefficient  is  simply  the  cosine  of  the  angle  between  the  abstract  vectors 
c4(x>  and  We  see  from  Equation  A-15  that  the  reconstruction  is  pro- 

portional to  the  square  of  the  sine  between  these  two  vectors. 

That  is,  it  the  reference  functions  are  orthogonal,  the  reconstruct  ion  will  be 
maximal.  Another  wav  of  interpreting  this  is  to  say  that  only  the  orthogonal 
(uncorrelatedl  portions  participate  as  references  in  the  reconstruction. 

The  process  outlined  above  is  analogous  to  the  Oram- Schmidt  orthogonal isat ion 
procedure  familiar  from  mathematics.  However,  there  is  the  substantial  differ- 
ence that  the  orthogonallzatlon  is  virtual.  The  reconstructed  patterns  do  not 
need  to  be  orthogonal  and  no  physical  rotations  of  the  reference  patterns  take 
place. 

The  ability  of  the  hologram  to  perform  orthogona Illations  on  sets  of  patterns 
shows  that  its  power  in  pattern  recognition  is  substantially  greater  than 
previously  suspected. 


APPENDIX  B 


EXISTENCE  PKOOK 


Our  purpose  here  is  to  show  chat  a multiplexed  hologram  may  be  formed  from  two 
separate  subject-reference  pairs  in  such  a wav  that  reconstruction  of  either 
subject  will  be  free  of  crosstalk  and  other  clutter.  In  particular,  it  will 
be  shown  that  this  is  possible  even  when  the  two  reference  patterns  are  highly 
correlated . 


Vie  start  the  proof  by  establishing  two  reference  patterns  Rj(.\l  and  K . 

Each  of  these  patterns  are  broken  up  into  two  subsets.  One  of  the  subsets  will 
be  used  to  reconstruct  the  associated  subject,  and  the  other  will  cancel  out 
the  crosstalk  from  the  other  hologram  of  the  ensemble.  Thus,  we  write: 


K.  - 


S * R1C 


R, 


- R' 


r;  + r,. 


ib-1) 


Here  we  have  designated  R{  as  the  subject  reconstructor  and  as  the  subset 
which  cancels  the  effects  of  partial  correlations  between  Rf  and  Rj.  Similar 
definitions  hold  for  the  division  ot  R>. 

In  forming  the  hologram  we  pair  Rj  with  a subject  Sf  and  Rj  with  Si.  We  assume 
adaptation  according  to  feedback  principles  discussed  in  the  main  text.  Thus 
the  terms  in  the  multiplexed  hologram  will  reach  an  equilibrium  state  resem- 
bling that  of  Equation  B-3.  We  may  write  the  multiplexed  hologram  as  the 
summat ion: 


s s2 
h - K;  + R, 

1 R2 

In  order  for  all  crosstalk  to  be  cancelled,  the  conditions 

V • S1 


and 


R,H  - S , 


i,B-  n 


ib-S'l 


must  be  satisfied. 

Then,  using  the  definition  R-l  in  Equation  R- 1,  we  get: 


r;s 

S1  + R ’ + Rlc” 


By  constraining  this  result  according  to  Equation  B-s  we  must  have 


r;s 

~rT  + RicH  ' 0 


By  a similar  argument  we  derive: 

r;s, 


+ R-  h * o 


R 


It  we  eliminate  H between  Equations  B-7  and  B-8,  we  find 

-i! 

> 

K’,‘ 


Rlcfl 

R,S, 


or : 


r’ 

K,  S, 

l 

Ic  l 

R: 

>/  RA-S2 

tB-71 


IB-81 


IB-10> 


Next  we  substitute  Equation  B-10  into  B-2  to  get: 


H - 


r: 


Factoring  terms,  we  have: 


H 


l 


K! 


1 V 


tb-l  n 


iB-U’> 


By  a symmetrical  derivation  or  by  employing  Equations  B-10  in  B- 1 , we  can  show 
that  H also  takes  the  form 


Wo  can  test  this  result  by  attempting  a reconstruction.  Playing  the  reference 
of  Kquation  B-l  through  the  hologram  of  B-12,  we  get 


let  us  examine  the  term  R,  H.  Bv  equation  B-7  we  have: 

Ic 


k;s 

R1CH  * * R* 


(.B-lsl 


V l 


tB-lnl 


U'-l 


Inserting  B-10  into  B-l'  we  have 


R,  \\  - - 
1 c 


\cSl$2 

K,  , 


tB-lSl 


Thus  we  see,  according  to  R- IS  that  the  last  tvv'  terms  in  B-lt>  cancel . This 
leaves  the  desired  result: 


KjH 


•Sl 


tB-l^ 


Bv  a similar  argument  we  can  also  show  that: 


K ,M 


S , 


tB-:oi 


What  we  have  done  here  is  to  demonstrate  that  there  are  no  pathologies  required 
to  form  a hologram  which  will  result  in  the  reconstruct  tons  ot  Kquat ions  B-l^ 
and  B-JO.  An  important  point  to  note  is  that  this  result  is  only  possible  it 
the  hologram  is  formed  from  ratios  of  subjects  and  references.  This  implies 
a new  type  of  ort hogona l itat ion  process  is  operative  here.  It  is  not  clear, 
at  this  time,  that  the  results  can  be  generalis'd  to  multiplex  ensembles  con- 
sisting of  three  or  more  suhholograms. 


B-  l 


t 


There  are  Interesting  symmetries  that  come  about  when  it  is  assumed  that  the 
cancellation  subsets  of  the  two  references  are  in  fact  the  mutually  correlated 
(or  common)  portions  of  these  two  patterns.  Then  we  have 


lc 


2 c 


R 

c 


(B-21 ) 


Then  the  first  simplification  comes  in  Equation  B-10,  which  produces,  after 
cancellation  of  terms: 


We  get  the  second  simplification  by  eliminating  Rt  between  Equations  B-7  and 
B-8.  The  simultaneous  solution  of  these  two  equations  leads  to 


R!  (R.  R.,  1 r - S.S.)  - 0 
1 lc  2c  12 


( B—  23) 


If  Rj  is  to  be  nonzero,  we  must  have 


/V  2 

>/  R1cR2c 


(B-24) 


In  the  case  where  R equals  R.,  , Equation  B-24  takes  the  form 


•JV2 


(B-25) 


In  the  very  special  case  where  S 


1 


S,  ■ S we  have: 


(B-26) 


Thus,  the  basic  ratio  form  is  preserved  throughout  all  these  processes. 


B-4 


